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1.0 EXECUTIVE SUMMARY 

Various aspects of the Shuttle Orbiter S-band network 
communication system, the S-band payload communication system, and the 
Ku-band communication system are Included In this final report. 

The first topic. Section 3.0, deals with a proposed method of 
obtaining more accurate S-band antenna patterns of the actual Shuttle 
Orbiter vehicle during flight. The preliminary antenna patterns using 
mock-ups are not realistic since they do not Include the effects of addi- 
tional appendages such as wings and tall structures. 

Section 4.0 discusses the Ku-band communication system, espe- 
cially the TORS antenna pointing accuracy with respect to the Orbiter 
and the modifications required and resulting performance characteristics 
of the convolutional 1y encoded high data rate return link to maintain 
bit synchronizer lock on the ground. Also Included Is the TORS user 
constraints on data bit clock jitter and data asymmetry on unbalanced 
QPSK with noisy phase references. 

Section 5.0 outlines the S-band payload communication system 
study, including the advantages and experimental results of a peak regu- 
lator design built and evaluated by Axiomatix for the bent-pipe link 
versus the existing RMS-type regulator. Also discussed is the selection 
of the nominal sweep rate for the deep-space transponder of 250 Hz/s, 
justifying this conclusion with a phase-plane analysis which Incorporates 
the transponder operating conditions and tolerances. And finally, the 
effects of phase noise on the performance of a communication system Is 
analyzed in great detail. The phase noise spectra is characterized using 
experimental results and included into models of the carrier tracking 
loops to determine the phase noise system errors of the overall system. 
Different types of two-way links are considered, including the phase 
noise effects of multiple internal control loops which are equivalent 
to three- and four-way links. Both in-band and out-of-band noise sources 
are also included. 

Eight appendices are included at the end of this final report 
which describe in detail the models used, the derivations and the con- 
clusions of the studies for those who are interested in the technical 
aspects of the investigations. 
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2.0 INTRODUCTION 

This final report summarizes the system Implementation studies 
of the S'Band Network Communication System, the Ku-Band Communication 
System, and the S*Band Payload Communication System. The specific areas 
of study are covered under Contracts NAS 9-15240E and NAS 9-15240F; how- 
ever, In the interests of continuity of conciseness, the results of the 
study will be reported as a unit. 
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3.0 S-BAND NETWORK COMMUNICATION ANTENNA CALIBRATION STUDY 

This section presents the results of the Investigation to date 
of an approach to calibrate the Orbiter S-band antennas (and UHF antenna 
also) during STS flights. The actual gain patterns of the flush-mounted 
S-band quadratures (PM, two on each side of the Orbiter) and the S«band 
hemi antennas (FM) are difficult, at best, to predict from the customary 
ground antenna range pattern measurements. This Is because, for flush- 
mounted antennas, the Orbiter structure Itself has a significant Influ- 
ence on the pattern. Consequently, the adjacent sections of the Orbiter 
structure In which the antennas are mounted have been mocked-up and, 
along with the antenna(s) mounted In them, used for pattern tests on the 
JSC antenna range. This technique typically yields fairly accurate gain 
measurements near the boresight of the antenna, I.e., normal to the sur- 
face of the structure. However, as the angle off boresight Increases, 
the accuracy In the gain pattern decreases due to edge effects of the 
limited structural mock-up and the missing influence of the more remote 
parts of the structure, such as the wings. 

Even though the operation of these antennas is expected to take 
place primarily in the central portion of the main beam, it is important 
to have accurate knowledge of the natural antenna gain pattern throughout 
the entire main beam and the si delobes. This is the case for the follow- 
ing reasons: 

(1) During the ascent and descent portions of flight, the 
range to the ground stations will be relatively short. As a result, the 
signal strength will be quite high. The configuration of the Orbiter and 
the ground station, however, may be such that communication will be via 
the low gain portion of the mean beam or even a sidelobe. 

During lift-off, the present geometric configuration is 
such that a weak portion of the S-band pattern will be in use. As a 
result, more precise information about the gain patterns is necessary In 
these regions. At such short ranges, even in the low gain portions of 
the pattern, the SNR should be very adequate, but the extent of such 
coverage is not really known. 
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(2) The long-term plan Is to replace the S-band quadratures 
with two-position switched beam phase-array antennas mounted In the same 
physical location. The S-band quadratures are presently planned for the 
first six to nine STS flights. The main problem of appendage blockage 
is reduced by the use of phased-array antennas. This would not be nec- 
essary, however. If by chance the S-band quadratures operated satisfac- 
torily in the fringe areas (i.e. , In the vicinity of -1 to +1 dB). 

(3) The long-term plan is also to transfer communication to 
the Shuttle/TDRSS links via the S-band quadratures or S-band phased 
arrays as soon as possible Into the STS mission. Because of the predic- 
ted narrow design margins for some of the S-band Shuttle/TDRSS links, it 
is important to show how large the spatial region is in which acceptable 
link performance can be expected. This can be determined with accurate 
knowledge of the S-band antenna patterns in the "fringe areas." A prime 
opportunity to accomplish this finer calibration exists during the STS 
phase because, during STS, communication with the Shuttle will Le via 
STON stations. While in orbit, the actual performance of the Shuttle 
antennas will be free of all the ground effects present on conventional 
antenna ranges. Thu':, with appropriate link calibration and data pro- 
cessing, the STS Orbiter/STDN link can function in the same manner as 

an antenna range. 

The remainder of this section describes a technical approach 
to be used during STS to obtain the finer gain calibration. The results 
described are preliminary, and further refinement of the approach is 
expected to continue on a subsequent contract. The emphasis is placed 
on the main proolem, namely, the changes in the S-band quadrature pat- 
terns due to appendages. 
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3.1 General Technical Approach 

The RF link Information which will be available during STS for 
purposes of Shuttle antenna calibration will fall Into two general cate> 
gorles. First, there Is the link Information which will * i obtained dur- 
ing communication while the Orbiter Is In an "as flown" mission. The other 
category of link Information Is that which Is obtained while the Orbiter's 
attitude Is determined and controlled specifically to support antenna cal- 
ibration. The FTR (Flight Test Requirenent) for this category Is to 
determine antenna coverage that cannot be verified on the ground. The 
only objective In these "special antenna-pattern tests" Is to utilize 
specific antenna coordinate profiles, or "passes," to obtain gain informa- 
tion for that portion of the antenna(s) pattern. This would be Implemented 
by having the Orbiter perform suitable roll maneuvers for two or three 
ground-station passes. The present plan is to perform a 720" roll at 2“/ 
sec, requiring approximately six minutes. During such a maneuver, only 
one of the S-band quadrature antennas would be operating, thereby provid- 
ing low gain and sidelobe information. The 720° roll would be performed 
twice, once for one of the lower quadrant S-band antennas and once for one 
of the upper quadrant antennas. 

The type of measurement information obtained from the programmed 
attitude profile will be superior to the "as flown" category. This Is 
because the portion of the pattern used in the "as flown" flights will 
necessarily concentrate on the main beam of one of the S-band quadratures. 
The antenna in use will be chosen so that the best signal strength is 
obtained. Tne flights programmed for antenna measurements will provide 
broader pattern coverage and will come closest to conventional antenna 
pattern measurement techniques; however, the "as flown' measurements will 
be the predominantly available Information due to the pressure of numer- 
ous mission requirements for STS flights. The signal processing and data 
reduction developed for the antenna measurements will be the same, however, 
for both categories. 

The general approach to updating the calibration of the Orbiter 
S-band antennas from the measurements described above is as follows; 

(1) Calibrate the STDN stations and Orbiter RF systems 
The best information we have at this time is: 
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(a) The Shuttle transmitter power can be determined to 
within ±0.1 dB. 

(b) NASA Goddard presently plans to accurately calibrate 
one of the ground stations. At this time, it Is not known whether this 
calibration will produce values that are known to within ±1 dB. 

(2) Record uplink and downlink received signal strength versus 
Orbiter structure. 

(3) JSC will provide BET (Best Estimate Trajectories) and atti- 
tudes of the flights. These will be sufficiently accurate that they can 
be assumed to be deterministic. 

(4) The data will be filtered to provide estimates of the actual 
antenna gain in the immediate vicinity of the location where the measure- 
ments were taken. The estimates can be updated as more STS flights are 
taken. For example, data from STS-1 may be replaced with actual measured 
points from STS-2 or used in addition to those f om STS-2 and succeeding 
flights. Part of the filtering process will consist of a "data good" 
algorithm to determine when data should be accepted or rejected. 

(5) In those regions where the orbital antenna measurements were 
not made and where the antenna pattern is relatively smooth, an interpola- 
tion algorithm will be developed to provide updated values of the antenna 
gain. For the sidelobes and/or in those regions where the antenna gain 
variation is very large with respect to either jlevation or azimuth, lit- 
tle if any interpolation is likely. In these regions, the average antenna 
pattern gain may be adjustable (biased) from that of the ground measure- 
ments. The overall antenna gain measurement process is enhanced because 
the gain of the main beam near boresight is expected to have no appreciable 
change due to the appendages. 

All aspects of the filtering algorithm will be fundamental in 
determining actual antenna coverage capability that could not be verified 
on the ground. All aspects of the filtering and interpolation process 
will also continuously compare measured orbital values with those obtained 
from the ground test-range patterns, particularly when confidence in the 
orbital measurements is low. 
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Figure 1 Illustrates the basic mathematical link relationship for 
determination of antenna gain from received signal power. From this fig- 
ure, It Is seen that. In addition to Orbiter attitude, the key parameters 
which must be known or measured are received signal power, calibrated STbN 
gain and circuit losses, Orbiter transmitted power and circuit losses, 
llne-of-sight range, and carrier frequency. This calculation Is not nor- 
mally encountered In conventional antenna measurement procedures because. 

In conventional procedures, a standard gain horn Is substituted for the 
antenna being measured and the received signal power Is measured. There- 
after, the dB difference In received signal power with the measured antenna 
Indicates the gain relative to the standard gain horn. 

Ideally, a completely calibrated STDN consists of precise knowl- 
edge of all gains In the power budget computation. Including antenna feed 
circuit losses, preamp gain and receiver gains. The signal strength is 
determined by monitoring the A6C levels. The relationship between A6C 
level and signal strength must therefore also be determined (calibrated). 

In addition, since it may be desirable to measure received signal power 
on the Orbiter, the STDN gains and losses and power in the transmit mode 
would ideally be accurately known. Similarly, the Orbiter transmit power 
and circuit losses In the transmit and receive modes would ideally be 
accurately measured. 

However, an error budget analysis will be carried out on all of 
the parameters necessary to normalize the raw data (AGC output) into a 
raw estimate of antenna gain. 

Parameter Bias RMS Error 

Ran-'e 
Attitude 
Circuit Losses 
Frequency 
Transmitter Power 
STDN Gain 
Gain Uncertainty 
AGC Quality 
AGC Nonlinearity Curve 
(dBm versus volts) 

Both bias and RMS error need to be estimated, as Indicated in the 
above table. With these estimates, the errors of th“ raw artenna gain 
data can be estimated. The range data is expected to be very accurate, 
based on Best Estimate Trajectory (BET). 



NOTE; All loss < 0 dB 


S-Band 

Antenna 



Downlink Communications Antenna 
Gain Determination 


Uplink Communications Antenna Gain 
Determination 


^RSTDN " '■^1 * *-R‘ SsiDN"^ ^^2 

“ ^RSTON” *^TU" '-^1 ■ '■R' SsTDN' ^’^2- 

where 

^TU = Shuttle transmit antenna gain 
Gru = Shuttle receive antenna gain 
Pjy = Shuttle transmitter output power 
Ppy = Shuttle received power at receiver 
LT^ = Shuttle transmit circuit loss 
LR^ = Shuttle receive circuit loss 

Lj^ = path loss 


’’rU * '^TSTDN * ‘-^2 * ®TSTDN * ‘’R’*’ ^RU " ‘■*^1 
Su * **RU ■ **TSTDN " *-^2 " ®TSTDN ‘ ‘’R " *-*^1 

^RSTON ~ transmit antenna gain 
GtSTON “ GTDN transmit antenna gain 
Pf^STON ” GTDN receive power 
Ptston “ GTDN transmitter power 
LR 2 * STDN receive circuit loss 
LT 2 = STDN transmit circuit loss ' 


Figure 1. Link Relationships for Determination of Shuttle Antenna Gain During OFT 
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Both uplink and downlink data will be available during the 
flights. The weighting of the uplink and downlink data needs to be 
determined. The SNR In each link can be estimated via the previous 
table for each link. One simplification of the overall data reduction 
process Is obtained by consolidating each uplink and downlink data point 
Into one data point. 

The STON stations, however, have heretofore had abundant design 
margins so that sucn accuracies have not been necessary. Therefore, most 
parameters In the power budget computations at most STDN stations are not 
accurately known. The most economical approach to circumvent this Is 
proposed as follows: 

(1) Accurately calibrate one STDN station and use this station 
to make absolute gain measuranents over a portion of the S-band antenna 
patterns. A high degree of confidence would then exist for the absolute 
gain in these portions of the patterns. 

(2) Use the information in (1) above to calibrate the STON sta- 
tions when the Orbiter passes by. Measurenents over the remainder of the 
patterns would then provide antenna gain values relative to those which 
are known absolutely. 

The best portion of the pattern in which to make these absolute 
measurements is the region of boresioht. This has the additional advan- 
tage that the pattern in the main beam should closely agree with the 
measurements made in the antenna ground test range. It is expected that 
the appendages will not appreciably affect the antenna pattern in the 
high-gain portion of the main beam. 

An alternate method of obtaining the absolute antenna gain mea- 
surements is by way of a calibrated standard gain horn at one STDN site. 
Whether or not this approach involves fewer calibration difficulties is 
yet to be determined. 



10 


3.2 Data Gathering Procedure 

The basic Information gathering and processing procedure to obtain 
the gain pattern Is shown In Figure 2. Provisions for gain calculations 
from both uplink and downlink communications are shown. The majority of 
Figure 2 depicts the 1nfom»t1on flow necessary to evaluate the equations 
shown In Figure 1. The box labeled ^^Unprocessed Gain Matrix" represents 
the three-dimensional matrix of measured antenna gain samples. This matrix 
Is J>;p1cted In Figure 3. The X and Y coordinates correspond to the quan- 
tized e.t look angle coordinates of the Orbiter antenna. The Z coordin- 
ate corresponds to the measured gain samples. There may be more than one 
sample for a given coordinate, resulting from the several passes by the 
sanie station or samples of the gain at that coordinate for several stations. 
In general, the elements of the three-dimensional gain matrix are Identi- 
fied by 

where 

G|^ * gain sample which Is a function of station, Orbiter antenna, 

station elevation angle, Orbiter lock angle, and Orbiter trans- 
mit circuit 

N » number of gain samples for the particular cell 

The matrix of GN„ . will. In general, have many cells for which there 
will be no values and others for which there will be several values. Thus, 
to arrive at the final antenna pattern. It Is necessary to pass this matrix 
through an appropriate filter/estimator. 



VEHICLE GROUND STATION DATA 



Figure 2. Information Flow for Gain Pattern Determination of Shuttle Antenna During STS 
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The drawing In Figure 3 should not lead to the conclusion that 
all sampies are separated by ae and In e and respectively. The 

actual locations of the measurements are not now presently known. It can 
be safely assumed that the measurement locations will be Irregularly 
spaced. The present plan for the "special antenna-pattern tests" Is to 
sample the AGC at 20 samples/sec during the roll maneuvers described above. 
This corresponds to 10 samples/degree • which Is quite adequate. 

After the signal processing has been completed, the desired output 
Is a grid of estimated antenna gain values. A satisfactory grid size 
appears to be 2" x 2". This Is certainly adequate over the main beam of 
the S-band quadratures. However, the size of the computational problem 
needs to be determined. The array of values for each of the quadrature 
antennas Is 90 x 45 based on the 2"-by-2" grid. This Is also assuming 
that each antenna pattern Is evaluated only over the region where It Is 
the primary antenna. This evaluation will be performed during the next 
phase of the Investigation. 
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3.3 Processing Algoritt»i> Considerations 

The signal processing for Shuttle antenna measurements can be 
partitioned Into two parts* as follows: 

(1) Preprocessing. The effects of range, frequency, circuit 
losses, biases, etc. are taken Into account. In addition, the necessary 
bookkeeping, such as which antenna Is In use, look angle, data source, 
etc.. Is carried out. The result Is a set of normalized antenna mea- 
surements which are used to update the appropriate antenna pattern. 

(2) The algorithm. The computational procedure employed to 
update the antenna gain pattern using the normalized antenna measurements 
described above. 

A recursive Ka1man-Bucy type estimation algorithm Is discussed 
below as a candidate suggested by Prof. B. Tranter, which has received 
extensive attention in the recent literature on signal processing. Is 
the Maximum Entropy approach. This method has been shown to have excel- 
lent performance In the estimation of power spectral densities when only 
sparse data is available. The limitation appears to be that highly 
reliable data is necessary or, equivalently, that the observations are 
at very high signal -to-noise ratios (SNR). This signal -processing tech- 
nique w.il be examined to see if it can be extended to spatial filtering 
at intermediate values of SNR. In this case, the dominant noise is the 
calibration error variance. 

If these techniques do not prove to be adequate, the class of the 
nearest neighbor algorithms will also be examined in the next phase of 
the investigation. 

We cannot determine which of these algorithms is optimum since 
essentially no statistical information is available. The real issue is 
which of the candidate algorithms will do the most satisfactory job. 
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3.3.1 Summary of Kalman-Bucy Til taring 

In this section, the basic Ideas of the discrete time Kalman-Bucy 
filter are presented. Restriction Is made to discrete time processing 
since that Is most amenable to digital computation and the raw data Is 
also expected to come to us In discrete time. 

The Kalman-Bucy filtering algorithm determines the best linear 
mean-square estimate of a signal Immersed In additive noise when the sig- 
nal and noise can be modeled In a certain manner. This basic signal model 
Is quite general and can be made to fit a wide variety of signal statistics. 

3.3.2 Basic Signal Model 

The basic signal model presented Is not the most general but Is 
more than adequate for the antenna calibration problem. The basic signal 
model 1s diagrammed In Figure 4 and Is described analytically by the fol- 
lowing equations: 


\ <» 

\ ( 2 ) 


wher. 


X. = a stochastic signal vector which is modeled in (1) as a 
first-order Markov process 

N|^ * measurement noise 

W|^ * input noise process which makes the signal stochastic 

Yk * Sfc + Nk is the measurement or observed process 

Sk = H'|( Xj^ is the system input process, a weighted value of the 
desired signal, X|(. It is the observed signal when there is 
no measurement noise. The notation' means transpose. 

k represents a discrete time argument, k 0, and without loss 
of generality, k takes on integer values only. 

We shall make the following assumptions with regard to the noise 

processes: 

1* {N|^} and {Wj^} are individual white processes. By white in 

discrete time, we mean that, for any k and a, k and are inde- 

pendently random vectors and Wj^ and are independent random vectors. 
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2. (N|^} and {W|^} are Individual zero mean, Gaussian random 
vector processes with known covariance matrices, namely, 

E ["k "»•] - “k «kt »•> 

E [«k «.•] = «k f’*’' 

where is the Kronecker delta. 

3. {N|^} and {Wj^} are independent random vector processes. 

In equations (1) and (2), {Fj^}, {Gj^}, and {H|^} are sequences of 
deterministic (generally known) matrices which are chosen so that the 
first-order and second-order statistics of {X|^} and {Yj^} agree as closely 
as possible with those of the real world system being modeled. 

3.3.3 Basic Kalman-Bucy Estimation Algorithms 

Given the basic signal model in equations (1) and (2) and in Fig- 
ure 4, the Kalman filter detennines the best mean-square estimate of the 

process {Xj^}, given the sequence of observations {Yj^}. The basic discrete 
time Kalman filtering problem can equivalently be stated as follows: 
Suppose that the initial state vector, Xq, of the desired 
signal is a Gaussian random vector with mean, E(Xg) = Xq 
and covariance matrix Pq, and independent of the noise pro- 
cesses {Nj^} and {Wj^}. 

Then determine: 


Vk-1 -E[*kl’', ■ 

\/k ^ E . * i k j 


h/k “ ^ ^ 1 J 

and the associated error covariance matrices, and respec- 
tively. In (4a), the best estimate is denoted by and represents 

a one-step predictor in that the estimate of {Xj^} is determined for time 
k given all of the observations up to time k-1, namely, {Yq, Y^ ,. . . ,Yj^_^}, 
In (4b), the estimate is and represents a filter estimate in that 
{X|^> is determined for time k given all observations up to and including 
time k, namely, {Yg, Y^,..., Yj^}. 
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We distinguish between the one-step predictor and the filter 
versions of tfe Kalman algorithm since the one-step predictor Is signifi- 
cantly simpler than the filtering algorithm. At this stage of the antenna 
calibration p'oblem development. It Is enough to consider only the one- 
step predictor version of the Kalman algorithm. 

Referring again to (4), the notation on the right side represents 
the conditloial expectations (ensemble average) of what Is to be estimated 
given what has been observed. Under very general conditions, this Is the 
best mean-square estimate. 

Witnout proof, we now present the one-step predictor version of 
the class of Kalman estimators. 

The Kalman one-step predictor comprises the system depicted In 
Figure 5 and is described for k ^ 0 by the equations 

\+l/k ' l^k " \/k-l * \ ^k 

with 

V -1 ^0 

The gain matrix Kj^ is determined from the error covariance matrix 
by 

" '"k ^k/k-1 \ [^k'^k/k-1 “k 


(5) 

( 6 ) 


assuming the inverse exist? (which is normally the case), 
error covariance matrix associated with the estimate 


The conditional 
is defined as 


L 


k/k-l 


+ E 


( 


\ " ^k/k-1 



1 < k-1 


(3) 


and is minimized by the estimate defined in (4a). This conditional error 
covariance matrix can be computed recursively by the discrete time Riccati 
equation, namely, 
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^k+l/k ' ^k^k/k-1 "^k/k-1 \^”k'^c/k-l \ * *^k^ 
“k'^k/k-1 j’^k' ^ ®k ^k ®k' 


(9) 


This equation Is Initialized by 


^0/-l ’ **C 


3.3.4 Properties and Advantages of the Kalman Filter 


Me hereby enumerate certain properties and advantages of the 
Kalman algorithm: 


1. The Kalman filter Is a linear, finite-dimensional, discrete- 
time system. This Is most advantageous In that it could have been nonlin- 
ear and/or infinite dimensional. 


2. The output of the system is dependent, of course, 

on the input The conditional error covariance matrix, on the other 

hand. Is Independent of the input sequence {Yj^}. As a result, the error 
covariance matrix sequence } and the gain matrix sequence {K|^} can 

be computed before the filter Is actually run. 


3. The conditional covariance matrix definition in (8) identi- 
fies as the covariance matrix associated with a particular esti- 
mate. Since is Independent of {Yj^}, however, we may take the 

expectation of both sides of (8) over all possible {Y|^} to conclude that 



( 10 ) 


This means that is also an unconditional error covariance matrix 
associated with the Kalman filter, i *6* is the covariance matrix 
associated with the estimator. 


4. The Kalman estimator is redrawn in Figure 6 to further empha- 
size its structure as a copy of the original signal model driven by the 
estimation error. The Kalman estimator can also be interpreted as develop- 
ing an estimate of the observation at time k given all observations up to 

^ A 

time k-1, namely as shown in Figure 6. The error (Y|^- -j). 
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is tile input to the gain matrix K|^ and is, essentially, what drives the 
estimator. If the error is zero, the copy of the basic signal model pres- 
ent in the estimator runs by itself and continues to generate additional 
estimates. 

5. Suppose the underlying signal model is time invariant and 
the input and output noise processes are stationary. Then, F|^, G|^, H|^, 

Q|^, and are all constant (independent of k). In general, however, 
^k/k-1 therefore K|^ will not be constant, so tiie Kalman estimator 
will normally still be time-varying despite time invariance and station- 
arity in the signal model. 

3.3.5 Application of Optimal Linear Filtering to the 

Antenna Calibration Problem 

In this final section, we present the first version of tailoring 
the Kalman-Bucy estimation algoritlvn to tiie antenna gain calibration prob- 
lem. As updated information is obtained about the accuracy of the various 
ground station and Shuttle parameters, the filtering algorithm will also 
be appropriately updated. We begin by making various observations: 

1. For the present, we restrict attention to the S-band quads. 

The des'red output is an updated antenna pattern which has a grid size of 
2®-by-2 . That is, the value for each 2®-by-2® component of the entire 
pattern will be an estimate of the average gain in that region or the 
gain at the center. There is no way to distinguish between which of the 
two values are being represented by the desired output. 

The specification of the angles for antenna pattern specifi- 
cation is shown in Figure 7. Based on a grid size of 2®-by-2®, there are 
4050 values necessary to specify the pattern for all four S-band quads. 

The antenna pattern will be specified by the vector at 
time k. Therefore, is a 4050-dimensional vector. This is too large. 

We reduce it by a factor of four by considering S-band quad antennas sep- 
arately. The dimensionality of X|^ is therefore reduced to 1012, which is 
still quite large. The dimensionality can be further reduced substantially 
by considering only those positions of the patterns which have a reasonable 
chance of being used. At the present time, it is not clear how much 
reduction in dimensionality is appropriate. 
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For each of the fcnir S-band quads, there are three Shuttle 
configurations [1]: 

a. Solid Rocket Bcrasters (Si®) External Tanks (ET) -*■ 
Orbiter. (ITils configuration lasts from liftoff to the first approximately 
120 seconds.) 

b. ET 0rt>1ter (first approximately 500 seconds) 

c. Ortiter only. 

Ideally, therefore, there will be 12 distinct antenna pat- 
terns and each data point which Is collected applies to only <me of the 
patterns. 

2. The raw data that Is collected are AGC voltage measurements 
onboard the Orbiter and In the various ground stations, as shown In Fig- 
ure 8. This data must be transformed Into average received power at the 
various locations. As shown In Figure 8, signal strength will be simul- 
taneously available from the Orbiter receive** and one of the ground sta- 
tions. To perform this transformation requires calibration of the AGC on 
the Orbiter and at the various ground stations. At present, we assume 
that this transformation Is carried out external to the filtering progrmn 
so that the inputs to the antenna calibration program consist of signal 
strength. 

3.3.6 Variable and Parameter Assignments for the Antenna Model 

The basic signal model In Figure 4 is tailored to represent the 
antenna pattern as a function of tl'^e. As indicated above, X|^ represents 
the vector of components which constitute the antenna pattern for one of 
the S-band quads In one of the configurations. To be specific, we assume 
that the Orbiter Is in orbit since most of the data generated will be for 
that configuration. We make the following assignments: 

Xq represents the antenna pattern measured In the antenna ground 
range and is assumed known 

X|^ , k ^ 1 represents the antenna pattern for the antenna mounted 
onboard the Orbiter at the kth tine period. More will be said 
subsequently about how k represents time. 



AGC Vol tage 
Data Onboard 
Shuttle 


Transformation to RCVD 
Average Power 
Signal Strength 


Average Received Power 
Onboard Shuttle 
Vs Time 


AGC Voltage 
Data at 

Transformation to RCVO 

Average Received Power 
at mtfi Ground Station 

Ground Station 

Average Power 

Vs Time 


Signal Strength at mt^ 



Ground Station 



Figure 8. Conversion of AGC Voltage to Received Average Power Signal Strength 
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Each Integer value of k represents that portion of one flight 
»^en coonunicatlon Is wlt^ one ground station. This Is best seen by 
Inspection of Figure 9. where an example Is presented. lYie time k * 1, 
for example, corresponds to the Orbiter conmunlcatlng with ground station 
BDA (see [1]) for designations, locations, etc.) during which measurements 
were taken at those antenna component Indices across from the vertical bar 
above k * 1 In Figure 9. 

The antenna pattern associated with is the modified version of 
Xq due to the appendages on the Shuttle and due to being mounted on the 
Orbiter itself. What we are after is an estimate of X^. It Is assumed 
that the antenna pattern remains constant throughout a given mission. 
Therefore, 


X, = X 2 = ... = = ••• = (11) 

vdiere represents the last time period of the first mission. The 
antenna pattern X.| is related to Xq by 

X, = F(,X„ (12) 

where the matrix Fq is unknown. Of course, our problem can be equivalently 
staff'd as finding Fq. We assume that Fq is diagonal, unknown and constant. 
The values of F|j are 

F|j = 0 , 1 1 k < - 1 (13) 

Equivalently stated, the matrix F|^ renains zero throughout the remainder 
of the first mission. 

Since the antenna pattern is assumed to remain constant, then 
G,^ = 0 , Q|c = 0 ®11 (14) 

At the beginning of the second mission, we assume that the antenna 
pattern is uniformly attenuated by the factor "a." This is modeled by 


( 15 ) 
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which Is equivalent to saying 


F,. = a I (16) 

where I Is an Identity matrix. The attenuation factor "a” Is Initially 
unknown and results from the TPS tiles being charred during reentry. The 
attenuation parameter "a" can be estimated* however, from signal strength 
measurement variations from one mission to the next. 


3.3.7 Variable and Parameter Assignments for the Observations 

The weighting between the value of antenna gain and the signal 
strength observation Is performed by the matrix Hj^. This matrix Is also 
diagonal. For a particular value of k. It has nonzero values only for 
those Indices which correspond to antenna positions where measurements 
were obtained, as Indicated In Figure 9. 

The nonzero values of Hj^ convert the antenna gain to signal 
strength. This Is accomplished by the communication link equation 


Pt ^T S X 


(17) 


For the downlink, the parameters In (17) are as follows: 

P-j- = Orbiter transmitted power 

Gj = antenna gain of Orbiter, at some designated angular coordin- 
ates (e, (|)) 

6j^ = Ground station antenna gain 
X = wavelength 

R = range between Orbiter and designated ground station. 

All measurements taken at a given ground station associated with 
a particular orbital pass will be associated with a particular value of k. 
For the downlink, each component of H|^ will contain 
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which relates the antenna gain onboard the Orbiter to the signal strength 
In the ground station receiver. 

3.3.8 Basic Flow Diagram for Antenna Calibration Estimation 

A basic flow diagram for antenna calibration estimation Is shown 
In Figure 10. The Kalman estimation algorithm Is shown In Figure 6» where 
a11 necessary matrices are assumed known. In Figure 10 and the previous 
sections, we have described the Inputs required to compute the matrices 
In Figure 6. 

As more specific Information Is obtained regarding the availabil- 
ity and accuracy of the necessary parameters, the basic flow diagram will 
be developed In further detail. 



Time 
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Figure 10. Basic Flow Diagram for Antenna Calibration Estimation 
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3.4 STDN Calibration Considerations 

In order to assess the magnitude of the problem of accurately 
calibrating the gain of the STDN stations, a preliminary Investigation 
of the techniques utilized by DSN to calibrate DSN stations at S-band has 
been undertaken. The Investigation has found that the DSN has the capa- 
bility of calibrating station antenna gain to a 3a accuracy of ±0.1 dB or 
better. This accuracy Is achieved as a result of making a series of very 
exact measuranents and utilizing averaging techniques on the data. On a 
one-shot basis, wherein one or very few measurements are used, a 3o accu- 
racy of ±0.3 dB can be expected. 

The technique utilized by DSN essentially utilizes accurately 
measured radio stars as RF sources. The flux density of several selected 
stars, listed In Table 1, was determined by means of an accurately cali- 
brated DSN station, DSS-13 (Goldstone 26M). This station was calibrated 
by means of a gain calibration transfer technique. The ultimate source of 
calibrated antenna gain was an S-band Standard Gain Horn developed and 
accurately calibrated by JPL specifically for this purpose. The major 
steps for calibrating DSN stations, via the radio star, beginning with 
the Standard Gain Horn, are shown In Figure 11. A block diagram of the 
DSN gain transfer system Is shown In Figure 12. The complexity and prob- 
lems associated with Implementing these types of calibration techniques at 
selected STDN stations will be investigated In the next phase of the study. 


Table 1. Absolute Flux Density and Antenna Temperatures at 2295 MHz 


Source 

Flux Density 
S ± la, JY 

100% Efficient Antenna 

26 m 

34 m 

64 m 

3C123 

31.0 

± 0.61 

5.92 

10.2 

36.2 

3C218 (Hydra A) 

26.7 

± 0.51 

5.10 

8.78 

31.2 

3C274 (Virgo A) 

136 

± 2.6 

25.9 

44.6 

158 

3C405 (Cygnus A) 

887 

± 15 

169 

292 

1034 

3C461 (Cassiopeia A)* 

1525 

± 27 

291 

502 

1777 


* 

Flux density known to decrease approximately 1% per year (observation 
epoch 1972.6) 











Figure 11. Flow Chart for Calibrating Arbitrary DSN Station 
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4.0 KU-BAND COMMUNICATION SYSTEM STUDY 

The Ku-band communication system study Investigated the 
accuracy with which the TORS Ku-band antenna could be, pointed at the 
Orbiter Ku-band antenna, as discussed in Section 4.1. The Orbiter 
Ku-band system high data rate convolutionally coded return link was 
modified to eliminate the need for a minimum transition density to 
keep the bit synchronizer on the ground in lock. Section 4.2 pre- 
sents the performance characteristics of the modification to the 
convolutional encoder. The performance effects due to TDRSS user 
constraints of data bit clock jitter and data asymmetry on unbalanced 
QPSK with noisy phase reference are discussed in Section 4.3. 
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4.1 Analysis of TDRS Ku^Band Pointing Accuracy 

In order to meet EIRP specifications, the TDRS must point Its 
high-gain Ku-band antenna to within 0.24** of the user. TRU has done 
extensive analysis and simulation of the open-loop pointing budget and 
has concluded [6] that the TDRS Is capable of meeting the requirements 
under worst-case conditions. This report Is a review of the TRW analysis, 
with emphasis on the Space Shuttle as a user. Certain parameters consid- 
ered by TRW Improve with respect to the Shuttle, and others degrade. The 
net effect, if the TRW analysis were flawless, would be that the TDRSS 
could easily service the Shuttle with the specified pointing accuracy. 
However, TRW assumptions as to residual errors after calibration appear 
to be optimistic. As a result, the ability of the TDRS to point to the 
Shuttle to within 0.24° Is marginal with more realistic assumptions as to 
residual callbratirn errors. 

In subsequent sections, we discuss the various parameters used to 
estimate the pointing accuracy, with emphasis on the more critical param- 
eters. Error budgets are given using our best estimates of attainable 
measurement errors. We consider errors with the Shuttle at 100-m11e and 
200-m11e orbits and the TRW worst-case atmospheric explorer orbit. 

4.1.1 TDRS Calibration Errors 

Item 17 of the attached pointing budgets appears to be the most 
controvc 'sial in the pointing budget. On orbit, the TDRS high-gain Ku- 
band antennas will be calibrated using the White Sands terminal. Specif- 
ically, the autotrack null will be calibrated using the autotrack error 
signal [7]. Prior to calibration, the autotrack null to mechanical bore- 
sight uncertainty Is greater than 1°. 

During calibration, the antenna Is swept past the White Sands ter- 
minal and the autotrack error signal is used to estimate the position of 
the null relative to the gimbal angle readout. Errors introduced by lags 
in the system are first calibrated out by offsetting the null and sweeping 
a fixed pattern. Local maxima and minima are correlated with gimbal angle 
readouts to estimate the system time lag. 

Table 2 of [6] Is a list of the errors Involved in the calibration 
procedure. Due to the complicated nature of the model, the resultant 
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residual calibration error has been determined by TRW using a Nonte Carlo 
sinRilation. The fixed error after calibration, itaii 17 of Table 3, is 
estinmted at 0.0418^ N-S and 0.0471” E>M. The autotrack noise (item 
of Table 2) used in the simulation is based on the worst-case minimum 
autotrack scale factor of 0.041 V/V/”. He feel tiiat It is probable that 
the null can be measured to tte accuracy claimed by TRH since worst-case 
estimates of the error terms were used in the model and the slope of the 
autotrack null is steep. 

TRM has used the same technique to estimate the residual calibra- 
tion errors of the transmit peak (boresight). The intent is to calibrate 
out the null-to-transmit boresight misalignment. This misalignment is 
composed of two parts: an electrical boresight term and a boresight null 
shift term. The electrical boresight term is specified at 0.06® ar.d the 
null shift term is specified at 0.032®, for a possible worst-case mis- 
alignment of 0.092®. TRM, however, has chosen to optimistically use an 
electrical misalignment of 0.013® and a null shift misalignment of 0.012®. 
Inasmuch as the contractor building the antenna will sign up to only the 
specified values (0.06* and 0.032®), we feel that the more optimistic 
assimiptions by TRW are not valid; consequently, we use the specification 
values in the estimates of the residual errors- In addition, these are 
not random errors— they are fixed quantities. Hence, they are added 
directly, net RSS'ed, to the worst-case E-M residual error of the auto- 
track null. This would not pose a problem if the transmit boresight could 
be calibrated; the misalignment error could be removed. 

The transmit boresight misalignment is calibrated by observing the 

output of a Hewlett-Packard Model 436A digital power meter, as the beam 

is swept past the ground station. The peak is tagged and correlated with 

the gind>al position readout. TRW has calculated that, at a received C/N 

of 45.5 dB, the peak can be read to an accuracy of 1.22x10"^ dB, equiva- 
-4* 

lent to 8.94x10 . Two problems arise: the power meter has a specified 

accuracy of ±0.02 dB and a least count (resolution) of 0.01 dB, and the 
fast variations (item 19 of Table 3) have a peak of 0.0876®, equivalent 
to a variation of 1.1 dB. Either item will preclude the claimed accuracy. 
The fast errors, consisting of quantization noise and mechanical noise, 
will appear in the power meter. Even in the absence of the fast varia- 
tions, the time constant of the power meter will prevent the measurement 
to even 0.01 dB. 
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(Max) 
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0.U66 

4.2 

0 . 0 & 

0.05 
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Variations 

(deg) 


0.026 


U.U5UU 

O.OlOU 

O.OISO 


l« Software 
Antenna/Booa 


2. Thermal (Boom) 

3. Thermal (Antenna) 

4. Elect/Nech Alignment Spacecraft Z 
to Misalignment 

Spacecraft Z to Z^ 

Measurement Uncer^lnty 

Z,^ (Mech) to A/T Elect 

Eartu Sensor 

5. Null Loc Alignment 

6. Mech Alignment 

7. litas Errors 

8. Power Supply 

9. Radiance Uncertainty 

10. Noise and Quantization 

Conti ol System 

11. Step Size 1 Dynamic Oscillations 

12. CDA Position Error 

13. Drive Axis Wobble 

14. Yaw Coupling Errors 

15. User Epheineris 

16. TDRS Position Determination 

nr r .r r nw . • /N-S 0.0418, 

17. Fueo Error After Calibration ojoy 

18. RSS of the Slow Varying Errors 
’9. i’?S of the Fast Varying Errors 

20. uncertainty Ellipse Major Axis (3o) 0.170 

21. Uncertainty Ellipse Minor Axis (3<>) 0.103 


22. Prooability of Pointing within 
U.24 degree 

23. Pointing Loss to tSA Forward 
ElKP (dB) at 0.24 degree 


Exceeds 0.U3 degrees for 20 sec following slew start or stop of other SA antenna. 
**ExcLeds O.Ub degrees for 15 sec following slew start or stop of other SA antenna. 


0.03 

(0.04)* 


0.06 

(0.08)** 

0.04 


0.0382 

0.113 


0.014 0.014 


0.0550 0.0836 0.120 


0.0876 



Tables. TDRS/Shuttle Transmit Pointing Budget, 0.06** Null-to-Boresight Error 
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In tlie next section, we discuss revised pointing budgets based on 
null-to-boresight misalignments of 0.06** and 0.092**. The 0.06** misalign- 
ment Is used as a likely value to be encountered with the actual antenna, 
and the 0.092** Is the worst-case specification. These are added directly 
to the null E-W calibration error to form the mean of the probability 
ellipse used In the calculation of the probability pointing within 
0.24** of the Shuttle. 

4.1.2 Open-Loop Pointing Budget 

Tables 3 and 4 are the open-loop pointing budgets, TDRS-to-Shuttle. 
Items 1-6 are taken directly from [6]. Items 1, 2 and 3 are beyond the 
scope of this analysis and are taken at face value. The values given are 
not unreasonable; hoi^ver, verification would require extensive analysis 
of the software algoritivns ana the mechanical model of the system. Items 
4, 5 and 6 are calibrated out by the procedure described earlier, with 
residual errors given in item 17. Itons 7, 8 and 9 remain unchanged in 
the N-S direction but, for the Shuttle, increase to 0.0494**, 0.0099** and 
0.0148**, respectively, in the E-W direction. They are given as the pitch 
error allocation (0.05**, 0.01** and C.015**) times cos (target elevation), 
with target elevation equal to 8.9** for the Shuttle. 

Items 10-13 are also taken directly from [6]. Again, verification 
of these quantities would require extensive analysis of the system hard- 
ware and structural design. 

Item 14, yaw coupling, is decreased from worst-case 0.129** for the 
atmospheric explorer to 0.0382** for the Shuttle. This error is given as 
0.25** sin (target elevation). 

Item 15, user ephemeris, is given as the 9-second uncertainty 
times the relative user angular velocity, 0 . With h the altitude of the 
user in miles, 

^ 1.)436x 10^^ ^ (h + 3959) . 

(6.371 X 10® + hx 1.609X 10^)^^^ (22284 - h) 

Thus, for h = 100 miles, e = 1 .25 x 10” ‘’/see and the uncertainty 
is 0.113**. For h = 200 miles, 0 = 1 .24 x 10”^**/sec and the uncertainty is 
0.112**. We use the worst-case value of 0.113** In the revised pointing 
budgets. 
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Item 16, TORS position determination. Is taken directly from [6]. 
The N-S error of Iten 17 (0.0418®), fixed error after calibration. Is 
taken from [6] and Is a result of TRW simulation of the effect of all 
anticipated error sources, as discussed In a prior paragraph. The E-W 
residual error Is the direct sum of the transmit peak-toHfiull boresight 
misalignment and the residual null calibration error (0.0471®). 

Items 18 and 19 are the RSS of the error In the respective col- 
imns. Item 20 Is the RSS of Items 18 E-H, 18 yaw and 19. Item 21 Is the 
RSS of Item 18 N-S and Item 19. Items 20 and 21 are termed the peak 
values of the errors and used as 3o values for probability calculations, 
with Item 17 N-S used as the mean of the minor axis and Item 17 E-W used 
as the mean of the major axis of a bivariate Gaussian distribution. 

Item 22 Is the Integral of the bivariate Gaussian over a circle 
of 0.24®, which gives a pointing loss of 9.6 dB. 

4.1.3 Results of Revised Pointing Budget An^^.lysis 

Two cases of null-to-boresight misalignment were analyzed for the 
Shuttle orbit. A nominal 0.06® misalignment is used to represent a com- 
promise between the optimistic 0.03® used by TRW and the worst-case 0.092® 
in the specification. The results of this analysis are presented in 
Table 3. The 0.06® case is equivalent to a pointing loss of 0.5 dB during 
calibration of the transmit boresight. This is well above the accuracy 
of the power neter, and may be attainable if not masked by the fast varia- 
tions extant in the hardware. The probability of pointing to within 0.24® 
(item 22) is 0.986 for this case. 

As a worst-case, a null-to-boresight misalignment of 0.092® is 
used. From Table 4, the probability of pointing to within 0.24® is 0.954. 

For purposes of comparison with the TRW worst-case (atmospheric 
explorer) budget. Table 5 gives the pointing budget based on the 0.137® 

E-W worst-case residual error. TRW calculates a probability of pointing 
within 0.24® of 0.993, while the worst-case estimate gives a probability 
of 0.898. 
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0ne«01nensional Errors 


Slow Variations (3<r) 


(deg) 


1. Software 
Antenna/Boom 




Fast 

Variations 

(deg) 


0.026 


0.013 


.0494 - 
.0099 - 


0.0150 0.0148 - 


2. Thermal (Boom) 

3« Thermal (Antenna) 

4 . Elect/Mech Alignment Spacecraft Z 
to Misalignment 

Spacecraft Z to Z^ 

Measurement Uncer^inty • 

Z^ (Mech) to A/T Elect 

Earth Sensor 

5. Null Loc Alignment 

6. Mech Alignment 

7. Bias Errors 

8. Power Supply 

9. Radiance Uncertainty 

10. Noise and Quantization 

i 

Control System | 

11. Step Size & Dynamic Oscillations 

12. GUA Position Error I 

13* Drive Axis Wobble 

14. Yaw Coupling Errors 

15. User Epheiueris 

16* TORS Position Determination 

nr .r riK C N~S 0.0418 

17. Fued Error After Calibration ^ o.l37 

18. RSS of the Slow Varying Errors 

19. RSS of the Fast Varying Errors 

20. Uncertainty Ellipse Major Axis (3o) 0.178 

21. Uncertainty Ellipse Minor Axis (3o) 0.103 


22. Probability of Pointing within | 
U.24 degree 

23. Pointing Loss to KSA Forward 
tlRP (dli) at 0.24 degree 


Exceeds 0.U3 degrees for 20 sec following slew start or stop of other SA antenna. 
Exceeds U.Ub degrees for 15 sec following slew start or stop of other SA antenna. 


0.03 

(0.04)* 


0.06 

(0.08)** 

0.04 


0.0382 

0.113 


0.014 0.014 


0.0550 0*0836 0.120 


0.0876 



Table 4. TORS/Shuttle Transmit Pointing Budget, 0.092“ NuU-to-Boresight Error 
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l« Software 
Antenna/Booai 

2* Thermal (Boor^ - 0.021 0.013 

3. Thermal (Antenna) 0.0 0.06 

4. £1ect/Hech Alignment Spacecraft Z 

to Z^ Misalignment 1.0 . 

Spacecraft Z to Z^ 

Measurement Uncer^inty 0.(166 

ZjA (Mech) to A/T Elect 0.2 

Earth Sensor 

5. Null Loc Alignment 0.05 

6. Mech Alignment 0.05 

7. Bias Errors 0.U500 0.0429 . 

8. Power Supply O.OlOO 0.0086 - 

9. Radiance Uncertainty 0.0150 0.0129 

10. Noise and Quantization 

Control Systoii 

11. Step Size & Dynamic Oscillations 

12. QUA Position Error 

13. Drive Axis Wobble 

14. Yaw Coupling Errors - - 0.129 

15* User Ephemeris - - 0.122 

16. TDKS Position Determination 0.014 0.014 

17. Fixed Error After Calibration o!l37^ 

18. RSS of the Slow Varying Errors 0.0550 0,795 0,1776 

19. kSS of the Fast Varying Errors 

20. Uncertainty Ellipse Major Axis (3o) 0.213 

21. Uncertainty Ellipse Minor Axis (3<>) 0.103 


22. Probability of Pointing within 
0.24 degree 

23. Pointing Loss to KSA Forward 
ElkP (dB) at 0.24 degree 


Exceeds 0.03 degrees for 20 sec following slew start or stop of other SA antenna. 
**Exceeds 0.0b degrees for 15 sec following slew start or stop of other SA antenna. 



0.129 

0.122 


0.014 0.014 


0.0550 0,795 0,1776 



Table 5. TDRS/Atmospheric Explorer Pointing Budget, Worst-Case 
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4.1.4 Conclusions 

A significant loophole has been found in the TRU analysis of the 
transmit pointing budget. While the pointing loss may be marginally 
acceptable for the Shuttle, other users could experience difficulty in 
acquisition. These results are predicated on the other errors being 
within budget. It should be emphasized that the remaining error budget 
items are TRW engineering estimates rather than hard specifications. 
Errors greater than the TRW estimates due to meclianical noise, quantiza- 
tion effects, thermal effects or larger than anticipated Earth sensor 
errors will adversely affect the ability of the Shuttle to acquire TORS. 

Should the null-to-boresight error, however, be larger than the 
worst-case 0.092**, it could probably be calibrated to this value using 
the techniques proposed by TRW. This error is equivalent to an error of 
1.25 dB at the power meter, and should be discernible if the other errors 
are within bounds. 
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4.2 Perfonnance Characteristics of the Ku*Band Channel 3 

Convolutional Code with G 2 Inversion 

The present implementation of the high data rate (mode 1. 
channel 3) error-correcting code can engender system performance degra- 
dation under certain input conditions. The data-derived clock at the 
receiver requires a minimum transition density, e.g., a transition-free 
string of data exceeding the maximum specified value of 64 s^bols can 
cause loss of lock at the bit synchronizer. Uncoded input data consist- 
ing of all zeros (not necessarily an unlikely event) will result in a 
coded output of all zeros and loss of bit synchronizer lock. In order 
to preclude this event, the code is being modified by inversion of the 
second encoded check bit ( 62 ) out of the encoder. This check bit is 
reinverted after demodulation and prior to decoding. 

Users should be made aware of two minor problems with the 
proposed method of 62 inversion. The encoded symbol stream can slightly 
exceed the maximum transition-free length under certain pathological (low 
probability) conditions and, with high probability, the decoder will out- 
put meaningless data if the uncoded input data is all zeros or all ones. 
The first problem is probably of no concern since the input data giving 
rise to tha long transition-free run of coded symbols has a low prob- 
abilicy of occurrence. The second problem may be of some concern since 
input data consisting of all zeros may be valid if the user ignores the 
constraints of at least 64 transitions in 512 bits and no more than 64 
bits without a transition. 



4.2.1 Run Length Performance of the Convolutional Code 

with Gg Inversion 

Figure 13 is a representation of the return link encoder with 
62 inversion. This code is derived from the NASA Planetary Standard 
rate one-half constraint length 7 code shown in Figure 14 by replacing 
g(x) with g(x^). The resultant code is an interleaved version of five 
constraint length 7 codes and is decoded with five parallel Viterbi 
decoders. The run length properties of the K = 7 code have been inves- 
tigated by Simon and Smith [8] who have shown that an input sequence of 
10101001011001 (read left to right) will result in an output sequence of 
1000000000000001 with G2 inversion. That is, the maximum run length is 
14. The results of Simon and Smith can easily be extended to the modi- 
fied K = 31 code by repeating each input symbol five times. Referring 
to Figure 1, by alternately concatenating five check symbols from the 
first n.jd 2 adder and five from the second, each symbol in the original 
K = 7 code will be replaced by five synbols. Thus, if each input bit is 
repeated five times, each check symbol will be repeated five times and 
the maximum run length is extended to 70 symbols. Note, however, that a 
unique 6?-bit input data stream (or its complement) is required. Assum- 
ing “random" input data at 50 Mbps, this event should occur once every 
several thousand years. 



INPUT 

UNCOOCO 

CHANNEL 


2 MBPS TO 
50 MBPS 
NRZ-L»M.S 


3i>STAGE SHIFT REGISTER 



NOOULO-2 ADDERS: NO. 2 


SYMBOL INVERTER 


5-STAGE SHIFT REGISTERS: 


ADDER CONNECTION COEFFICIENTS 


1. I0204I0000I 

2. l0004llOOOAl 


/ f 

f 9 9 i i 
9 9 i t t 

III 1 1 1 1 irn— 

tO-BIT OUTPUT BUFFER 
LOADED AT EVERY FIFTH 
UNCODED INPUT BIT 


OUTPUT 

ENCODED 

CHANNEL 


k HOPS TO 
100 HOPS 
NRZ-L 


4 ^ 

cn 


Figure 13. Rate One-Half Encoder Representation, Ku-Band Return Link, 


Mode 1 



INPUT 



Figure 14. The Rate One-Half, Constraint Length 7 NASA Planetary Standard Code 
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4.2.2 Code Behavior with DC Input 

With all zeros (or all ones) as Input to the encoder* the 
resultant coded sequence will consist of alternating strings of five 
zeros and five ones, e.g., ...111110000011111..., which Is the vdiole 
point of Implementing G 2 Inversion. In the event that the five decoders 
are not synchronized to the correct phase of the check bit pairs by prior 
transmission of nonzero data, each decoder will arbitrarily get alternat- 
ing ones and zeros as Input. After reinversion, a decoder will attempt 
to decode either all zeros or all ones. Since this Is a transparent 
code, the all ones "encoded" sequence will produce a valid output data 
stream of all ones. This output sequence will be Interleaved with the 
output of the other decoders, producing a cyclic string of ones and 
zeros. This Is because a decoder cannot acquire branch synchronization 
on all ones or all zeros. In fact, under these conditions, the decoder 
win output erroneous "data" 90% of the time. This can be Illustrated by 
examining Figure 15. The relative origin time of the alternating sequence 
can take on one of 10 values with respect to the first of the five 
decoders. In this example, decoder 1 (D^ ) takes as Input check bit 2, 
number 3 inverted at the transmitter and check bit 1, number 3 as the 
receive. Thus, decoder 1 sees a string of all ones, a valid code word, 
and outputs all ones. By contrast, decoder 4 sees all zeros and outputs 
all zeros. The only case wherein all zeros will be output is wnen 
decoder 1 starts decoding check bit 1, number 1 and check bit 2, number 1, 
reinverted at the decoder. Assuming that the decoder starts up at a 
random point on the encoder sequence, this will happen 10% of the time. 

One of the remaining nine incorrect decoded sequences is the all-one 
sequence, an inverted version of the actual data. This, too, presumably 
could be acceptable to the user by virtue of the inverted frame synchron- 
ization. It should be pointed out that periodic frame sychronization bits 
in the all -zero sequence will appear as noise to & decoder which is not 
correctly branch synchronized. On the average, each decoder will see 
five "noise" symbols— not enough to effect branch synchronization. 

The question arises as to the number, if any, of additional 
sequences which have the property that branch synchronization is ambig- 
uous. This can be analyzed formally by examining the response of one 
constraint length 7 encoder. The bit inversion can be ignored since this 




Input Data: 00000 

Rncoded Symbols (with inversion): 000001 1 1 1 1000001 1 

Decoded Data: 1110011 10011100 


Figure 15. Effect of Zero Inp'it Data on Decoder Output 
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has the effect of complementing the encoded symbol sequence if the phase 
of the check bits Is slipped by one symbol and the complement of a code 
word Is also a code word. Let C^(x) and C 2 (x) be the polynomial repre- 
sentation of the encoder output, e.g., 

C-|(x) - C'jQ + ^ ^12^ *** ~ 9i(^) 1(x)* 

~ ^20 ^ ^ 21 ^ ^22 ^ ~ 

with g^(x) the generator polynomial for the first check symbol and I(x) 
the polynomial representation of the Information sequence. If the branch 
synchronization Is ambiguous, the symbol pairs ^20 Sr Si S2* S 2 Ss*-*- 
will form a code word. This can be written in polynomial form as: 

So ^ S2^ ■*■... — 9i(x) I*(x), 

Cii + - 92 (x) I*(x), 

with I*(x) the alternate information sequence. Thus, 

92(x) I(x) = g^(x) I*(x), and 
9|(x) I(x) 

I*(x). 

y 

This pair of equations can be solved to determine the input conditions 
under which branch synchronization is ambiguous; e.g., if €^€2 represents 
a code word, so does C 2 C^. 

Solving for I(x) and I*(x), we find that 



I*(x)jxg2^(x) + g^^(x) = 0. 


and 
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For the NASA constraint length 7 code* 

g^(x) » 1 + + X® + X®, 

and 

g^^(x) * 1 + X + x^ + x^ + x^. 

Thus* 

xg2^(x) + g^^(x) * x^^ + x^^ + + x^ + x® + X® + x^ + x^ + 1 

= (x+l)(x^^ + x^ + X® + x^ + x^ + x^ + 1) 

* P^(x)P2(x) 

factored Into Irreducible polynomials P|(x) and P 2 (x). 

If 

I(x) = Iq ^ + ... + 1 ] 3 X^^ + ..., 

the criterion that I(x)P|(x)p 2 (x) = 0 implies that 

^0 ^ ^ ^3 ^ ^6 ^ " ^8 ^ ^9 ^ ^10 ^ ^12 * ^13 = ® 

or 

^13 ^0 * ^3 ^ *6 ^ *7 ^ ^8 ^9 ^ ho ^ hz ' 

That is, given 13 specific input bits as initial conditions, if the four- 
teenth bit ( 1 ^ 3 ) satisfies the above equation, the output sequence is 
ambiguous. Note, however, in order for this sequence to sustain itself, 
each subsequent input bit must also satisfy the relationship. That is. 


h4 ^ h h ^ h ^ ^7 ^ *8 ^9 ^ ho * hi * hs* 
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The sequence can be modeled the output of a linear feedback shift 

register, as shown in Figure 16, with taps appropriate to P^(x)P 2 (x). The 

properties of these sequences are well known. P^(x) = 1 ■<- x corresponds 

to the all-one sequence and its complement, the all-zero sequence, which 

we know are ambiguous. From [9], me find that p^Cx) is primitive, and 

^ 12 

the shift register model of the sequence will have period 2 - 1. 

Depending on the initial conditions of the register, the sequence gener- 
ated will be all zeros, all ones, or one of two complementary periodic 

12 

sequences, each with period 2 - 1. The probability of any data source 

generating this specific “infinite" sequence with period 4095 is presum- 
ably rather low; hence, we have shown that branch synchronization with 
this code will not be a problem with real data. In any case, if 
I(x)p^(x)P 2 (x) = 0, then I*(x)p^(x)p 2 (x) = 0, and the decoded I*(x) will 
merely by a oelayed verison of I(x). This would certainly not be a prob- 
lan with a single decoder; the only poblem is in conjunction with the 
five interleaved decoders, which is shown to be a result of pathological 
input data. 

The onl> remaining caveat is that the input data to the en;oder 
should not have a period of five, as this would result in each 'h tne 
five decoders seeing a constant string of ones or zeros. 



(a) Source Model for Ambiguous 
Check Symbol Sequence 
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4.2.3 Sunmary of Performance 

With certain input data, the mode 1, channel 3, encoder/decoder 
with 62 inversion will not perform as specified. These conditions are 
delineated below. 

• With very low probability, the encoded symbol stream 
will have transition-free runs longer than 64 symbols. 

• With high probability, the decoder will not attain branch 
synchronization if the user constraints are ignored and 
all zeros are input. 

• With high probability, the decoder will not attain branch 
synchronization if the input data is cyclic with a period 
of five bits. 

• No aperiodic data stream exists which will preclude 
branch sychronization. 

Problem area 1, the possible output of transition-free runs of 
encoded s.nibols greater than 64, is a result of unique input data and 
s!x)uld not be of concern. Problem area 2, the Inability of the decoder 
to synch 'onizer on all -zero data, will be of concern only if a user 
attempts to initiate transmission with data consisting of all zeros (or 
all ones) in violation of the user constraints or if the input data is 
cyclic with a five-bit period. 



55 


4.3 Effects Due To Specific TDRSS User Constraint Parameters 

In the previous study (NAS 9*15240E [10]), all but two Ku-band 
TDRSS user constraints were evaluated and compared to Hughes Aircraft 
Company specifications. The cumulative performance degradation effects 
of the TDRSS user constraints were analyzed and summarized. In this sec- 
tion, the performance degradation of the remaining two user constraints, 
data bit clock jitter and data asymmetry on unbalanced QPSK, are 
presented. 

4.3.1 Ku-6and Bit Synchronizer BER Degradation 

due to Data Bit Clock Jitter 

Appendix I determines the timing error variance and resulting 
bit error rate degradation in a Ku-band Shuttle-compatible second-order 
bit synchronizer ^«ie to input data clock jitter modeled as either a sin- 
usoidal signal or a spectrally flat random process. Both NRZ and Man- 
chester data are considered. 

Results are obtained for the bit synchronizer tracking error 
for both the sinusoidal and random jitter cases. Next, the bit error 
probability for both sinusoidal and Gaussian data clock jitter is derived 
and plotted. Using the timing variance, it is then possible to compute 
the bit error probability and the associated degradation from the jitter- 
free case. 

Finally, the three bit synchronizer types used with Ku-band 
data are evaluated to determine their respective bit error rate degra- 
dation. Depending on which specification, modulation format and syn- 
chronizer are considered, bit error rate degradations vary from less 
than 0.1 dB to over 10 dB. The degradations tend to be small (0.1 dB) 
when the frequency deviation and deviation rate are specified at 0.10% 
bu quite large when specified at 0.1% R^, where R^ Is the symbol rate. 

To evaluate the BER degradation, it is necessary to know the 
bit synchronizer loop bandwidths. Based on a telephone conversation with 
John Roach of Harris Corporation (Melbourne, Florida), the following bit 
synchronizer specifications were obtained: 

1. Ultra-High Data Rate Bit Synchronizer (UHDR) 

R^ = 75-150 Msps 

f^ s 60 kHz (B^ = 200 Hz) 
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Manufacturer: Motorola 

Jitter problems comments: Yes* since they could not make 

the loop bandwidth narrow enough 
due to Internal delays 

2. High Rate Bit Synchronizer (HDR) 

Rj = 10-75 Msps (will operate down to 4 Msps) 
a 0,08% R^ (B^ = 0.267% R^) 

Manufacturer: Harris Corporation 
Jitter problem: No. 

3. Low-Medium Data Bit Synchronizer (LMDR) 

R^ = 1 ksps to 10 Msps 

fn = 0.1% Rg, 0.3% R^, 1% R^ and 2% R^ 

Manufacturer: Aydin Monitor 
Jitter problem: No. 

First we consider the sinuscidc. jitter case. We consider two 
subcases; the first is for Af = 0.01% and f^^^ = 0.01% we ^all this 
case (a). The second is for Af = 0.1% R^ and f^^^ = 0.1% R^, which we call 
case (b). The BER degradations are shown in Table 6. 

For the random case using the same two subcases* we obtain the 
results for BER degradation which are shown in Table 7. 

In conclusion, we see that, when the frequency deviation and 
the deviation rate are specified at 0.01% R^, the degradation does not 
exceed 0.2 dB but, with the specifications at 0.1% R^, the degradation 
can be over 10 dB. Also from Tables 6 and 7, it can be seen that Man- 
chester data causes more degradation than NRZ data, as Is well known. 
Fina.ly, sine wave jitter and flat random noise jitter cause roughly 
the same level of degradation based on the Ku-band Shutter jitter 
specification. 

4.3.2 Effect of Data Asymmetry on Unbalanced QPSK Signals 

with Noisy Phase Reference 

In Appendix II, an analysis of the impact of data asymmetry on 
the bit error rate performance of a QPSK signaling scheme used to trans- 
mit two data streams with different rates and different powers is 
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Table 6. Degradation for the Three Bit Synchronizer Classes and 
the Two Jitter Specifications [(a) and (b)] for Sine 
Wave Jitter (Uncoded Data) 


Bit 

Synchronizer 

Class 

Degradation 


NRZ Data 

Bi-^-L Data 

UHDR (a) 

0.1 dB 

0.2 dB 

0.25 

UHDR (b) 

2.5 dB 

= 10.0 dB 

2.5 

HDR (a) 

<0.1 dB 

< 0.1 dB 

0.125 

HDR (b) 

2.0 dB 

“ 8.0 dB 

1.25 

LMDR (a) 

<0.1 dB 

< 0.1 dB 

0.1 

LMDR (b) 

1.7 dB 

4.5 dB 

1.0 


Table 7. Degradation for the Three Bit Synchronizer Classes and 
the Two Jitter Specifications [(a) and (b)] 

Random Noise Jitter (Uncoded Data) 


Bit 

Synchronizer 

Class 

Degradation 


NRZ Data 

B1-^-L Data 

UHDR (a) 

0.03 dB 

0.08 dB 

0.25 

UHDR (b) 

2.6 

dB 

> 10.0 dB 

2.5 

HDR (a) 

< 0.1 

dB 

<0.1 dB 

0.125 

HDR (b) 

1.6 

dB 

9.0 dB 

1.25 

LMDR (a) 

< 0.1 

dB 

<0.1 dB 

0.1 

LMDR (b) 

0.,’ 

dB 

2.3 dB 

1.0 


r 
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presented. Such a model arises* for example* In the QDSB system employed 
in the Shuttle Ku-band Mode 1 return link [11]. which Is essentially at 
cascade of tmo such unbalanced QPSK systems. Of primary Interest here 
will be the subcarrier QPSK system and* In particular* the version Imple- 
mented via a digital phase-shift-modulated square-wave subcarrier [12]. 

The performance of the recently popular unbalanced QPSK (UQPSK) 
systems In the presence of reference phase noise has been examined by a 
number of authors (see [13] to [17] and references therein). The general 
approach Is to evaluate the bit error probability conditioned on the data 
and the value ^ of the phase error* then average over these variables. 

The marginal probability density function (pdf) of ^ strongly depends 
on the employed tracking loop and is. In general* difficult to accurately 
evaluate for sophisticated receivers. Here we follow a truncated Taylor 
series expansion approach* Indicated in [15], which provides credible 
results* assuming high loop signal -to-noise ratio (SNR). Regardless of 
the particular method used, the bit error probability (BEP) depends on 
the powers and data rates of the individual channels. 

The i>-conditioned BEP (which is a reasonable estimate of BEP 
for very high SNR) is an increasing function of (j», the reason being that 
an imper'^ect phase lock loop attenuates the demodulated signal power and 
at the same time increases interchannel cross talk. 

Data asymmetry is a potential source of performance degradation 
in any digital transmission system. It arises whenever the modulator 
spends more time at one amplitude state than the other as a result of 
the misalignment of a threshold device. It is usually defined as the 
difference between the elongated pulse length and the shortened pulse 
length normalized by the nominal length. References [18] and [19] per- 
tain to previous work on the issue of determining the effect of data 
asymmetry on the BEP of BPSK signaling ([18] for a variety of coded and 
uncoded data) and UQPSK with equal data rates [19]. Since no cross talk 
was considered in [19], the two QPSK channels are essentially two inde- 
pendent identical BPSK channels with different power allocations. 

Appendix II addresses the problem of evaluating the perfor- 
mance degradation of a UQPSK system due to both data asymmetry and noisy 
reference. As we shall see, the degradation of each channel can quali- 
tatively be perceived as consisting of a "self-degradation" term plus a 


A/t 
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"cross-degradation" term. The former Is the superposition of a power 
attenuation (due to noisy reference) along with the impairment caused by 
asymmetry. The latter (which is actually the interchannel interference) 
exists because of the presence of noisy reference and also depends on the 
asymmetry of the companion channel data stream. Hence, the results to be 
derived here should coincide with those in [18] and [19], if the phase 
error is assumed to be zero. 

The parameter of interest here is the additional signal -to-noise 

ratio aSNR (d6) required to compensate for the losses due to asymmetry and 

phase noise for each channel. It is defined as the difference between the 

-5 

SNR required to achieve a BEP of 10 under the aforementioned conditions 
and the nominal value of 9.6 dB for an ideal BPSK system without impair- 
ments other than Gaussian noise. 

We have assumed a power ratio of 4 and a data rate ratio of 10 
for the two channels. These are typical values for the I and Q components 
of the subcarrier QPSK system of the Orbiter's Ku-band, which operate at 
2 Mbps and 0.2 Mbps, respectively. 

According to the approach taken in Appendix II, the phase error 

2 

impact is sufficiently characterized by the mean m. and the variance a. 

<p <p 

of the tracking jitter (j>. When a biphase Costas tracking loop is employed 
at the receiver, references [20] and [21] provide some (although compli- 
cated) analytic results for j^. We felt than an extension of these 

•p 

results in the presence of asymmetry is a prohibitively complicated task; 
therefore, we can only postulate some values for m. and a.. However, it 
is intuitively appealing to assume that, in the presence of small asym- 
metry, typical values for will be slightly higher than what appears in 
the figures of [20] and [21], As such, we have selected to examine the 
cases of m^ = 0°, 3°, 5“ and a.^ = 0,4,9 (degrees)^, 

9 9 

Before commenting on the computer results, let us focus on the 
assumptions made in the course of the analysis. First, we consider NRZ-L 
data for both channels. An analysis for biphase-L data would follow on 
the same lines. Second, we have assumed perfect symbol synchronization 
for each individual channel. Third, we have assumed that the data rate 
ratio N is much greater than one, in general (for the specific applica- 
tion where N= 10, such an assumption is well justified). The above 
assumptions suffice in order to derive a good estimate of the bit error 
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probability of the high data rate channel. For the low data rate channel, 
two more assumptions are involved: first, that N is an integer and, sec- 
ond, that the two channels are aligned, i.e., no epoch difference exists. 
Although the second assionption is perhaps rarely met in practice (since 
there are two different data clocks for the two channels, they most likely 
are independent, therefore not aligned), one can argue that the signifi- 
cance of the epoch difference diminishes as N increases. Similarly, if N 
is not an integer, one can consider the integer part of it and still get 
credible results. In summary, the cardinal assumptions for the analysis 
to be reasonably precise are that (1) each channel is perfectly synchron- 
ized, and (2) N is sufficiently greater than one. The low data rate 

The low data rate channel (Q) is found to be more sensitive to 
parameter variation than the high data rate channel (I). In general, the 
SNR losses increase exponentially with the amount of asymmetry present 
for both channels. For the values considered. Figure 7 in Appendix II 
shows the range from 0 dB to 1.4 dB for the I channel and from 0 dB to 
1.9 dB for the Q channel. As shown in Figure 8 of Appendix II, the 
effects of phase mean and phase variance are effectively decoupled. How- 
ever, as expected, as increases, the effect of phase offset becomes 
more and more drastic. 
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5.0 S-BAND PAYLOAD COMMUNICATION SYSTEM STUDY 

Three major system interface Issues in the payload cofnmuni cation 
system required a detailed investigation. These interfaces are the pay- 
load interrogator (PI) and Ku-band signal processor (KuSP)* the sweep 
acquisition of the payload by the PI, and two-way phase noise generated 
in the PI and turned around by a coherent transpon ter. 

5.1 Bent-Pipe Performance with the Current PI/KuSP 

RMS Regulator Capability 

The maximum allowable deviation of the Ku-band FM transmitter 
is ±11 MHz. Optimum bent-pipe performance is obtained whenever the modu- 
lating waveform is allowed to produce this maximum deviation in accord 
with its peak values. For this reason, Axiomatix has proposed that a 
peak type rather than an RMS type of AGC regulating loop be employed. A 
comparison between the RMS and peak regulator performance, as measured by 
FM transmitter mean deviation, is sumnarized in Table 8. 


Table 8. Bent-Pipe FM Transmitter Mean Deviations (MHz) 


Waveform 

RMS Regulator 

Peak Regulator 

One Sinusoid 

3.7 

7.8 

Two Sinusoids 

3.7 

5.5 

Three Sinusoids 

3.7 

4.4 

Four Sinusoids 

3.7 

3.9 

Square 

3.7 

11.0 

Gaussian 

3.7 

3.7 


The peak regulator optimizes the FM deviation for all waveforms, 
providing maximum bent-pipe SNR performance for all waveform conditions. 
The RMS regulator, on the other hand, may have its output scaled to pro- 
vide optimum deviation for any chosen waveform, but the deviation for 
all other waveforms will be suboptimum. If, for example, the RMS regu- 
lator output is optimized for the Gaussian waveform (characteristic of 
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PI low SNR conditions), the transmitter will be underdeviated for all 
other waveforms. As part of the development and specification of the 
PI/KuSP bent-pipe Interface, Axiomatlx undertook the design and evalua- 
tion of a signal waveform peak regulator. A final report on the entire 
activity appears In Aopendlx III. 

Now suppose that the RMS regulator Is used and Its output Is 
scaled to provide optimum deviation for a single-sinusoidal waveform. 

From Table 8, the mean deviation will be 7.8 MHz. If a Gaussian waveform 
rather than the single sinusoid then appears at the Input to the RMS reg- 
ulator, the RMS regulator will automatically scale the Gaussian waveform 
to cause a mean deviation of 7.8 MHz. But the peak-to-peak RMS ratio of 
a Gaussian waveform Is on the order of 3:1; therefore, the peak deviation 
will be 23.4 MHz, or more than two times larger than the maximum deviation 
limit. 

To prevent overdeviation, an amplitude clipper Is usually 
employed at the Input to the FM transmitter. Thus, for the Gaussian wave- 
form example just cited, it will be clipped at its 1 .4o level, causing 
extreme distortion and SNR loss. 

For the bent-pipe, the waveform condition for which the RMS 
regulator output should be scaled to provide optimum deviation is an open 
issue. The final choice will have to be a value judgment. An expedient 
solution may be to optimize the bent-pipe performance for the two-sinusoid 
case (as this may be the most likely bent-pipe signal). Thus, the mean 
deviation would be set at 5.5 MHz and amplitude clipping would keep any 
other waveform, such as the Gaussian, from exce'^Jing the 11 MHz limit. 

The peak-to-RMS ratio would therefore become 2. Underdeviation would 
occur for single-sinusoid and square waveforms, while the Gaussian wave- 
form would be clipped at its 2a level. It is interesting that the cur- 
rent TRW RMS regulator and output circuit design provides for a 1.5:1 
peak clipping-to-RMS ratio. TRW, however, would like to change the max- 
imum to 8V peak-to-peak, achieving the 2:1 ratio. 

SNR performance estimates for the underdeviation cases may be 
easily calculated. The peak clipping cases, however, are not readily 
amenable to analysis (although performance bounds may be obtained by 
analysis). Simulation or actual hardware measurements must therefore 
be used to obtain quantitative figures for the clipping cases. 
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5.2 Sweep Acquisition of Deep-Space (PS) Payloads by the PI 

On June 22, Axiomatlx completed the effort to assess the proper 
sweep rate for the DS transponder. A formal viewgraph presentation was 
subsequently prepared and oresented at RI on June 26. This meeting was 
attended by RI, NASA and Axiomatlx personnel. 

The viewgraphs used at the June 26 meeting appear as Appendix IV 
to this report. It is noted, however, that an error appeared on the last 
page of the presentation, resulting in a too-fast sweep rate of 400 Hz/s 
being agreed upon at the June 26 meeting. As a part cf the resolution 
process. It was decided by all present that the DS transponder would never 
be operated In conjunction with the PI transmitter +4 dBm power level at 
a range of 10 nml. Thus, the sweep rate corresponding to the +27 dBm 
level was selected. 

At the TRW monthly program review on July 12, Axiomatlx corrected 
the aforementioned error. TRW also estimated that, within the neighbor- 
hood of the proper sweep rate, a tolerance as large as ±30% (rather than 
±20%, as assumed on June 26) of the nominal sweep rate value could be 
expected. The final report, therefore. Is that the nominal sweep rate 
must be set at 250 Hz/s. 

Appendix IV to this report Is a set of viewgraphs which define 
the problei'i. Indicate the nature of the problem solution, show the solu- 
tion as a function of transponder operating conditions and tolerances, 
and present recommended sweep rate limits. The approach employed to ob- 
tain the critical sweep rate values is a modified method of phase-plane 
analysis that apepears in Appendix V. Since the calculations are based 
on an ''Ideal' second-orde** PLL transfer function, a 20% backoff has been 
made to allow for the mechanization imperfecuions of the DS transponder. 

5.3 Two-Way Phase Noise Performance for PI 

and Conerent Payload Transponders 

In order to obtain a realistic characterization of system per- 
formance, it is essential to analyze all of the many error sources, in- 
cluding phase noise. Phase noise is genera^^ed in the transmitter refer- 
ence oscillators and the receiver local oscillators. Other sources which 
are not as obvious are in the mixer (multiplier) logic noise, AM-to-PM 
conversion due to link nonlinearities, and possible vibration effects. 
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The result of this phase noise Is degradation of the carrier tracking 
loop performance, and this contributes to the overall system error prob- 
ability performance. 

Appendix VI contains a detailed discussion of the effects of 
phase noise using mathamatical models of the carrier tracking lOops and 
Introducing the appropriate power spectral densities of the various phase 
noise Instabilities. Section 2.0 of this appendix outlines the system 
model of a two-way coherent communication system and describes the basic 
system configurations. Section 3.0 calculates the mean-square phase 
noise for a one-way coherent link, while Section 4.0 extends this elabo- 
rate analysis to the two-way link. Incorporating both In-band and out-of- 
band phase noise sources. For the special situation of a close range 
communication link. Section 5.0 deals with the essentially coherent case 
of a common transmitter and receiver reference oscillator, where the cor- 
relation time of the phase noise Is short compared t'j the round-trip delay 
times. The conclusion, as might be expected. Is a slight modification of 
the noncoherent two-way link. Section 6.0 expands the analysis to Include 
other Internal frequency synthesizer oscillator sources by modeling and 
Incorporacing these synthesizer phase noises. Section 7.0 extend? the 
?tudy to the three- and four-way mean-square phase noise situations. 

next appendix. Appendix VII, describes a method of charac- 
terizing the phase noise spectra by a mean-square error approximation 
using measured laboratory results. Two examples of phase noise modeling 
for specific actual phase noise measurements are outlined, and compari- 
sons are Included in both cases which show; the derived and expected 
results. 

Finally, Appendix VIII calculates the expected phase noise 
RMS deviation upper bound of 9® for both the SGLS and STDN modes using 
the mean-square phase noise relationships of Appendix Vi and the mean- 
square error approximation of Appendix VII. 
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APPENDIX I 

KU-BAND BIT SYNCHRONIZER BER DEGRADATION 
DUE TO DATA BIT CLOCK JITTER 


By 

Jack K. Holmes 



1.0 SUMMARY 

This report determines the timing error variance and resulting 
bit error rate degradation in a Ku-band Shutt1e*compatib1e second-order 
bit synchronizer due to input data clock jitter modeled as either a sinu- 
soidal signal or a spectrally flat random process. Both NRZ and Manchester 
data are considered. 

Results are obtained for the bit synchronizer tracking error 
for both the sinusoidal and random jitter cases. Next, the bit error 
probability for both sinusoidal and Gaussian data clock jitter is derived 
and plotted. Using the timing variance, it is then possible to compute 
the bit error probability and the associated degradation from the jitter- 
free case. 


Finally, the three bit synchronizer types used with Ku-band data 
are evaluated to determine their respective bit error rate degradation. 
Depending on which specification, modulation format and synchronizer are 
considered, bit error rate degradations vary from less than 0.1 dB to 
over 10 dB. The degradations tend to be small (0.1 dB) when the frequency 
deviation and deviation rate are specified at 0.012 but quite large 
when specified at 0.12 where R^ is the symbol rate. 

2.0 PROBLEM STATEMENT 

The problon this report addresses is the determination of the 
bit error rate degradation caused by data clock timing errors that are 
either sinusoidal or a random process. This specification appears In basi- 
cally similar forms in numerous documents. From the Ku-band specifica- 
tion [1] (Appendix VII, p. 456), we have: 

"70.3.3.2.2 Data and Clock Jitter . Data and clock jitter 
shall be taken as the time-varying component of data and/or 
clock transitions compared to a jitter- free data stream or 
clock. The definition contains both a jitter magnitude and 
a rate in the following form: 

(a) Jitter magnitude shall be taken as a percentage of the 
data symbol clock rate. The resulting magnitude shall be 
the peak magnitude for sinusoidal jitter or the 3-sigma 
value for random jitter. 

(b) Jitter rate shall be taken as a percentage of the data 
symbol clock rate. The resulting jitter rate shall be the 
peak rate for sinusoidal jitter or the 3-sigma value when 
the jitter is random. 
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When the data are coded, bit Jitter shall refer to symlK>1 
Jitter. Also when a clock accompanies the data, the Jittar 
shall refer to clock Jitter as well as data Jitter." 

Also from [1] (p. 147), the high data rate (¥DR) Input to the 
SPA has a Jitter specification: 

“HDR Payload Digital Data Input to SPA 

h. Frequenc;y Jitter: plus or minus a frequency deviation 
of 0.01X rms of Uie data rate at a rate of 0.01X rms 
of the data rate." 

Other specifications on tills topic Indicate a 0.1X of the data 
rate deviation and a 0.1X of the data rate deviation rate. 

We will address the sinusoidal jitter case first, then consi >r 
the “random" case. 

3.0 JITTER MOOa DEDUCED FROM SINE WAVE SPECIFICATION REQUIREMENTS 
Consider a phase^modulated sine wave clock modeled by 

y(t) * i/?A sin (u^t e(t)) (1) 

where A is the rms signal amplitude, <dq Is the angular clock frequency 
(rad/sec) and e(t) is the phase modulation process which causes the jitter. 
In Figure 1, a typical NRZ data stream and the associated (unmodulated) 
sinusoidal clock are illustrated. Since the phase modulation process, 
e(t), is assumed to be sinusoidal, we obtain 

e(t) = 0Q sin (u^t) (2) 

where is the peak phase modulation of the clock and is the frequency 
of the phase deviation. Now, since the jitter magnitude is specified as 
a percentage of the data symbol clock rate (frequency an^litude), we have 

o(t) = 9 q 0 !^ cos (oi^t) (3) 

and the peak is some percentage of the data symbol clock rate, i.e.. 


Bit Period = Clock Period 



Figure 1. NRZ Bit Sequence and the Associated Clock 
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0^ . 2. if - MR, (4) 

where D Is the relative frequency deviation relative to the data rate. 

From Uie HDR specification [1] (Section II), we deduce that 0 ° 

0.0001 (0.01%). The jitter rate specification implies that 

“ta • (5> 

where is such that = F^R^ and is the relative modulating fre- 

quency (relative to R^). The jitter rate specification specifies the rate 
of the sinusoidal jitter. Now, using (4) and (5) in (2) yields 

0(t) = sin {bi^t) = sin (2wF|^R^t) (6) 

in m 

When the relative jitter magnitude (D) and relative jitter rate (F^jj) 
are equal, we obtain 

0(t) = sin (2ir y) (7) 


which inJicates a peak one radian phase deviation or 1/2ir = 16% of a bit 
period timing jitter (see Figure 2). Other values of D and F^^^ yield dif- 
ferent phase deviations. In the derivation of (6) and (7), the peak jit- 
ter magnitude (frequency deviation) was used but, for the jitter rate, 
only one frequency exists so that the meaning of peak rate is not clear. 
It is felt that (2) is the only meaningful model of sinusoidal jitter. 

The clock is therefore modeled as 


y(t) 


/Z A sin 


(“o‘ ♦ ^ (2” v>)- 


( 8 ) 


If the jitter frequency, F^j^R^, is very low compared to the bit 
synchronizer natural frequency, f^j, all the "jitter" will be tracked out. 
On the other hand, if the jitter frequency Is high compared to the natural 
frequency of the bit synchronizer, then all the "jitter" will be passed 
through, with a resulting high BER degradation. 
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H K 


Clock Deviation 


r ■ " - " - - T “i 

I I 

I i 



Figure 2 . Illustration of One-Radian Clock Deviation 
and the 16^ Bit Time Deviation 
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We now determine the response of a second>order bit synchronizer 
with a one radian peak phase modulation on the Input. In general, H(S), 
the closed-loop transfer function. Is given by 


H(S) 




S +25 


(9) 


where Is the loop natural frequency (rad/sec) and 5 Is the damping 
factor. We wish to evaluate the peak squared tracking error expressed 
by 



|H(f)|^ 


f = f. 


(10) 


where S * 2n1f, (|>^ is the input peak phase error and ♦p is the residual 
peak phase out of the bit synchronization loop. If we select 5 * 
a commonly used value, then It Is easy to show that 




f^ 

-3 T 

^ * V 


(11) 


where f^^ = natural loop frequency expressed In Hz. From 

(10), using the Input phase modulation process of (6), we have 


or 


f' 

7J fUf* 


m 


f = F R 
m s 


2 ^ D' 


m 


m s 

F ^ R ^ + f ^ 
ms n 


( 12 ) 


(13) 


or 
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p 2 2 

. 2 _ « 2.x2 '^in '^s 

♦p 'p'’ rrrr- 


F ’ R ’ + f 
ms n 


Now we fix r , R- and D, and let vary* so that using F^ R^ = f 
ns m'' ^rasm 


-X . — 

R.V + f.^ 


m n 


If we differentiate the right-hand side of (14) and set It equal to zero, 

we obtain the solution, f^ * f . It therefore follows that the maximum 

m n 

error occurs at f„ « f when the deviation &f = R.D is fixed and the 
m n s 

error at that frequency is given by [using (2), (4) and (15)] 


where 


R V 


P f =f 
m n 


It therefore follows that 


* ~ 0n f„ = f„; &f fixed 

Tn — n m n 


Whenever f„^f, the peak value is less for fixed Af. Figure 3 illus- 

m n 5 2 

trates the response <l'p^/Af as a function of f^^^ for fixed af. This curve 
is related to BER degradation since a large value of peak phase error, ^p. 
causes a large BER degradation. Blyth and Carpenter [2] have summarized 


We are fixing the ^^requency deviation and clock nominal 
frequency and varying the phase-modulating ^^requency. 
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measured BER degradation performance of a bit synchronizer* due to jitter 

and pointed out under a constant value of Af the degradation peaks very 

near f • f ! 
m n 

It should be further pointed out that the peak phase error. 
for a fixed input phase modulation deviation, 6q, Is maximum for frequen- 
cies much greater than f^^^ • f^. This fact is easily demonstrated by 
starting with (14) and letting f^^^ = Kf^^, so that we have 



(19) 


so that we obtain 


♦ 


2 

P 




+ 1 


, 2 _JL 
0 


( 20 ) 


2 2 

Hence, as K = ■*■ “» i'p " ®o ’ 

In Figure 4, the peak output phase error is plotted versus 
W^n ^ second-order (high gain) bit synchronizer loop with 5 = 0.7G7 
which has a unit input peak phase deviation. 

In conclusion, the phase modulation of the tone having a period 
equal to the duration of the bit time that satisfies the jitter magnitude 
and jitter rate specifications is given by 


e(t) = sin (u^t) 
m 


( 21 ) 


where Af = DR is the peak jitter magnitude and f = o)„/2ti = F R is the 
s rw j mm ms 

jitter rate, with R^ the symbol rate, 0 the '•atio of frequency deviation 


tc data rate and F^ the ratio of modulation frequency to data rate. 


Hence, 
is the peak 


when 0 and F„ are obtained from the specification's = D/F„ 
m 0 2 2 *” 

phase iTiodulation, and Figure 4 allows one to compute <I'p^/ 6 q^ (f ^>6 ratio 


Based on the available literature, a common test procedure is 
to vary Af for fixed f^ and measure the BER degradation. 
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of bit synchronizer loop error squared to Input phase squared) when 
^n/^n specified. Recall that f^ Is the bit synchronizer loop natural 
frequency (rad/sec). 

4.0 PROBABILITY OF ERROR FOR NRZ AND MANCHESTER DATA VERSUS 

SINUSOIDAL PEAK TIMING ERROR 

We now determine the probability of error versus Ej,/Nq for sev- 
eral values of peak timing error based on a sinusoidal timing error vari- 
ation with time. The motivation for developing the following results Is 
based on appraising the effects of the jitter magnitude and jitter rate 
when the jitter Is sinusoidal, as specified In the Ku-band specifications. 

It will be assumed that the bit clock jitter rate Is slow com- 
pared to the data rate so that the average bit error rate can be described 
as 


PE = f P(E/t) p(t) dt 
•00 


( 22 ) 


where P(E/t) is the bit error rate conditioned on the value of the timing 
error, t, a.id p{x) is the probability density function of the error timing 
error r. First consider P(E/x) the conditional probability of bit error 
given .. parcicular value of timing error x for NRZ data. For uncoded NRZ 
signals, there are two cases for the bit error probability. In the first 
case, when no transition in the bit stream occurs, the probability of 
error is simply 



where 

7 , - t ^/2 

Qv“) = / — ® 


( 24 ) 
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When a transition does occur, the bit error probability is no longer 
expressed by (23) but, rather, is modified to 



The reason for the factor (1 - 21t|) is due to the previous bit (or fol- 
lowing bit) being of the opposite sign and being included in the integra- 
tion of the current bit duration of T seconds (see Figure 5). The matched 
filter integrates over T seconds. As the bit timing error increases, the 
integration includes the present bit plus the adjacent bit. When thore 
is no transition, this error produces no decrease in the integrated sig- 
nal voltage. However, when there is a transition, there is a decrease in 
the integrated voltage during one bit time. In fact, with a liitle reflec- 
tion, it Is clear that the decrease in integrated voltage is linear and 
goes to zero when x = ±T/2. It follows that the integrated voltage con- 
tains the factor (1 - 21x1) for |xl <1. 

In conclusion, if we denote by Pj the transition probability, 
the conditional probability of error conditioned on x is given by 



1x1 <0.5 (26b) 


where the Manchester case can be derived in a manner similar to the NRZ 
case. 

Next, we obtain the probability density function of the timing 
error, p(t). We descriue the equivalent phase modulation by 


"{•(t) = <|>p sin (oj^jjt + 9 q) 


(27) 



13 


(a) No Transition 




Figure 5. The Bit Integration Region in the Case of 
No Transition and a Transition with Timing Error t 
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where (denoting the data rate by R), is the peak phase error (rau), 
is the modulation frequency « 2 wR^) and 0 q is a random variable uni- 
form over (-ir,v). 

Since, for any t, + 3 q is uniformly distributed over (-w,*), 
we model 4(t) as far as the distribution of ♦(t) goes by 

♦(t) » *p sin 0Q . (28) 

where Sq is a uniform random variable, uniform on (-n,ii). To obtain the 
density function of we have [3], 

P(^) = ^29) 

.1 2 J y /. 2 J 

2tt Jp - 2ir/»p - ^ 

SO that 

p(^) = J _ (30) 

Now since 4>p/2ii = Xp, the peak timing error, we have, with another change 
of variablt?s. 


p(t) 



-■fp 1 ^ 1 ^ 


P 


(31) 


Hence, the average probability of error is given by 



and 
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MANCHESTER: PE 



Equations (32a) and (32b) were programned on a digital computer, with the 
results Illustrated In Figures 6 and 7. Notice that a 52 (0.05) peak sine 
wave timing error causes about 0.5 dB degradation at PE - 10*^ for NRZ data. 
For Manchester data, a 5% peak causes about 1 dB degradation. 

Although It Is not obvious from the curves, there Is no 1rreduc> 
Ible error In this sinusoidal timing error case as there exists In the 
Gaussian error case. 

A comparison of the Mancl^ster case can be made with sane experi- 
mental measuren^nts of J. Rivers [5] of the NASA Johnson Space Center. 

In Figure 8, Rivers-measured BER degradation and the theoretical values 
obtained from Figure 7 are compared for Manchester data at PE = 10~^. As 
can be seen from Figure 8, the theory Is quite accurate, having at most a 
0-3 dB error. 


5.0 JIHER MODEL DEDUCED FROM THE 

GAUSSIAN SPECIFICATION REQUIREMENTS 

Now we consider the Gaussian modulation case and attempt to 
determine a model of the jitter process. Again consider the sinusoidal 
clock having a period equal to the bit duration. We have for the clock 
signal 


y(t) = A sin (oigt + e(t)j (33) 

where e(t) is assumed to be a stationary Gaussian random process. Since 
the 3-sigraa jitter is specified to be a percentage of the data rate, we 
have that 


3^ f ^e{f) df = (DR^)^ 


(34) 


Unless the phase process has a specified spectral density, it is impossible 
to find the 3-sigma modulation frequency. In order to proceed, we will 






2 3 4 5 6 7 8 9 10 11 12 13 14 

R(dB) 

Figure 7. PE versus Peak Sine Wave Timing Error--Sinusoidal 
Bi-i|i-L Data, Bit Synchronizer Tracking 





Tp (fraction of a bit) 

Figure 8. BER Degradation versus tp 
Manchester Data, Sinusoidal Jitter 
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assume a phase process spectral shape that Is flat In frequency. It 
should be pointed out that this Is a “reasonable" assumption, but Is not 
based on measurements. Figure 9 Illustrates the phase process spectral 
density. We have assumed that all the phase spectra Is contained In 
Hz. Since o>^(f) Is the spectra of the phase rate or frequency 
(I.e., Jitter), we have 


(2ir)2 0 ^/ * f e/*e(f) df - f i/jjeif) df 


(35) 


where Is the rms frequency deviation (Hz). Evaluating (34) and (35) 
with the phase process model of Figure 9 yields 


ill s 


9 . (2ir)^ = 9 . / (2ir)2 f^ df = (2ii)^ (DR^)^ (36) 

m s 


or 


2 (2ir)2 9 Nq ^^'" 3 - -^ - = (2ir)2 (DRj.)^ 


(37) 


Solving for Ng yields 


6(F R ) F 
' m s^ ro 


0 ^ 

* 7T * 


-DR^ < f < DR^ (38) 


Now that N- the spectra! density of the. phase proces<5. is determined, we 
“ * 

can compute nrs phase deviation: 


DR. 


1 




( 39 ) 


s ' T s' \m 



Figure 9. Model of the Phase Process 
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or 


a 


2 

8 



r«ids^ 


(40) 


Notice again* when 2 > F , o » AJi radians. However, the Input phase 

IH d 

deviation Is not nearly as Important as the phase error variance of the 
bit synchronizer given by 

/ |1 - H(f)|2 i2{(f) df « a 2 (41) 


where H(f) Is the closed-loop transfer function of the bit synchronizer. 
Evaluating (39) for a second-order bit synchronizer with 5 * 0.707, we 
have 



F R 



^ m s 

f^ 

1 

-F R 
m s 





df 


(42) 


Note that 





df 


2f 


m 




f /f 
m' n 


dx 


1 + x" 


(43) 


So from (42), using (43), we arrive at 

(44) 


The integral of (44) was evaluated by computer since the tables 
appeared to give an incorrect result for values of f^/fp larger than one. 
The curve is shown in Figure 10. 



0 


m n 


dx 

1 + x^ 



1.1 
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For values of f^f,, < 0.5, (44) can be well approximated by 


2 




(45) 


2 

Hence, when and D are known, can be evaluated from (45) 
or (44) along with Figure 10. Note, to convert to o^ (timing error), 
we have the simple relationship 


® ^ (fractions of a bit). (46) 

T fclT 

6.0 PROBABILITY OF ERROR FOR NRZ AND MANCHESTER DATA VERSUS 
THE RMS TIMING ERROR FOR A GAUSSIAN PHASE PROCESS 

Since the rms timing error has been computed in (44) and (45) 
using (46), we can now determine the probability of error assuming thdt 
the phase process maximum modulating frequency is much less than the 
data rate, i.e.. 


m 


m s 


(47) 


Just as in Section 4.0, it can be shown that the conditional 
error probability P(E/x) is given by 


NRZ: P(E/t) 

MANCHESl ER: 





I'll’ 



P(E/t) • (1-Pt) Ql (i-4|i|)) + p^q( ^ (1-2|t| 


T 1/2 (48) 


where x is the timing error, Pj is the transition probability and Q(x) is 
defined by 
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Q(x) • f 


■t^/2 


e 


dt 


(49) 


The probability density function of the timing error Is given 
by 

P(*f) • exp 

/2ito_ 

T 

It follows from (48) and (50) that» when Py » 1/2; 



(52) 

These functions are plotted In Figures 11 and 12 for the NRZ 


and Manchester (Biphase-L) symbol format case. Note that a 2.5% (0.025) 
Gaussian timing error produces 0.4 dB degradation for NRZ and about 1 dB 
degradation for Manchester data. It can be shown that irreduciD- error 
probabilities occur in the Gaussian timing error case, causing very high 
degradations at larger timing errors. 

7.0 EVALUATION OF THE BER DEGRADATION FOR THE SHUTTLE KU-BAND 

RETURN LINK BIT SYNCHRONIZER 

In this section, we apply the results of previous sections in 
order to estimate the bit error rate (BER) degradation due to either sinu- 
soidal or a flat Gaussian phase modulation process. To evaluate the BER 
degradation, it is necessary to k'^ow the bit synchronizer loop bandwidths. 







Figure 12. PE versus RMS Timing Error — Gaussian Error 
Bi-(j>-L Data, Bit Synchronizer Tracking 
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Based on a telephone conversation with John Roach of Harris 
Corporation (Helbourne* Florida), tiie following bit synchlK>n1zer specifi- 
cations were obtained: 

1. Ultra-High Data Rate Bit Synchronizer (UWWO 

- 75-150 Msps 
f„ 5 60 kHz (B|^ = 200 Hz) 

Manufacturer: Motorola 

Jitter problems comnents: Yes, since they could not 

make the loop bandwidth narrow 
enough due to Internal delays 

2. High Rate Bit Synchronizer (HDR) 

R * 10-75 Msps (will operate down to 4 Msps) 
f^ = 0.08* Rj (B^ * 0.267* R^) 

Manufacturer: Harris Corporation 
Jitter Problem: No 

3. Low-Medium Data Bit Synchronizer (LMDR) 

Rj = 1 ksps to 10 Msps 

f = 0.1* R., 0.3* R_. 1* R., 2* R 

II 9 9 w w 

Manufacturer: Aydin Monitor 
Jitter Problem: No. 

First we consider the sinusoidal jitter case. We consider two 

subcases; the first Is for af = 0.01* R^ and f^ = 0.01* R^, we call this 

case (a). The second Is for Af = 0.1* R^ and f^^^ = 0.1* R^, which we call 

case (b). The BER degradations are shown In Table 1. 

For the random case using the same two subcases, we obtain the 
results for BER degradation which are shown In Table 2. 

In conclusion, we see that, when the frequency deviation and the 
deviation rate are specified at 0.01* R^, the degradation does not exceed 
0.2 dB but, with the specifications at 0.1% R^, the degradation can be 
over 10 dB. Also from Tables 1 and 2, It can be seen that Manchester 
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Table 1. Degradation for the Three Bit Siynchronlzer Classes and 
the Two Jitter Specifications tU) and (b)3 
for Sine Wave Jitter (Uncoded Data) 


m 

Synchronizer 

Class 

degradation 


NRZ Data 

B1-4*L Data 

UHDR (a) 

0.1 dB 

0.2 dB 

0.25 

UHDR (b) 

2.5 dB 

II 

O 

2.5 

HDR (a) 

<0.1 dB 

< 0.1 dB 

0.125 

HDR (b) 

2.0 dB 

= 8.0 dB 

1.25 

LMDR (a) 

<0.1 dB 

< 0.1 dB 

0.1 

LMDR (b) 

1.7 dB 

4.5 dB 

1.0 


Table 2. Degradation for the Three Bit Synchronizer Classes and 
the Two Jitter Specifications t(a) and (b)] 

Random Noise Jitter (Uncoded Data) 


Bft 

Synchronizer 

Class 

Degradation 

y^n 

NRZ Data 

Bi-4-L Data 

UHDR (a) 

0.03 dB 

0.08 dB 

0.25 

UHDR (b) 

2.6 dB 

> 10 dB 

2.5 

HDR (a) 

< 0.1 dB 

< 0.1 dB 

0.125 

HDR (b) 

1.6 dB 

9.0 dB 

1.25 

LMDR (a) 

<0.1 dB 

< 0.1 dB 

0.1 

LMDR (b) 

0.7 dB 

2.3 dB 

1.0 
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data causes more degradation than NRZ data, as Is Mell-known. Finally, 
sine wave jitter and flat random noise jitter cause roughly the same 
level of degradation based on the Ku-band Shuttle jitter specification. 
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1.0 INTRODUCTION 

In this report, an analysis is attempted of the impact of data 
asymmetry on the bit error rate performance of a QPSK signaling scheme 
used to transmit two data streams with different rates and different 
powers. Such a model arises, for example, in the QDSB system employed 
in the Shuttle Ku-band Mode 1 return link [1], which is essentially a 
cascade of two such unbalanced QPSK systems. Of primary interest here 
will be the subcarrier QPSK system and, in particular, the version impl 
mented via a digital phase-shift-modulated square-wave subcarrier [2]. 

The performance of the recently popular unbalanced QPSK (UQPSK) 
systems in the presence of reference phase noise has been examined by a 
number of authors (see [3]- [7] and references therein). The general 
approach is to evaluate the bit error probability conditioned on the 
data and the value (}> of the phase error, and then average over ^hese 
variables. The marginal probability density function (p.d.f.) of (J) 
strongly depends on the employed tracking loop and is, in general, 
difficult to accurately evaluate for sophisticated receivers. Here 
we follow a truncated Taylor series expansion approach, indicated in 
[5], which provides credible results, assuning high loop signal-to- 
noise ratio (SNR). Regardless of the particular method used, the bit 
error probability (BEP) depends on the powers and data rates of the 
individual channels. 

The ({i-conditioned BEP (which is a reasonable estimate of BEP 
for very high SNR) is an increasing function of 4>, the reason being 
that an imperfect phase lock loop attenuates the demodulated signal 
power and at the same time increases interchannel cross talk. 

Data asymmetry is a potential source of performance degradation 
in any digital transmission system. It arises whenever the modulator 
spends more time at one amplitude state than the other, as a result of 
the misalignment of a threshold device. It is usually defined as the 
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difference between the elongated pulse length and the shortened pulse 
length, normalized by the nominal length. References [8] and [9] per> 
tain to previous work on the issue of determining the effect of data 
asymmetry on the BEP of BPSK signaling ([8], for a variety of coded and 
uncoded data) and UQPSK with equal data rates [9]. Since no cross talk 
was considered in [9], the two QPSK channels are essentially two inde- 
pendent identical BPSK channels with different power allocation. 

The present work addresses the problem of evaluating the per- 
formance degradation of a UQPSK system due to both data asymmetry and 
noisy reference. As we shall see, the degradation of each channel can 
qualitatively be perceived as consisting of a "self-degradation" term 
plus a “cross-degradation" tenn. The former is the superposition of 
a power attenuation (due to noisy reference) along with the impairment 
caused by asyinnetry. The latter (which is actually the interchannel 
interference) exists because of the presence of noisy reference and 
also depends on the asymmetry of the companion channel data stream. 
Hence, the results to be derived here should coincide with those in 
[8] and [9], if the phase error is assumed zero. 

Section 2.0 sumnarizes the results of this current research 
and emphasizes the assumptions involved in deriving them, since these 
assumptions weigh heavily on the validity of the conclusions. In 
Section 3.0, we elaborate on the system model assuned by the analysis 
and, in Section 4.0, we proceed to evaluate the BEP for Channels 2 and 3. 


2.0 SUMMARY 

The parameter of interest here is the additional signal-to- 

noise ratio ASNR (dB) required to compensate for the losses due to 

asymmetry and phase noise for each channel. It is defined as the 

-5 

difference between the SNR re'^uired to achieve a BEP of 10 under the 
aforementioned conditions -<«d the nominal value of 9.6 dB for an ideal 
BPSK system without impairments other than Gaussian noise. 

We have assumed a power ratio of 4 and a data rate ratio of 10 
for the two channels. These are typical values for the I and Q com- 
ponents of the subcarrier QPSK system of the Orbiter's Ku-band, which 
operate at 2 Mbps and 0.2 Mbps, respectively. 
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According to the approach taken In this report, the phase 
error Impact Is sufficiently characterized by the mean m^ and the 
variance of the tracking Jitter (j). When a biphase Costas tracking 
loop Is employed at the receiver, references [10] and [11] provide 
some UHhough complicated) analytic results for o^. We feel that an 
extension of these results In the presence of asymmetry Is a prohibi- 
tively complicated task; therefore, we can only postulate some values 
for and a^. However, It is Intuitively appealing to assume that. 

In the presence of small asymmetry, typical values for will be 
slightly higher than what appears In the figures of [10] and [11]. 

As such, we have selected to examine the cases of m. = 0‘’,3°,5° and 

2 2 * 

= 0, 4, 9 [degrees) . The analytic results are summarized In (27a), 

(32), (39), and (48) for the high data rate channel 2 and In (27b), 

(32), [43), (58), (62), and [63) for the low data rate channel 3. 

Before commenting on the computer results, let us focus on the assimp- 

tlons made In the course of the analysis. 

First we considered NR2-L data for both channels. An analysis 
for biphase-L data would follow on the same lines. Second, we have 
assumed perfect symbol synchronization for each individual channel. 
Third, we have assumed that the data rate ratio N is much greater than 
one, in general (for the specific application where N = 10, such 
an assumption is well justified). The above assunptions suffice in 
order to derive a good estimate of the bit error probability of the 
high data rate channel. For the low data rate channel, two more 
assumptions are involved: first, that N is an integer, and second, 

that the two channels are aligned, i.e., no epoch difference exists. 
Although the second assunption is perhaps rarely met in practice 
[since there are two diffe’ jnt data clocks for the two channels, they 
most likely are independent, therefore not aligned), one can argue 
that the significance of the epoch difference diminishes as N increases. 
Similarly, if N is not an integer, one can consider the integer part 
of it and still get credible results. In summary, the cardinal 
assumptions for the analysis to be reasonably precise are that (1) each 
channel is perfectly synchronized, and (2) N is sufficiently greater 
than one. 
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Figures 7(a-c] and 8(a-c) show the plots of ASNR as a function 
of the percent asymmetry A. with m^ and var (|) as parameters. To gain 
some confidence about their correctness, let us notice that, for var (^-0 
and A»0 0"^ Is the single parameter then), the results coincide with 
the results of Figures 9 and 10 of [6] for both channels. Also, for 
the high data rate channel with m^»0, o^*0, and different degrees 
of asymnetry, the results coincide with Table 1 of [8]. 

A study of Figures 7 and 8 shows that the low data rate channel 
(Q) Is more sensitive to parameter variation than the high data rate 
channel (I). In general, the SNR losses Increase exponentially with 
the amount of asymmetry present for both channels. For the values 
considered, they range from 0 dd to 1.4 dB for the I channel and from 
0 dB to 1.9 dB for the Q channel. We also notice that the plots are 
virtually parallel (for a fixed o^), which means that the presence of 
a mean phase error shifts the curves upwards. In other words, the 
effects of phase mean and phase variance are effectively decoupled. 
However, as expected, as increases, the effect of phase offset 
becomes more and more drastic. 

3.0 ANALYTIC MODEL 

Vjo adept in this report the following modulation and demodu- 
lation schemes. 

3.1 Modulator 

In Figure 1, the three-channel interplex modulator (Hughes 
Aircraft Company version) is shown. The waveform C(t) is given by 

C(t) = [wg^.t +e(t)] , (1) 

where 

- 1 \ 

0(t) = m2(t)m3(t) 0Q -( -‘'-2 h. (2) 

fhe angle 9(t) is shown in Figure 2. The functions m 2 (t) and m^ represent 
the digital data streams for the two channels, respectively. 
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The total three*channe1 signal is 


s(t) = ^ {CCt) cos o)Qt + S(t) sin u^t} (3) 

where wq is the carrier radian frequency and 

S(t) = v^m^(t) 

Assuming that ni 2 (t), perfectly ±1 V t , then further analysis 

of (1) shows that [2]: 

C(t) = /P7m2(t) g (“sc^‘®o) ^ S (“sc^’^o) 

where 

S (‘"sc^’^o) ■ I '5in(“sc^^®o) ^ I" ®in (“sc^-®o) 

s' (“sc*’®o) ~ ?^in (“sc^'*‘®o) ■ I^in (“sc^’®o) 

Waveforms shown in Figure 3. 

The power latio 



It can be further shown (always assuming thatmi(t), m 2 (t) take on 
only values ±1) that, for this implementation in Mode 1, the bandpass 
hard limiter has no effect on s(t), so that z(t) = s(t). This is 
because 

V(t) = Jp^ m^^(t) + 5in^(“sc^ ^ ^ constant 


Note that, in the presence of distu<^Dances affecting the amplitude of 
either m^(t) or S^h (•) o*' both, this will no longer be true and the 
redistribution of power and the triple cross -modulation term appearance 
will have to be more closely examined. 
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3.2 Demodul ator/Detector 

Figure 4 shows the assumed form of the demodulator/detector 
for channels 2 and 3. A similar loop Is onployed for tracking the 
carrier and demodulating the high data rate channel 1. We assume that 
this loop has no effect on the subcarrier demodulation loop, i.e., 
there Is no phase-tracking error In the carrier demodulation loop. 

We will have to relax this assumption later, when we will examine the 
effect of parameter variations in the high data race channel 1 since 
these variations introduce a phase error (t>^ In the carrier loop. For 
the time being, we concentrate on the effect of parameter variations 
in channels 2 and 3 only. 

Note the presence of the bandpass filter before the subcarrier 
loop. The output xtt) of the bandpass filter can be written as 

x(t) = C^(t) + ngp(t) (8) 

where C^(t) is the bandpass filtered version of C(t) (where only the 
fundamental is retained) and ngp(t) is bandpass noise represented by 

ngpU) = .^|N^(t) cos u)g^t • Ng(t) sin (9) 

In (9), N (t) and N U) are approximately statistically independent 
white Ccompared to the data bandwidth) Gaussian noise processes with 
single-sided spectral density Ng W/Hz and single-sided bandwidth 
Bm < Hz. In the absence of the bandpass filter, such a repre- 

n SC 

sentation would '/e of questionable validity. 

If we expand Sin (•) into a Fourier series and retain only 
the first harmonic, (1) gives [2]: 

= - x/P cos e(t) sin ^ ^t + -JPZ sin 9{t) cos u t 
It V s sc It s sc 

or, using (2) (always with the assumption that m^ (t) = + 1 , i = 2,3^, 
we get 
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( 10 ) 
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C^(t) = m 2 (t) cos 0 Q sin nijCt) sin cos 


(I) t 

sc 


The povMT ratio at the receiver is 



2 

cos 6, 




sin 0, 


cot^ 0, 


as can be seen from (11). Comparing (12) with (7)> we observe that 
the BPF rearranges the power allocation. Define 



1 

2 


1 

2 


(if e„ 

(if 


P 


s 


t 


s 


so that P2r^^3r ” (^ BPF). Then 

(10) is written as 


C^(t) = JZP^^ m^(t) sip +^/ 2 ^iTi 3 (t) cos u^^t 

We assune switching type phase detectors (i.e., square wave reference 
signals) for the coherent demodulation of the signal. Let us assume 
that 


r5(t) = S in [o,^^t - ^ 
r^(t) = Sos [L^^t - ♦] 

where ^{.t) is the phase error process. In (16), we have assumed, 
for simplicity, unit reference gains. 

Expanding (16) into a Fourier series and noticing that only 
the first harmonics 4 /tt sin (w t-<|>) and 4 /tt cos (oj t-(()) correlate 
with the quadrature signals, respectively, we get 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16a) 

(16b) 
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ej(t) " |c,(t) +ngp{t)J • J sin 
e^U) • |c,(t) + njp(t)]. I cos {»spt-»)|LP 


Waveforms e^U) and ej.(t) are the inputs to the chanrel data detectors 
which we assume to be matched filters. The outputs of these integrate- 
and-dump circuits at the end of the kth and Kh signaling interval 
of each channel are given by 

kT +e 

°2 ' h •(. iitI P2(t-(K-') Tj-CjJdt (18a) 

I lk-l)T2+e2 


“3 ' C^(t) Pjft-U-l) T3-C3)dt (18b) 

where P^(t); i = 2,3 is the basic unit power symbol pulse in the Uh data 
stream m^(t) and is defined to be nonzero only in the interval (0,T-)- 
Hence, we have assumed that the data stream is ideally of the form 


o& 

mi(t) = ^ a^ P.(t-nT. -e.) ; i = 2,3 


(19) 


where 


n=-oo 


(T) {^n^U i = 2,3 are independent ±1 sequences with the properties 


1 ] 


(1 

ffl = n 



= 6 

(0 

m n 


E»a a ! 

("i M 

1 = 0, 

all 

n,n and 

i f j 


(20) 


(2) i = 2,3 is the arbitrary pulse epoch in these same modulations. 
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As we shall later see, due to different asynmetrles, the series 
representation of the m^(t) as a function of an^ In (19) Is no longer 
correct. In the following, we shall derive an expression for the prob- 
ability of error as a function of m^(t) rather than an^(t). We shall 
also assume that ^(t) Is constant over the baud of either channel. Sub 
stituting (17) Into (18) and dropping the factor 2*^/u, we get 


^2 '^T *"22 ”32 * ''* 2 ^^^ 

D3 = it>33 cos^ + m 23 sln^ + N3(t) 


(21a) 

(21b) 


where 


m. 


'22 


tn 


m 


m 


j Pp(t-(k-1) T„-C.,)dt 

‘2 (k-l)T 2 +e 2 ^ ^ ^ ^ 


32 


33 


■2 '(k-l)T2+ 2 

tT +£3 

1 < J J 


23 


N2(t) 




. 'V ^3 


(22a) 


n,3(t) P2(t-(k-l) l2-e2)dt (22b) 


m 3 (t) P 3 (t- (fc- 1 ) T 3 - 3 )dt ( 22 c) 


m2(t) P3(t- (e-1) T3-£3)dt (22d) 


N3(t) 


3 (t-l)T3+«3 

cos4> / N^(t) P2(t-(k-l) T2-t.2)dt 

(k-1 )T2'*‘C2 

kl2+t\j 

+ sin>^ / N (t) P,(t - (k-l)T. - eJ [dt (22e) 

(k-l)T2^£2 

- cos^i / N (t) P,(t - (t-1 )T, - e,)dt 


1 

T- 


+ sin^ 


in^ / 




(e,-l)T3>c3 


^(t) P3(t- (4-1)13 -£3)dt 


(22f) 
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Notice that, in the absence of any distortion, (22a) and (22c) yield 

Hk-1),2 ^(t-l).3* ’^spectively. 

From (22e) and (22f), it is easy to show that 



The decision scheme (a zero-level threshold device) outputs 


= sgn{D^} ; i = 2,3 


(24) 


at the ^th signaling interval of either channel. 

Conditioned on ^ 22 * '^32’ *^33* '^23’ random variables 

D 2 and are Gaussian, with mean 


E{D 2 /(rti 2 ; i=2'.3),*} = ^22 cos* - «n 32 sin* (25a) 

(25b) 


E{02/(ifi.2i i~3,2) ,*) ^33 cos* + vP^,. sin* 


and variance = Nq/ 2T^; i = 2,3. 


IS 


H'-nce, the conditional probability of error for each channel 

/ 


P(e2/(fli^2’’"^’^^ = Q 


/ _ /R3 \ 

sgn{rti22}‘^\/R2 *^22 

P(e3/(ffl^.3;i=3,2) ,*) = Q|^9n{tii33) '^33 cos* + ni23 sin*|J 

channel) 


where 

and 


(27) 


(28) 


Ny = T 3 /T 2 (»1) 

is the data rates ratio. (We have arbitrarily assumed that Channel 2 has 
a higher data rate than channel 3) 


and 


Q(x) . J. «pj- 


dx 


(29) 


Finally, sgn(-) stands for the sign function. 
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4.0 AVERAGE ERROR PROBABILITY 


To obtain the average error probability* one has to remove 

the conditioning on the random variables and It Is clear that 

both random variables are data dependent In the sense that apy 

parameter variations In the data stream affect the performance of 

the Costas loop and hence the p.d.f. of 4). Unfortunately, we cannot 

easily obtain the p.d.f. P^(d>); hence, we are facing two possible 

<P 

alternatives: 


(1) Evaluate the 4>-conditioned P(e./<j») of either channel just 
by removing the conditioning on m-j. In the case where the variance 
Is very small (which Implies high SNR in the loop), then P(e-/4)) stands 
as a reasonable estimate of P(e^). The result of such an analysis will 
be a set of curves of P(e^/(p) as a function of SNR, parameterized by 
the mean value of 4> and the degree of asynmetry. 


(2) Expand P(e^/m^ j ,4) into a Taylor series in 4) around m^ 
and retain only the first three terms (a similar approach has been 
taken in [5]). The results are credible if one again assunes small 
a . . Then , 

<p 


dP(e./m,m ) 

P(e^-/m,*) = P(e^./rt,m^) + ^ ^ • ((^ - m4>) 

d2p(e/m,mj 

— ^7^* 


(30) 


so that. 


if P^(^) is 


the marginal p.d.f. of 4, we get 



/p(e./n,4) 

P(e.S,m^) 


P (.t) d4 

d^P(e^./m,m^) 



( 31 ) 


where m s (m^^.m.^); i,j * 2,3;(i ^ j). 

Equation (31) is parameterized by m^ and o^. Removing the 

9 9 

m-conditioning from each of the three terms in (31) will yield P(e-). 
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The second approach takes us one step further in approximating 
P(e-) closer because of the presence of the correcting second term. 
From (26a) and (26b), we have 


dP(e^/m,(|i) 


d<t> 




exp 


z.^) d2^ 

~l\w 


(32) 


where z^(i>); i=2,3 is the argument of Q(*)- But 


dz 


2 _ 


d(» 

dz 


/ ~ ^3 \ 

sgn{m22}-l,^ m22 sin# + 

— = - sgnfmj^} m^j sin<> ~ *”23 


(33a) 

(33b) 


so that 
d 


/S- 


/ - - ^ ^ 
l\^2 "*22 ~ f N~ ”’ 3 '^ 


• exp 


— )-sgn{iti 22 f (34a) 


and 


d 

d^' 


,'an 


P(e3/m,$)= 11133 sin«j» - "*23 

j ( v/r3^33 cos, 

I ) 


' exp 


7 

d"P(e2/(n,,‘') S9n{m22’r ^ 




df 


/2„ ■ ' i 


(35) 


where 


dz. 


si 


j^~- m32 cos^ ~ "if ' sgn{m22l 


(36) 


and 


dz. 


^2 sqn 


(3’) 
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From (35), (36) and (37), it follows that 


d"^P(e2/m,(fr} 
-T 


sgn{ni22} 




2 

Z2-sgn{m22l‘exp |- 

K) H sgn{"’22})‘^"Pr 


^ ^2 


I ^2 
y T( 




exp 


. |l + (2*)^ 


or 


d^P(e«/m,^) sgn{m95,} / 

5 = • 1 vR, m, 

d^ ^ \ ^ ' 


/R, 


^22 "* / N ^32 

V T / , 


I ! ~ ' 1 

^/R2m22 sin^i + m22COS(j>] 

T 


?| /RjlIijjCOS* -/ 


J 


Similarly, from (34a), 



- ■ V 


(38) 


d^P(e3/m,«{.) sgnlm^^} ^ 


L iLv 


(39) 




where 


dz. 


25 - yi<3 «33 5 fn* - ">23 cos* = - -g^ sgnjnijj} 


(40) 


and 


★ 

d(}> 


= Z3 sgn{m33l 


(41) 


The result is 


d‘^P(e3/m,4>) sgn{iii33} 


• *^33 c°5()s + ">23 sin(|ij 

^^17)33 cos 4 i+^/R^ni 23 sin(t.j 3 | 

I 


d4> 

'H#3”33 m 23 C0S(j>j 


( 42 ) 
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Equations (38) and C42) are explicit functions of A and ()>. Once 
these are determined (see next section for m), the second derivatives 
can be evaluated. Note the presence of the factor sgn(*)> as a conse- 
quence of (26a) and (26b). 

Me proceed now to evaluate In (23). Several assumptions 
will be Imposed, and the analysis will advance from the simplest case 
to more complicated cases by relaxing some of the assumptions. 

The parameter variation to be examined is data asymmetry. 

We shall assume that both channels suffer from data asymmetry of 
degree A, which Is presumably the worst case. Define 

A = ^ (length of long pulse - length of short pulse) t (nominal pulse length) 

(43) 

See Figure 5. Let us first assume NRZ-L format for both data streams. 
Parameter a of Figure 5, pertaining to the amount of asymmetry occurring 
at each transition, relates to A as a = A/2. Let us also first calculate 
11122 ^ 33 ' these quantities measure the “self-degradation" 

of each channel due to asymmetry; hence, the epochs play no role 
(we assume perfect sync*'’^)nization for each individual channel although 
the QPSK scheme will most likely be staggered because of the different 
data rates between the two channels), so that no assumption for the 's 
is necessary up to this point. 

The values of ^22 and m^^ pertaining to an individual symbol 


depend on 

the adjacent symbols. Table 1 

sunfvnarizes 

the results. 


Table 

1 • ^^22 ^33 ^ 

(NRZ/NRZ) 



(i)(i)(i) 

Ci)(i)(-i) 

{>1)(1)(1) 

(-1)(1)(-1) 

Prob 

1/8 

1/8 

1/8 

1/8 

^22 * 
■^33 

1 

1 

1 

1 


(i)(-n(i) 

(i)(-u(-i) 

(-1)(-1)(1) 

(-1)(-1)(-1) 

Prob 

1/8 

1/8 

1/8 

1/8 

^ 22 ' 

■"33 

-1 + 2A 

-1 + A 

-1 +A 

-1 
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4.1 High Data Rate Channel 

Having assuned ^ 3 >T 2 » we now proceed to evaluate 1n^2* 

Because Tj>T2t within an interval of length T2 one of the following 
exclusive events can happen: 

3 

= the event of having Channel 3 at the state in the 

interval [Ck-l)T2.kT2] 

= the event of having Channel 3 at the state aj^3=-l in the 
interval [(k-l)T2,kT2] 

E^^ = the event in which Channel 3 switches state in [(k-l)T2»kT2] 
with corresponding probabilities P(E^,), P(E^i), and P(E^,,) such that 

o ' I “ I 5W 

P(eJ^) + P(e:^) + P(E^^) = 1. 

Assigning probabilities to the aforementioned events is a 
rather difficult task, since these events depend on the ratio T2^T3» 
the amount of asymmetry present, the epochs statistics, etc. However, 
things are simplified completely if we make the assumption that 
T3/T2 » 1 . This assumption is justified in our case, where 13/12^10 
(see [ 1 ]). Then, as T 3 /T 2 "^"’ ^sw probability measure zero, 

so that 

P(E+i) + P(E?i) = 1 (44) 

3 3 

In the absence of asymmetry, it would be P(E^-j) = P(E.^) = 1 / 2 . However, 
because of the asymmetry A present, each of the two probabilities equals 
the average percentage of time alloted to each symbol correspondingly, 
i . e . , 


^ ' *1 • ( 45 ) 


(See Table 1 also. ) 



16 


Combining (23b) and (45) and the fact that both p^(t) and 
P 2 (t) are NRZ-L data, we form Table 2: 


Table 2. (m^ 2 ) 


*k3 

+1 

-1 

Prob 

^'* 7 > 

^ 1 -^) 
2 U 2 ^ 

*32 

1 

-1 


We are now in a position to write an expression for the average error 
probability in the high data rate Channel 2: 

7 (' *1) '’(VC.')) + 1 (' -|) P(e2/(1,-1)) 

* (’ *?) * 7 (' -?) '■(V*-'**.-'))] 

* sQ (' P(V(-'t2A.1)) +}('-?) P(e2/(-U2A.-l))] 

* 7[7('*7) '’(VCl.-D) 1 1 -|) P(e 2 /(-l.-l))] (47) 

where the notation 

P{e 2 /a, 6 ] - P ( 62/^22 ^ 3, m 22 = b) (48) 

and the corresponding expression for P(e 2 /(m 22 »'^ 32 )) given by (31) in 
conjunction with (38). (Case I is a special case of II when = 0.) 

We now sunmarij? the conclusion of this section as follows- 
Equation (47) yields the avetcg® hit error rate for the high data rate 
Channel 2. Assumptions used to derive (47) are sunmarized below. 

(1) Phase-error variance is small, so that the truncated 
Taylor series (30) gives reasonably good results. 

(2) Both channels suffer from the same degree of asymmetry A 
as defined in (43). The case of different degrees of asymmetry for 
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each channel is a trivial modification of the above results. Indeed, 
assume that each channel is subjected to as^etry i ~ 2,3 and 
A2/A3. Then each parameter A appearing as an argument of P(e2/m) 
in C47) should be replaced by A2, and each parameter A appearing in 
the multiplicative factors of the form (l+A/2) or (1 - A/2) should 
be replaced by A^. 

(3) Data are NRZ-L for both channels. Cases where one or 
both channels have biphase-L data can be handled similarly. 

(4) Perfect synchronization is assumed for each individual 
channel . 

(5) The ratio T3/T2 »1. This enabled us to induce assunption 
(.44) which led to forming Table 2. 

(6) No assumption about the epoch statistics is necessary. 

4.2 Low Data Rate Channel 

We now proceed to evaluate 0123, and subsequently P(e3). If 
E means "expectation with respect to x," then 

X 

^ ^m33 I ^m23 {*’^V"’33’"’23^}| 

We shall first evaluate P(e3/(m33,m23)) . The analysis will be carried 
out under the assumption that the two channels are synchronized, i.e., 
the epoch difference is identically zero [18]. This means that, 
within a pulse P3(t), there exists exactly N(NeZ^) pulses P2(t) 
(nominally, i.e., in the absence of asymmetry). It is intuitive, how- 
ever, to conclude that, for large rate ratios N, the presence or 
absence of epoch difference is not of great practical import. 

In Figure 6, we have shown a possible sequence of N Channel 2 
pulses (in the absence of asymmetry) within a T3-sec time interval. 

The pulses have been indexed by (i = 0,1 ,. . .N+1 ) according to their 
position within the interval. Note that, because of the existing 
asynmetry, two more pulses (kg and kj^^^) must be take, into account 
since their presence affects the results. We shall argue later that 


★ ^ 

Z denotes the set of positive integers. 
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for N large, their contribution is insignificant, which will result 
in certain simplifications of the final expression. It is the multi- 
plicity of the possible patterns of the N successive pulses, combined 
with the presence of asymmetry, that necessitates the following 
analysis. We make the following assumptions: 

'* y£ [o'n+ 1] ~ event that a pulse occupying the Yth 

position has a +1 or -1 value 

Y»5e [o'n+ 1] “ joint event that the Yth pulse has an 
i sign and the 6^ pulse has a j sign 

k = number of positive pulses in the interval [1,N] 
ti.e., excluding pulses kQ and k|^+-|) 

i - number of transitions from a positive to a 
negative pulse or vice versa in the interval 
[1,N] (i.e., excluding the transitions between 
the boundary pairs (0,1) and (N,N+1) 

*'eff * oumber of effective transitions in the interval 
[1 ,N] Ci*e-« with the boundary transitions 
included, whenever they occur). 


With the above definitions available, it is a question of exhaustive 
search to show that the following conditional relationships hold: 


"eff 


i + 1 


1 + 2 


if 


if 


if 


or 


or 


or 


or 


or 


or 


■"O.N+I 

"0,N+1 

‘'O.N+V 

E^’" 

"0,N+1 

E^’” 

‘'0,N+1 

E“ »'*’ 

‘"0,N-1 

‘"0,N+1 

''0,N+1 

■"O.N+l 


n E 
n E 


HE- 


3.1 
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Note that min k » 0 and max k = N. We now examine the limits £ . and 

jj min 

0^ t as functions of k and Eii. 
max IN 


k * 0 

If ' or 

k » N 


If / or 


e|’" n(l <k<N-l) 
E“|J n(l <k<N-l) 

E|’Jn(2<k<N-l) 
E|’iJ n(N/2 <k<N-2) 


For we have to distinguish between the cases and 

(a) N_ : 


2k - 2 . . . 

. . . if 

E|’Jn(2<k<N/2) 

2(N-k) . . . 

- . . if 

E^’Jl n(N/2< k<N) 
E^;j;;n(llk<N/2) 

2k - 1 . . . 

. . . if or 



E^;Jn{i <k<N/2) 


I e|;J; n(N/2< k<N-l) 

2(N-k) - 1 .... if I or ’ 

( E';J nei/2<k<N-l) 


2k 

2(N-k) - 2 . 
N-2 . . . . 


if E"’" n(0 <k <N/2) 

if E’’" n(N/2 <k<N-2) 

if 
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‘^odd* 

e|:J nz._k<J^-p 

n(N-l)/2<k<N 

e|;‘ m<k<^-i)/2 

E-;J ni<k^N-i)/2 

3.3b 

e|** n (n-i)/ 2 < k < N-1 
Ej’jJ n(N-l)/2 <k <N-1 

E‘;‘ n 0<k<(N-l)/2 
Ej’" n(N-l)/2 <k <N-2 

We now proceed to evaluate F = Pfeym^^)* IF ^ define F(*) to 
mean "F given (•)>" we can apply the law of total probability to get 






p.q=+.- 
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where P|^(A/B) denotes the probability of the event A or that random 
variable A takes the value A, given the event B or that the random 
variable B assumes the value B* In the presence of N Channel 2 pulses 
In seconds. 

The quantities Involved In (50) can be evaluated as follows 
Obviously, the random variable k follows a binomial distribution with 
parameters p = q=l/2, so that 

= it) • (if iir = ^ • (u) 

Because of the Independence of the k. and k- pulses (k. ^k-), we have 

• J I J 

that 


’’n(^o,n+i ~ '*n(^o!n+i) 4 ’ 


The quantity Pfj(*'/CE]’^),k) has been evaluated In Appendix A and is 
given by (A-24), CA-25), and (A-34). 

Observe that the limits t and £ depend on both k and 
(E^*^) through Tables 3.2 and 3.3. The quantities P|^(E^’^/k) have 
been evaluated in Appendix B and are given by (B-2) through (B-4) and 
i.id CB-5) through (B-11). 

Finally, we have to evaluate ’^;k;£). Note that 

the triple CCEQ*^+^)lCEi *^);£] virtually defines through Table 3.1 
With k and as ' parameters, it is easy to check from (22d) that 


m23(^»^»^»£gff) “ ~ (N”k)) T£ A * ^2 


or 


m 23 (N,A,k,£g^^) 


2k - N + A *)i 


eff 


N 


since T 2 /T^ = 1/N. So, simplifying the notation, we have that 


(^O^N+1 ’^1 ,’n**^’^) “ 


2k-N+A*£ 


/"’ 33 ’’" 23 ~ 


eff 


N 


(51) 


(52) 


(53a) 


(53b) 


(54) 
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where is determined as a function of (Eq*^^)» *• 

Table 3.1 

In principle, the problem has been solved. In the following, 
we shall first try to perform certain simplifications that result from 
the formulas and., at a second stage, we shall try to proceed to some 
approximations which make the final result look more tractable without 
sacrificing significant accuracy. 

First, let us consider only the conditioning on k and write 

N*2 

F » P^(0) ‘ F(0) + Pj^(N) . F(N) + P,^(l) • F(l) + P^(N-l) • F(N-l) + £ P„(k) F^ 

k“2 

and, applying (50), (51), (52) and (54), Appendices A and B, and 
Tables 3.2 and 3.3, we find that 


:P.q 


F- ' 

1 


p.q 


F((p,q);(-1,-1);0;0) + F((p,q);(l ,1 ) ,N;0) 

+ F((p,q);(l,-l);l;l) + F((p,q);(-l,l);l;l) 


+ F((p,q);(l,-D;N-l;l) + F((p,q) ;(-l ,1 );N-1 ;l) 


+ (N-2)-F((p,q);(-l,-l);l;2) + (N-2)*F((p,q);(l ,1);N-1;2) 
N-2 


. y 

k-2 


£ (k 


t) 


£=2;even 


max 




£=1 ;odd 


max 


'’(( i-fU) ((^-1 5/2) • f ( (p • ‘) 


£=1 ; odd 




£=2;even 


( 55 ) 
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Note that (55) Is accurate, I.e. , only some 'rangement of terms has 
taken place. 

Equation (55) can be programmed as It Is to yield an exact result. 
If Tables 3.0 (for tg^^) and 3.3 (for t^j^) along nrlth (53) are taken Into 
account. 

4.3 Approximations to F-Equatlon 

(1) We first note that the coefficients of the first eight terms 
within the -t. expression are either 1/(2*'**^) or (N-2)/(2**^^). For 

N = 10, these coefficients are of the order '10 and can therefore be 
neglected. 

(2) We assume that the "edge effects" are negligible, I.e., the 

contribution of In f Is not significant. Indeed, these factors 

will contribute some amount of asymmetry In (0,Tj) during * (1/4 • 3/4 

+ 2/4 • 1/2 + 1/4 • 0)*100% = 50X of the time, and the contributions will, 
on the average, be minimal. In other words, is 

Q 

for ^ -100% of the time 
for ^ *100% of the time 
for ^ *100% of the time 



We can then say that (9* «.+(t+l )*6+(^+2)*l y,6 = t + 1/2, and we 

can use this value independently of Table 3.1. The combination of (1) 
and (2) greatly simplifies F, which new takes the form (also removing the 
factor 1/4 since the conditioning on Eq^n+i removed): 


F 



N-2 


I 



t-2; even 



+ 2 



t=l ; odd 


•( 


N-1-k \ 
(«-l)/2/ 




/ k-1 \/N-l-k\ p* 

Z 

\(i-2)/2/\ t/2 / 

2; even 



( 56 ) 
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where we have defined 

* / 2k - N + A . (t + i)\ 

-) < 57 ) 

and we have also exploited the syranetry between and . (We 

always assune N even.) Furthermore, we can break the sunmatlon In (57) 
Into three disjoint sums (2<_k<N/2), (k*N/2),(N/2<k<N-2) to get from 
(57) and Table 3.3a: 





Observing certain symmetries, (58) is reduced to 
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Finally, if we use the identity 



setting a = k-1, b = N-k-1 and x = £/2, we get 
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it 

Equation C61) is the final expression for ^*^( 63 /^ 33 )* where ^ is 
given by CS7) and P(e 2 /m^*iii 23 ) is given by the combination of (26a), 
(31), and (42). The unconditional Pg(fi 3 ) will be derived by averaging 
F over according to Table 1, i.e., 

^(63) • j P(e3/«33* 1) + j Plej/iSjj* -HA) 




Figure 1. Three-Channel Interplex Modulator 








Figure 4, Subcarrier Demodulator/Detector 
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Figure 5. Data Asymmetry Model 






Figure 6. A Possible Sequence of N Channel-2 Pulses Within One Channel-3 Pulse 
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GLOSSARY OF SYMBOLS 


^(t) « Tracking phase error process 



i.j = 2.3 



{W.Nii]} ■ 

i.j-i.-i \ • 

Y»56[0,n+l]| 


k = 
l = 


^eft ^ 


J> j> = 

min’ max 

F(x) ^ P(e 3 /m 33 ,x) = 


T = 



Mean of ^ 

Variance of ^ 

Data rates ratio 

Error probability of channel i conditioned on 


The event that the ^th pulse equals 1 

The joint event that the yth and 6th pulses equal 
i and j, respectively. Note: Whenever there is 
a possibility of misunderstanding, it is replaced 
by (i.j). 

Number of positive pulses in [l.N] 

Number of transitions from +1 to -1 or vice versa. 
Within the internal [l.N], with the boundary tran- 
sition excluded. 

Number of effective transitions in [l.N], with the 
boundary transitions included, whenever they occur. 

Limits of i 


Probability of error for Channel 3 conditioned on 
m 33 and other events k. 

Number of pairs of adjacent positive pulses in 
[l.N], given in e]’^ and k. 

Number of combinations of m objects taken every n 
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ADDENDUM A 


EVALUATION OF THE CONDITIONAL PROBABILITY DISTRIBUTION 
OF THE RANDOM VARIABLE l GIVEN THE EVENT eJ | j AND k 


In this appendix, we evaluate the probability distribution of 
the random varlpJle £, conditioned on k and Let us first 

define P^U/tL||-J).k) ^ Pa/Cl.j).k): 


(1) If k = 0. then P(0/C-1 .-D.O)* 1 (A-1) 

and every other cond. prob. Is zero. 

§ 

(2) If k=N, then P(0/1 .1 ),N) = 1 (A-2) 

and every other cond. prob. is zero. 

(3) If k = l, then P(l/a.-l),l) = P(l/(-l.l).l) =1 (A-3) 

and P(2/C-l.-l).l) =1. P(t/{1,1),1) =0 (A-4) 

(4) If k = N-l, then P(l/Cl,-l),i1-l) =p(l/(-l.l),N-l) =1 (A-5) 

and Pt2/(l.l).N-l) = l, P(il/(-l,-l),N-l) =0 (A-6) 


Let us assume 2<k< N-2, and let us define the random variable 
T as the number of pairs of adjacent positive pulses (kp,kp^^), given 
and k . Then it is easy to verify that 


2(k-t) 

2(k-t) 

2(k-x) 

2{k-x) 


2 

1 

1 


i lc M2 

K 2 

k - ^ 

K 2 


if 



if 



if 



if 


E 


1,N 


(A-7) 


Let us now consider the set S of all possible pairs of adjacent 
(positive) pulses. Then, 
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of cardinality k>1. The way we select t of these pairs does depend 
on hence, we distinguish the following cases. 

A.1 

In this case, 1, kj^^N. The k pulses lie In an Interval 
of N' = N-2 slots. Then, t of the pairs In S can be selected In 
ways, and the possible selections are all equiprobable. Hence, we 
can write 

P(T/(-l,-l).k) * (•';’) PCEj/H.-D.k) (A-8) 

where the event In CA-8) stands for any particular selection of 
the T pairs. We shall select the following particular pattern: 

^s ~ * ^^2**^3^ '** ^^T*^T+1^/^1 (A-9) 

I.e., Eg stands for the collection of all patterns having the first 
(t+1) positive pulses adjacent, and the rest P=k-(T+1) nonadjacent. 

The above hold for x > 0. If x=0, then 

P(x=0/(-l,-l),k) = P^no adjacent pulses/(-l ,-l ),k] (A-10) 

and will be derived later In this appendix. 

We shall try now to find the probability of the event Eg In 
(A-9). We write 

P{Eg/(-l ,-l ) ,k} = 5(t,N' ,k)*P{any Individual pattern/N' ,k} (A-11) 

where H(t,N',k) Is the number of patterns within the set Eg, and 

P I any Individual pattem/Nik|= r-j-x (A-12) 
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Next, we evaluate H(T.N',k). We note that 

h(t.NM<) - H(T,N',k;k^) + s(T.N'k;k2) +...+ 5(T,N'.k;k^) (A-13) 

where h(t,N' , k;k^), ® 1 »• • • »k^) stands for the number of patterns 

belonging to whose first positive pulse occurs at the position k^. . 

The upper limit k^^ is easily shown to obey the equation 

or 2[k - (t+1)] + (t+ 1) + k^ = N 

k *N'-2k + T + 3 (A-14) 

m 

We can then combine (A-13) and (A-14) into 

N'-2k+T+3 

i(x,N',k) • ^ =(^,N',k;kj) (A-15) 

ki=2 

in the following, we proceed to evaluate 5(x,N' ,k;k^ ) . From 
the setn ' of the problem, it follows that this number equals the 
number ,p' of all the possible patterns containing P = k-(x+l) posi- 
tive nonadjacent pulses within an interval of = N' - [k^. -2 + (x+1)] 

== N' - Ck^ + x) + 1 slots. From (A-14), it follows that tt^. takes on the 
values e [(2k - 2x - 2), CN' - X - 1 )]. An additional constraint should 
be added here, namely, that the first slot within the places should 
not contain a positive pulse (since, otherwise, we would have x+2 
pairs of adjacent positive pulses, contradicting the problem formula- 
tion). An additional constraint is that = 2k-l-N', in the 

case k > N/2. 

One of the following two ways could be used for solving this 
combinatorial problem. The first is to formulate a difference equa- 
tion which the variable cl! should satisfy. Such an equation arises 
by splitting all these patterns enumerated by ^ into the ones con- 
taining a positive pulse at the second slot within the ir-interval, 
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plus the ones which do not. Hence, we can write 


TF ^ ir-1 

with initial conditions 

1 » 


+ I 


*-2 

P-1 


; IT >3 


^ 


. = 0 


(A-16) 


= 1 . 




= 0 
IT 


and try to solve CA-16). 

The other approach is to consider that we have a collection 
of objects of two kinds: 


First kind: 


irP- i.e., a negative pulse followed by a positive 

pulse. This assumes that no t>ro positive pulses 
will be adjacent. There are p of these objects. 


Second kind: ■ ^ y » i.e., a negative pulse which does not accompany 

a positive pulse. There are ir-2p of these objects. 

Hence, the number of all distinguishable permutations containing these 
two objects is 


= (V) 

If we substitute (A-17) into (A-16) and use the identity (JJ) = + (JJl]) 

we indeed verify that (A-17) is the solution of (A-16). 

Combining (A-8), (A-11), (A-12), (A-15), and (A-17), we con- 
clude that 

P(T/(-l,-l),k) = >^2\ ^ \ k-T-1 ) 

(k/ ^-2(k-T-1) 

or, changing variables ((j) * ir^.-k+t+l ) , we have that 
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(k-T\ N-2-k 

• TjjV- "2 
V k / ^®k-l-T 


(A-18) 


The constraints for CA-18) are as in Table 3.3a. Also, we agree that 

( 0 ] 


6 = 1- 


Finally, the desired quantity P(t/(-l,-l),k) is obtained from 
CA-18) and tA-7): 


/ k-1 \ N-2-k 

P(t/(-l,-l),k) = 2 (t/2-l)* ^*even>2 

\ k j ' ' 


(A-19) 


♦=£/ 2-1 


with constraints 


/ 2 ^ k < N/2\ / N/2 < k < N-2 \ / k = N/2 \ 

\2<£<2k/ \2 < £ < 2(N-k)-2/ \2 < £ < N-2/ 


For example, P|^(£=2/(-l ,-l ) ,N-2) = 1/1 


[q^ = 1 , which is correct. 

We now examine P(.t= 0/C-1 ,-l),k). Obviously, P(x=0/(-l ,-l),k) 
if k^N/2. Hence, we examine only the P(t= 0/(-1 ,-l ),k) with k<N/2. 
Observe that, if k = N/2, t> 1, and (A-19) is applicable as it is. 
After some thought, we can see that 


= 0 


P(x=0/(-l,-l),k;2 < k < N/2) = 


.N'+l 

k u rN'+l 

; where 6,, 




(A-20) 


However, we can show that (A-20) is virtually incorporated in (A-19). 
' deed, for s, = 2k (.which is equivalent to t = 0), (A-19) gives 

N-2-k 

I 

P(£=2k/(-l,-l);k;2 < k < N-2) = 

( k 


' * ) 
,k-i; 


(A-21) 


cind to verify the claim, we have to show [combining (A-20), (A-21)] that 
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rr) ■ i^v ) ' T (.^) 

(fr=k-l 


To prove this, we use the identity ([12], problem 2-2): 




Putting r=k-l, r+n = N-k-2 in (A-23), we verify (A-22) directly. 

The final result for the casj is summarized in the 

following (after minor manipulations and using (A-22) again): 


P(0/(-1,-l);0) - 1 
P(t=2/(-l,-l);l) = 1 


P(£/(-l.-l);k) 

where 



i = even l 2; 2 ^ k N-2; £ obeys restrictions of Table 3. 3. a 
Finally, P(£/(-l ,-l );k) = 0, elsewhere. 


A.2 ce;:;) 

This case is exactly symmetrical to A.l, if one substitutes 
the k positive pulses with N-k negative ones. Hence, the result 
can be directly derived from (A-24): 

P(0/(1,1);N) = 1 , P(£=2/(1,1);N-1) = 1 

/k-l\ /N-k-l\ 

P(£/(i,!);k) = — 

\k-2/ 

where i = even ^2; 2 ^ k £ N-2; i obeys restrictions of Table 
Final ly,P(£ /(l,l);k) = 0, elswhere 



(A-22) 

(A-23) 


(A-24) 


(A-25) 

3.3a. 
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Note that, in order to derive CA-25) from CA-24), we have used the 
1 denti t1 es 

/N-2\ /h-K-A 

Vk-2/’ \N-k-f/ \|-l/ 

A. 3 Ce|’|,) : 

These two cases are identical, so we treat the first of 
them. To find PCVU »-i)«k). we use total probability to write 



P(t/(l.-l).k) = P(t/(l,-l);k;k2=2)-P(k2=2/(l,-D.k) 

+ P(il/(l,-l);k;k2>2)*P(k2>2/(l,-l),k (A-26) 


We 


first find Ppj(k 2 * 2 /Cl ,-l ) ,k) and Pnj(k 2 » 2/(1 ,-l),k) as follows 


P(k2=2/(l,-l).k) 


^k2=2;(k-l) positive pulses in [2,N-1]/(1 ,-1 ) 
pj^(k-l) positive pulses in [2,N-1]/(1 ,-l 






(A-27) 


Then 


k2 > 


2/(1, -l);k] = 1 . p[k2=2/(l.-l);k] = 



(A-28) 
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Second, we note that 


P t/(l ,-l);k;k2*2 * P ;(k-l ) positive pulses in C2.N-1] (A-29) 


and 


P t/(l ,-l );k;k 2 > 2 = P ;(k-l ) positive pulses in [2.N-1] (A-30) 


If we substitute CA-27) through (A- 30) into (A-26), a recur- 
sive relationship appears for evaluating P(2./(l,-l).'k), since the iH 
quantity in (A-30) is of the form (A-24) for the CE^’j^) case. To 
proceed further, we split again the expression in (A-29) in the same % 

format as the original CA-26), replacing N-»-N-l and k-*-k-l. T^i s 
can be repeated n steps, and the result will be of the f orm ^ffter 
some manipulation): 



P[ii/(l,-l);k] 



(A-31) 


The r,v. n accounts for the number of recursions needed to bring 
he quantity P|y|EVCE|’j[j,f,).tk-n)] to zero. A question that arises 
now regards the value of n„^^. From Table 3-3a, we see that 

t 2k if 1 < k < N/2 

^max ^ 2{N-k) - 2 if N/2 < k < N-2 

we are referring to the easel 
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Since both N and k are reduced by n at each step, the second 
limitation is always satisfied. Hence, n (i.e., the maximum number 

hiqX 

of permissible steps) will be determined by forcing 


or 


max 


max 


. » - 1 . 2(k-n^^) 


t-1 


- k 


From (A-24), we have that 






k-n 


so that (A-31) canbined with (A-33) yields 


P(V(l,-l);k) 


/N-k-l\ 

(^) 

w 


k--^i / k-n-1 

I 1 

n=l 


kT- 


- li 


or, using identity (A-23), we finally arrive at 


P(ii-V(l,-l);lJ = 1 , P(£=l/(l,-l);N-l) = 1 


(A-32) 


(A-33) 



(A-34) 


voere i ~ odd> 1; 2^k£N-2 and i obeys the restrictions of Table 3. 3. a. 
Finally, ( (2./(1 -1 ) ;k) = 0, elsewhere. An identical result holds for the 



ADDENDUM B 


EVALUATION OF THE CONDITIONAL PROBABILITIES OF 
THE EVENTS e| GIVEN k POSITIVE PULSES IN N SLOTS 

In this appendix, we evaluate the probabilities " 

p((i.j)/kj as follows: 

?((l jj/k) = . p ((l , l);(k-2) positive pulses in [2,N-1] ) 


and, since the two events in the numerator are independent, we get 


P((l.i)/k) 


_ ^ Pf(k-2) positive^^ulses in [2,N-1]) 


or, finally. 



P((l,l)/k) 



(B-2) 


and similarly. 


/N-2) 

P((l,-l)/k) = P((-l.l)/k) = 

(k) 


(B-3) 


and 


P((-l,-l)/k) 



The above were derived assuming 2 ^ k < N-2. Obviously, 


(B-4) 



B2 


P((l.l)/0)» P((l.-l)/0) » P((-l.l)/0) * 0 (B-5) 

P({-1.-l)/N) * P((U-1)/N) * P((-M)/N) = 0 (B-6) 

P(0.1)/l) = 0 (B-7) 

P((-l.-l)/N-l) = 0 (B-8) 

P((-T,'l)/0) = P((1J)/N) = 1 (B-9) 

P[{1,-1)/1) = P((-M)/l) = P[(l,-1)/N-1) = P((-l,l)/N-l] = i (B-10) 


P((-l.-l)/l) = P((1J)/N-1] 


N-2 

IT 


(B-m 



APPENDIX III 

PEAK REGULATOR DESIGN AND BREADBOARD EVALUATION 



APPENDIX A 


PEAK REGULATOR DESIGN AND BREADBOARD EVALUATION 


1.0 INTRODUCTION 

As a part of the development and specification of the PI/KuSP 
bent-pipe interface, Axiomatix undertook to design and evaluate a signal 
waveform peak regulator. Such a regulator permits the largest possible 
peak frequency deviation of the FM transmitter for any waveform shape, 
thus maximizing the mean-square deviation and FN link SNR. 

The purpose for this activity was twofold: 

(1) To demonstrate the simplicity of implementation using readily 
available integrated circuits 

(2) To show that the loop will be stable and perform to expecta- 
tions for all input waveforms. 

2.0 FUNCTIONAL DESIGN 

Figure 1 is a functional diagram for the peak regulator loor. 

The inp t, V , is first passed through a two-pole LPF having a 3 dB fre- 
quency of 4.5 IWz. This LPF is redly not a part of the loop proper but 
is necessary to simulate the output bandwidth characteristics of the PI 
receiver, following the LPF is a gain-controllable amplifier (GCA) whose 
gain, K, is voltage controllable by a reference voltage, The nominal 
open-loop gain of this amplifier is +16 dB and is variable by the refer- 
ence from about +36 dB to -34 dB. Prior to the GCA is a 21 dB attenuator 
and following the GCA is a +20 dB fixed gain amplifier, so that the over- 
all nominal input/output gain is +15 dB (a voltage gain of 5.6). Since 
the regulated output voltage is specified to be ±2V p-p, the nominal 
input voltage is 357 mV p-p. 

Output of the regulator, V^, is input to an absolute value cir- 
cuit which is mechanized in the form of a half-wave linear rectifier. 

(Note: The earlier functional design showed a full-wave linear rectifier 

but, because of waveform symmetry, the half-wave circuit was judged to be 
sufficient and simpler to implement.) A gated reset peak detector following 













A-3 


the rectifier functions to store or hold the largest voltage peak observed 
over a 0.99 ms time period. At the end of each 1 ms sample period, the 
peak value is transferred to a sample-and-hold amplifier where it is 
retained while the peak detector memory (capacitor) is reset and the next 
peak detection operation takes place. 

Designating the sample-and-hold output by the symbol V, the loop 
instantaneous error, e, is formed by subtracting a IV reference from V. 
Thus, e = V-1. An LPF follows (time constant = 1 sec) which averages 
over a large number of the peak detected error values that may change 
each 1 ms. The LPF output is denoted by e. Following the LPF is a 36 dB 
gain an^lifier which produces at its output 61e subject to a maximum con- 
straint of ±0.«V. This constraint is imposed to keep the reference volt- 
age to the gain-controllable amplifier within a linear operating range, 
precluding the possibility of an unstable loop condition. The amplifier 
voltage reference, Vj^, is fonned by offsetting the error 61 e by -3.6V. 

The original specification of the loop regulation performance 
was to be ±1% on the output for a ±20 dB input variation. Thus, the 
equation of regulation would be 



5.6 V.p X 


ioo(v^^-n 

10 op 


where V^p and V^.p are, respectively, the regulator waveform peak output 
and input voltages. As will be seen from the actual circuit performance 
discussed in the following section, this capability was achieved over 
only a 16.5 dB range. 

The peak error sampling rate was chosen as 1 kHz, and the peak 
detector averaging time per sample is 0.99 ms. Thus, for any of the var- 
ious waveform shapes considered, where the lowest subcarrier (sinusoid) 
frequency is expected to be about 30 kHz and random noise occupies the 
full 4.5 MHz bandwidth, each peak sample should be very close to the true 
peak value of the waveform. Averaging over a thousand or so error samples 
also provides the loop with a reasonably rapid response to dynamic input 
level changes, but is sufficiently long to obviate response to very short 
signal transients or the possibility of an occasional impulse noise burst 
due to EMI or other sources. 
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3.0 REGULATOR CIRCUIT DESIGN AND PERFORMANCE 

Figure 2 is a circuit diagram for the Axiomatix peak regulator 
breadboard, and Figure 3 shows the companion timing logic circuits. As 
most of the regulator design is based upon operational amplifier configu- 
rations, a correspondence between the circuit diagram and the functional 
diagram of Figure 1 may be easily established. The amplifiers having the 
AD prefix are Analog Device types, and AD583 is the sample/hold amplifier. 
Block AH0152 is a FET switch used to discharge the peak detector capacitor 
(100 pf). Clamp for the error voltage amplifier output is provided by 
the pair of reversed IN457 duo-diode groups. Timing waveforms are pro- 
duced by monostable multivibrators. 

The RCA CA3002 was originally selected because of its large gain 
control capability— up to 70 dB for a 1.5V control voltage differential 
range. When the amplifier was tested to ascertain all of its operating 
characterisitics, it was discovered that the maximum input voltage had to 
be limited to 150 mv p-p. Above this value, virtually independent of the 
gain control bias, the amplifier output exhibited a voltage saturation 
(compression) condition. The result was that the amplifier could not be 
driven to the levels intended in the original design. This, therefore, 
necessittted an adjustment of the intended operating point and resulted 
In a regulating range of 10 dB below the nominal input and 6.5 dB above 
the nominal input (rather than ±20 dB). It was decided, however, that 
this would not compromise the prime reason for the breadboard evaluation 
of demonstrating excellent pcak-to-peak regulation as a function of a 
variety of complex waveforms plus 4.5 MHz lowpass noise. 

The initial regulator breadboard was constructed with a full-wave 
linear rectifier (see Figure 34 of the Final Report for 1978 Contract 
NAS 9-15240D for the circuit). This circuit performed very well as a 
function of input frequencies up to 2.5 MHz. Above 2.5 MHz, however, 
unsymmetrical phase shifts ippaared between the waveforms produced by 
the individual half-wa... rectifier outputs. Further, the overall fre- 
quency response was limited to about 3 MHz due to the inability of the 
AD507 amplifiers to preserve the harmonic structure of the full-wave 
waveform above this frequency. Rather than redesign the full -wave rec- 
tifier using different (wider bandwidth) amplifiers, it was decided to 




too 


2AK 



tor Circuit Diagram 










Figure 3. Peak Regulator Timing Logic Diagram 
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use a linear half-wave circuit Instead, as all the contemplated Input 
waveforms are amplitude symmetrical. This change eliminated two A0507 
amplifiers, with the result that the half-wave rectifier response was 
found to be adequate up to 4.2 MHz utilizing IN914 diodes. 

All the reii'.a1n1ng circuits performed essentially as expected. 

A 470 n resistor had to be Installed between the FD300 diode and 100 pf 
capacitor In the peak detector circuit In order to suppress an undesirable 
switch-off transient from the AH0152 FET discharge switch. 

4.0 DEMONSTRATION AND TEST RESULTS 

A measurement of the loop regulation capability was made In the 
following manner. Figure 4 shows the configuration used to obtain the 
necessary data. The approach taken was necessitated because of limited 
equipment but proved to work reasonably well. 

The regulation error was measured In a differential fashion. As 
shown In Figure 4, an oscilloscope was employed to measure the AV^ as a 
function of V^. The measurement was calibrated for V^. = 178 mV-p and 
a 1 v-p at a frequency of 500 kHz. (A 500 kHz frequency was selected 
for this measurement so as to minimize any possible phase shift through 
the regulating circuits. Some perceptible phase shift at f = 2 MHz had 
been observed through the CA3002 amplifier as a function of its gain.) 

The scope, operating in its differential input mode (Y^ - Yg), was cali- 
brated on the ^2 Input, while the Y^ variable attenuator was adjusted to 
produce an amplitude null on the scope screen. Voltage, V, was held con- 
stant while was changed using the operational adjustable attenuator 
and measured using a high-frequency RMS-to-DC converter in conjunction 
with a DVM. The result is that regulation error, AVq, as a function of 
, could be read directly on the scope. The only problem encountered 
was that the wideband noise generated in the regulation loop amplifier 
circuits made it difficult to measure error voltages below 5 mV. 

Table 1 summarizes the regulation measurements. As can be seen, 
the regulation range of -10 dB to +6.5 dB about the nominal (calibration) 
point was achieved. 

The throughput frequency response of the regulator was measured 
exclusive of the LPF, and the 3 dB frequency was found to be about 5.2 MHz, 
the main contributor to the roll-off being attributed to the +20 dB output 




Figure 4. Loop Regulation Measurement Configuration 
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Table 1. Measured Regulator Performance 


*1 

*0 

% Error 

25 

860 

4.0 

50 

970 

3.0 

75 

992 

0.8 

100 

996 

0.4 

178 

1000 

0.0 

300 

1002 

0.2 

400 

1006 

0.6 

500 

1009 

0.9 

600 

1011 

1.1 

700 

1018 

1.8 

800 

1028 

2.8 

900 

1050 

5.0 

1000 

1100 

10.0 


Exceeds Regulator Specification 


Calibration 


Exceeds Regulator Specification 


CA3002 Amplifier in Compression 


All voltages are in mV -peak 
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amplifier. With the input LPF connected, the throughput 3 dB frequency 
was 3.9 MHz. 

The remainder of the measurements made on the regulator were 
rather qualitative and consisted of mixtures of signal waveforms and 
noise. A 1.024 fWz sinewave (representative of a subcarrier), a 200 kHz 
square wave (representative of NRZ data), and random noise were combined 
at varying levels and the regulator input and output observed on the 
scope. The performance was judged to be as expected. It was estimated 
that the regulator was able to hold regulation on the 4.5 MHz noise peaks 
at about 2.75a. 

The regulator was never observed to exhibit instability, and was 
able to accommodate input on/off step transients. 

5.0 CONCLUSIONS 

Except for the somewhat limited regulator dynamic range of 16.5 dB, 
all other aspects of the peak regulator design and breadboard evaluation 
met the objectives set for the effort. If, in the future, it is desired 
to increase the dynamic range to the original design goal of 40 dB, this 
may be done by replacing the CA3002 amplifier by a FET or diode current 
controlled attenuator configuration. 
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DETERMINATION OF DON TRANSPONDER ACQUISITION SWEEP RATE 


CORRECTION NOTICE: THE LAST PAGE HAS BEEN CORRECTED AND REVISED 

ORIGINAL: 


6/26/79 


CORRECTED: 7/17/79 



THE bAblC PROBLEM 


Axiomafix 


DSN TRANSPONDERS WILL BE SWEEP ACQUIRED BY THE PI UNDER CONDITIONS OF MODERATELY 
STRONG SIGNAL LEVEL 

AT MODERATELY STRONG SIGNAL LEVELS, THE DSN TRANSPONDER PLL IS CHARACTERIZED BY A 
TRACKING BANDWITH (26^) BETWEEN 150-225 Hz AND A DAMPING FACTOR (^) AS LARGE AS THREE 

THE SWEPT PLL MAXIMUM SWEEP RATE LIMIT FOUND IN THE LITERATURE (VITERBI, ET AL) IS 
FOR A DAMPING FACTOR 0F\/2/2 « 0.707. THIS RESULT IS THEREFORE NOT APPLICABLE TO 


THE DSN TRANSPONDER 


PRO^.l-; :1 SOLUTION 


^5 Axiomatix 


• OBTAIN THE ABSOLUTE MAXIMUM SWEEP RATE FOR THE DSN TRANSPONDER PLL TRANSFER FUNCTION 

• THE PROCEDURE FOR A SPECIFIC TRANSFER FUNCTION IS TO FIND THE VALUE OF SWEEP FREQUENCY 

* A^/27tTHAT results IN A SADDLE-POINT B AND \ ‘ “ PLL NATURAL FREQUENCY, 

RAD/SEC) 





nSN TRANSPONDER CHARACTERISTICS 


PLL TWO-SIDED THRESHOLD BW ( 28 ^ 0 ) 

PLL THRESHOLD DAMPING FACTOR 
MAXIMUM SUPPRESSION FACTOR (tt^) 

PLL STRONG SIGNAL BW 

PLL STRONG SIGNAL DAMPING FACTOR ) 

^max 

STRONG SIGNAL SUPPRESSION FACTOR 
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/! 


NOMINAL TOLERANCE MINIMUM 

18 Hz 13 Hz 

0.707 

0.0547 


225 Hz 

163 Hz 

3.03 

3.03 (By Definition) 

1.0 

1,0 


t 




ACQUISITION SWEEP RA^E (Hz/SEC.) 


DSN TRANSPONDER ACQUISITION SWEEP RATE 
VERSUS 

LIMITER SUPPRESSION 


Axiomatix 




LIMITER SUPPRESSION FACTOR 
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MINIMUM DSN PAYLOAD AND ORBITER LINK CONDITION^ 
FOR 10 NMI RANGE 
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PI_ 

TRANSMITTER POWER 



37 dBm 

_27 dBm 

4 dBm 

PAYLOAD RECEIVED SIGNAL LEVEL 

-92 dBm 

-102 dBm 

-125 dfim 

DSN TRANSPONDER NEGATIVE TOLERANCE ALLOWANCE 

-4 dB 

-4 dB 

-4 dB 

LIMITER SUPPRESSION FACTOR 

1.0 

1.0 

0.79 


• FOR VERY NEAR-RANGE, THE LIMITER SUPPRESSION FACTOR IS 1.0 FOR ALL PI 
TRANSMITTER POWER LEVELS 


:> 

f 
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PI TRANSMITTER MAXIMUM NOMINAL SWEEP RATES 
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• 20% BACKOFF FROM ABSOLUTE MAXIMUM PHASE PLANE LIMIT IS ASSUMED 

• +30% SWEEP RATE TOLERANCE ON PI ASSUMED 


PI TRANSMITTER POWER 



37 dBtn 

27 dBm 

4 dBfn 

Wj_Q = 18 Hz 
(NOMINAL) 

470 Hz/s 

470 Hz/s 

365 Hz/s 

Wj_Q * 13 Hz 

250 Hz/s 

250 Hz/s 

190 Hz/s 

(TOLERANCE MINIMUM) 

i 


M 


SWEEP RATE 
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APPENDIX V 

MAXIMUM SWEEP FREQUENCY FOR PHASE>LOCKEO LOOPS 
WITH LARGE DAMPING FAQORS 

By 

Marvin K. Simon 



1.0 


INTRODUCTION 


The classic paper by Vlterbi [1] Investigates the swept acquisition 
capability of a phase lock loop (PLL) using the phase plane approach. 

Therein Is was show) that. In the absence of noise, a second-order loop with 
damping factor c « 0.707 could always* He brought to a locked condition as 
long as the sweep rate (rad/sec) was less than half the square of the loop 

V 2 2 

natural frequency p. For values of 0 /?< Aq< $ , the loop could not be 

guranteed to lock. This Is equivalent tc saying that some of the phase plane 

trajectories corresponding to Initial conditions In certain regions of the 

plane would not approach the equilibrium (lock) point In the steady state but, 

rather, go on to a stable limit cycle. If the Initial conditions were such 

as to bring the loop Into lock, however. It would continue to track. For 
2 

Aq> 6 , the loop would not lock at all. Moreover, If It ever reached the 
lock point. It would not continue to track. 

Extension of these results to other values of loop damping was con- 
sidered by Dye [2] and Tauseworthe [3]. In particular. Dye showed that, for 

2 2 

5=0.5, lock-on is guaranteed for lu < 0.387 0 , marginal for 0.387 0 < 

2 2 ^ 

Aq< 0 , and not possible for Ag> 0 . Tauseworthe verified these results and 
generalized them by investigating the maximum guaranteed lock-on capability 
for a range of damping factors between 0.25 and 1.6. 

The purpose of this report is twofold. First, we wish to extend 
Tauseworthe' s results to larger damping factors. Second, we shall consider 
the situation where the PLL is preceded by a narrowband bandpass limiter (BPL) 
which is typical of the deep-space network (DSN) transponder receiver. Actu- 
ally, these two generalizations will not be performed independently. Rather, 
the increase in damping factor will come about because of the increase in 
SNR in the BPL input bandwidth, thus raising the limiter suppression a from 
its threshold value oq toward its strong signal value. The damping factor 
(relative to its threshold value Cg, assumed here to be 0.707) is directly 
proportional to the square root of the ratio of a to og. 

The approach taken in this report will once again make use of the 
phase plane. A simple method for obtaining phase plane trajectories on a 
digital computer was discussed by the author in [4]. The approach used 
there will prove useful in obtaining the desired results described above. 


★ 

The word "always" as used herein implies an initial frequency offset 
between the input and loop's VCO reference of any magnitude and a sign in the 
direction of the sweep. 
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2.0 THE LOOP EQUATION 

A second-order PLL with perfect Integrating loop filter 


F(s) 


^ ^ V 


Is governed by the differential equation of operation [3] 


*0 ■ ♦ * ’z 

dt 


( 1 ) 


(2) 


where ^ Is the loop phase error. The loop damping c can be related to B 
and t 2 by 




s 


Bt2 

1" 


Thus, defining the normalized time variable t = 2cBt, (2) can be rewritten 
in the form 



2 2 

Further, letting v = d(j/dT, and equivalently, d ^/dx = vdv/d<|), (3) can be 
put in a form desirable for using the phase plane approach of [4], namely. 


dv 

3* 


+ cos ^ 


-4- - sin 

4cS 


(4) 


3.0 THE PHASE PLANE TRAJECTORY APPROACH 

In [4], it was shown that the parameter arc length defined by 

ds^ = dv^ + d4>^ (5) 


was a desirable independent parameter for constructing a plot of versus 
V (a phase plane) on a digital computer. Thus, our goal (as in [4]) is to 
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find a pair of coupled first-order differential equations for ^ and v, 
i.e. , 

and solve these iteratively. To see how this is done, we first observe 
from (5) that 



Substituting (4) into (7) and (8) gives the desired set of differential 
equations in the form of (6), namely. 
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4.0 DETERMINATION OF THE SADDLE POINT PHASE PLANE TRAJECTORY 

The phase planes originally determined by Vlterbi [1] contained 
a trajectory (herein referred to as the saddle point trajectory) which 
passed through the saddle point I * 0, ♦ » ir - sin" (Aq/P ). The sig- 
nificance of this trajectory depends on the value of loop damping relative 

2 

to the normalized sweep rate Aq/p . For example. Figure 1 Illustrates a 
phase plane plot for c = 0.707 and 0 < Aq < p^/2. The saddle point Is 
denoted by point B and the saddle point trajectory corresponds to the 
trajectory passing through points A, B and C. In particular, we note 
that the value of | at point A = -ir) herein denoted by Is greater 

than the value of ^ at point C =ir) herein denoted by Thus, any^ 

phase plane trajectory with ^(0) < 0 (^(0) denotes the Initial condition 
on ^(t)l will stay below this saddle point trajectory and eventually 
reach the equilibrium point D = sin" (Aq/p^), = O). This Is con- 

sistent with our previous statement that the loop will "always" reach lock 
2 

when A«< p /2. 

u 2 2 

If c = 0.707 and p /2 < Aq< p , the saddle point trajectory will 

be such that either or point A disappears entirely. In the lat- 

ter case (see Figure 2), for example, we observe that there exists a cor- 
ridor in the phase plane bounded by the trajectories E - E' and F - F' 
or G - G' and F - F' , for which any trajectory originating in this band 
will be whipped past the equilibrium point and on to a stable limit cycle. 
Thus, despite the fact that Aq > 0 and ^(0) < 0, the loop will not come to 
lock if ^(0) is of such a magnitude as to fall into this critical corridor. 

From these examples, it is easy to conclude that the critical 
value of Aq, namely, pv2, for which all trajectories with ^(0) < 0 reach 
the equilibrium point, corresponds to a saddle point trajectory for which 
^A ~ Figure 3). This, then, is our task, namely, for each value 

of loop damping c, we search for the largest value of Aq/P^ which gives a 
saddle point trajectory with = Jq. This value of Aq, namely, 
is then the maximum sweep rate (in rad/sec) for which the loop is guaran- 
teed to acquire. 

Using the pair of coupled differential equations of (9), we start 
at (in the neighborhood of) point B (whose coordinates are known, namely, 

V = 0, $ = TT - sin" (Aq/p )) and numerically integrate forward (ds 



A Typical Phase Plane Trajectory Showing Solutions 
to the Loop Equation; ? = 0.707, Aq<6'^/2 
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These trajectories 
are in corridor and 
don't achieve lock 


Figure 2. A ' ;ical Phase Plane Trajectory Showing Solution. 

to the Loop Equation 

95 
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Figure 3. Limiting Trajectory for which Loci with 4t(-ir) 
< Converge and |(-tt) > Diverge 
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positive) to point C (^ « v) and backward (ds negative) to point A * 9 ). 
2 

The value of Aq/ 8 at which Is then the desired result. 

5.0 DETERMINATION OF THE MAXIMUM SHEEP RATE IN TERMS OF 
THE DESI'^N POINT LOOP PARAMETERS 

Letting and 6 q, respectively, denote the design point (thresh* 
old) values of loop damping and natural frequency, then, for any operating 
point SNR in the BPL input bandwidth, the corresponding values are given 
by 


B = 



( 10 , 


In (10), Oq is the threshold limiter suppression factor and o is the cor- 
responding suppression factor at the operating SNP. Substituting (10) 
into (9) enables one to determine the critical saddle po^nt trajectory as 
a function of <* for fixed values of Oq, 6q and Cq, In particular, for 
ag = 0.0547 a.id Sg = 0.707, Figure 4 is a plot of 

versus « for values of $g corresponding to two-sided threshold loop band- 
widths Wj^g of 13 and 18 Hz. The relationship between 6g and W^^g is given 
by 


6 


0 



( 11 ) 


Table 1 tabulates some of these values along with the corresponding values 
of B and c as computed from (10) • 



800 



(D0S/ZH) 
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Table 1. PLL Parameters as a Function of Limiter Suppression 


a 

6 

C 

"l 

J 

^max 

0.1000 

16.5719 

0.9559 

20.1773 

0.6140 

26.8371 

0.2000 

23.4363 

1.3519 

36.0213 

0.7320 

63.9894 

0.3000 

28.7034 

1.6557 

51.8653 

0.7890 

103.4583 

0.4000 

33.1439 

1.9119 

67.7093 

0.8260 

144.4132 

0.5000 

37.0560 

2.1375 

83.5533 

0.8510 

185.9801 

0.6000 

40.5928 

2.3415 

99.3973 

0.8690 

227.8967 

0.7000 

43.8452 

2.5292 

115.2413 

0.8840 

270.4689 

0.8000 

46.8725 

2.7038 

131.0853 

0.8890 

310.8556 

0.9000 

49.7158 

2.8678 

146.9293 

0.9040 

355.6132 

1.0000 

52.4051 

3.0229 

162.7733 

0.9120 

398.6225 

r. = 0.707 
0 

“lo = 

13.0000 

3q = 12.2565 

“0 = 

0.0547 
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APPENDIX VI 


THE THEORY OF MEAN-SQUARE PHASE NOISE PERFORMANCE 
OF ONE/TWO-WAY COHERENT COMMUNICATION LINKS 


1.0 INTRODUCTION 

In assessing the performance of one/two way coherent communication 
links, one must consider the many sources of random fluctuations present 
in the overall system and their individual degrading effects. While it 
is most common to consider only those degradations due to the additive 
thermal noise on the link(s), a complete characterization of system per- 
formance should also include the effects of such phase noise sources as 
transmitter and receiver local and reference oscillators, mixer (multi- 
plier) logic noise, /W-to-PM conversion due to link nonlinearities, and 
possible vibration effects. All of these phase noise components degrade 
carrier tracking loop performance and those which are not tracked by the 
loop(s) have an additional harmful effect on overall system error prob- 
ability performance. 

Regardless of the analysis technique employed, the first step 
in assessing tracking loop and bit error performance degradations due to 
phase noise effects is to introduce appropriate f '•hematical models for 
the carrier (discrete or suppressed) tracking Ic a) which include the 
various sources of phase noise referred to above. These models ordinar- 
ily take the form of equivalent linear loop block diagrams with phase 
noise sources described by power spectral densities characterized by 
inverse power law behavior as a function of frequency away from some 
nominal value. A complete discussion of the various phase and frequency 
instabilities present in precision frequency sources, along with their 
characterization, is given in [1]. For our purposes, it is convenient 
(as was done in [2]) to represent the equivalent two-sided lowpass power 
spectral density S^(o)) of the various phase noise instabilities by the 
(inverse) power series* 

We shall ignore any discrete frequency components in the 
spectrum such as those caused by spurs or vibration effects. 
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k»0 

The |u|.^ terms corresponds to white phase noise, |o>p^ corresponds to 
flicker phase noise, |otr^ corresponds to white frequency noise, while 
|u)|” corresponds to flicker frequency noise. 

In the past, the calculation of mean-square phase error in the 
presence of phase noise has been based on an equivalent linear loop model 
for a one-way link, with specific closed-form results obtained only for 
the case of 0.707 loop damping (Eq. (11) of [2]). While the abstract and 
introduction of [2] imply that these theoretical results can be applied 
to assess the tracking loop and bit error probability performances of 
various Shuttle S-band links (e.g., payload/Shuttle and TDRSS/Shuttle), 
in reality, this is not true since these links are indeed two-way coher- 
ent communication systens. In fact, the theory developed in the first 
part of [2] is never used in assessing the performance of these links; 
rather, estimates of the actual power spectral densities of the many 
phase noise sources combined with piecewise linear approximations of the 
various loop transfer functions allow for numerical integration computa- 
tion of mean- square phase error. 

In this report, we begin by generalizing the specific one-way 
results given in [2] to the case of arbitrary loop damping. Following 
this, we present a theory for computing mean-square phase error due to 
phase noise sources in two-way coherent communication links. Here again, 
the two loops are allowed to have arbitrary damping factors and natural 
f’^equencies. Several different two-way transfer functions characteristic 
of payload/Shuttle links (e.g.. Inertial Upper Stage (JUS)/Payload Inter- 
rogator (PI)), are considered and, in each case, expressions for the 
mean- ""’are phase error compone.ts corresponding to each term in the 
power spectral density oi (1) are derived. These results, . in combined 
with ’•imi''^r >jsults available tor mean-square phase error due to addi- 
tive thermal noise [3], then allow for assessment of total mean- square 
phase error performance and the co> responding degradation in bit erro** 
nrobal ility performance. 
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2.0 SYSTEM MODEL 

Consider the simple model of a two-way coherent communication 
system Illustrated In Figure 1. As such, the reference system might rep- 
resent the PI aboard the Shuttle while the vehicle transponder might rep- 
resent the lUS/SGLS transponder.* When the switch Is In position 1, the 
system of Figure 1 operates In a two-way coherent mode, I.e., the uplink 
signal to the vehicle Is coherently demodulated and a frequency- translated 
version of It Is retransmitted on the downlink, which Is then tracked by 
the reference tracking loop. When the switch Is In position 2, a free- 
running auxiliary oscillator aboard the vehicle provides the downlink 
carrier and, as such, the system operates as two Independent one-way 
1 inks. 

An equivalent linear baseband model of the system, including 
various phase noise sources, is illustrated In Figure 2. Using simple 
control system theory, the block diagram of Figure 2 can be further sim- 
plified to Figure 3. From Figure 3, we Immediately observe that the 
phase error (j) in the reference system is given in terms of the various 
phase noise sources by ** 

*^0 6i(s) G2(s) t|».j.j((s) 

+ [1 - H^(s)][l - H|j(s)] G^(s) G2 (s) ^^(s) 

- Cl - H'(s)] 6^(s) G2 (s) 

- [1 - H^(s)] 62 (s) 4<r(s) 

- [1 - H'(s)] ^^(s) (2) 

Assuming for simplicity that the IF amplifiers in the reference system 
are identical ideal "brickwall" bandpass filters, with 6^(s) and G2(s) 
their equivalent lowpass transfer functions, i.e.. 


In reality, the TCXO-1 inear phase modulator combination is 
only a part of a more general-purpose frequency synthesizer aboard the 
Shuttle PI which contains other internal loops and phase noise sources. 
For the purpose of Sections 2.0 through 5.0, we consider only the direct 
fuedthrough paths of the TCXO to the transmit and receive output termin- 
als of the synthesizer. Later on, we generalize the model to include the 
other internal frequency synthesizer oscillator sources and other TCXO 
paths. 

We assume here that the switch iu the vehicle system is in 

position 1. 




Figure 1. A Simple Block Diagram of a Two-Way Coherent Communication System 








Figure 2. An Equivalent Linear Baseband Model of the System 


Reference System 


ing Oscillator Phase Noise Sources 
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!“l 1 “h 

el sewhere 


then» using the notation 

♦(s) « G(s) *(s) , 

Equation ( simplifies* to 

♦g(s) » 6^Hy(s) [1 - Hf^(s)] ♦y.(s) 

+ Gf [1 - Hy(s)][l - Hr(s)] .jiyCs) 

- [1 - H^(s)][$rx(s) + ^r(s) + <l»^(s)] 


The corresponding mean-square phase noise Is obtained from 


j' 


• 5^ / *e<*> “i* 


t 


where the overbar denotes statistical expectation. Letting 

S^(s) 4 

denote the power spectral density (in Laplace transform notation) 
phase noise process iii, then (6) is evaluated as 

H,(s) H„(-s)|l-H,(s)][l-H„(-..)] Sj^^(s) ds 


* W Q- i - "r<-*] * V<^'] 


We further assume that 1« = 1 - Hj^(s), where 


. , A 

Hg(s) = s 'l^' F (s) ’ yivi^a-n^wi 

R R 


respect te !p(s), 


( 3 ) 

(4: 


( 5 ) 

( 6 ) 
( 7 ) 

of the 

s) ds 
Is 18) 



or letting s * j«. 


+ ^/ G4r^ll-H„(o«){^ |i-H|j(j«)|^ s* (*>) du 

•CD 

^ h C [*♦«,<“> * %<“> * %/“'] <*> 

From (9), we observe that three different types of integrals 
occur; the first two involve two-way comnuni cation (i.e., two loops in 
cascade) and the third* one-way communication (i.e., a single loop). 

In the following section, we derive closed-form results for ttose terms 
when S^{(d) is given by (1) and the loop transfer functions correspond to 
second-order loops with arbitrary loop danq>ings and natural frequencies. 

3.0 CALCULATION OF MEAN SQUARE PHASE NOISE 
FOR ONE-WAY COHERENT COMMUNICATION LINKS 

A second-order carrier-tracking loop has the in-band transfer 
function [4] 



where r = 4c is the loop damping factor (c is the loop damping) and 
is the loop natural frequency in rad/s. The corresponding out-of-band 
transfer function is clearly 



Finally, the magnitude-squared versions of (10) and (11) are respectively 
given by 
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|H(W|^ 



|1 - 


b> 




Now consider the evaluation of 


’hi I’ - *• 

*00 

where jl - H{j«>)|^ Is given by (12), S.{«) by (1), and 

i S*(“) 1 1»1 i “H 

•>iVw) ~ ) 

^ 0 ; elsewheri. 

Substituting (1) combined with (14) into (13) yields 


( 12 ) 


(13) 


(14) 


2 

0 


1-way 


4 



1-way 


(15) 


where 


2 

‘’k 


1-way 


1 /"" |1 - H(W|2 dw 



(16) 


For simplicity of the resulting expressions and with little 
loss in accuracy, we shall, wherever possible, allow the upper limit of 
the integral in (16) to be infinity rather than Clearly, this is 
possible (i.e., the integral converges) for k * 2,3 and is justified by 
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the already assumed fact that We begin by evaluating (16) for 

these two values of k before going (mi to the more complicated cases cor- 
responding to k ■ 0»1. 

Integrals of the form required in (16) can be found In both 
their definite and Indefinite forms In a standard table of Integrals [5]. 
With regard to [5], we shall herein adopt the notation [5, p.xxx, Eq.(yy)] 
to Identify equation (yy) on page xxx of this reference. 

a. (k»3) 

From "[5, p. 67. Eq. (2.161-2)]. 


I- 


1 5 — 

X + ax + b 


dx = 




; h‘ > 0 


2 

2q sina 


tan 


•1 .2 < 0 

\ q^ sino / 


(17) 


where 




4b ; f = i 9 


1/4 

q = b ; cosa = - 


Zt/E 


a+h 

2 


(18) 


Letting 


a - (r-2) - 01 ^ R j b = 


(19) 


we have 




/r^- 4 ; f = ^ R - /r^ 


(». 


; 9 = 




q = (i»^ ; coso = 


R 

I 


sina 


/4-R^ 


( 20 ) 
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and 



Evaluating (21) between zero and infinity gives the definite integral 



Finally, letting ® in (16), as previously agreed upon, we immediately 
have from (22) 



/ 
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From C5. p. 67, Eq. (2.161-3)], 

^ fp-.f f-^ ; (24) 

J x^xHb " J x'^+g " J x^+f 

where f, g, and h are defined as in (18). Since b * > 0 in our case, 

then h < a, and both f and g are positive. Hence 



Since (24) results merely from factorization of the denominator of the 

integrand and a partial -fraction expansion, we can apply a similar fac- 
2 

torization for h <0. To obtain a closed-form solution analogous to 
(25), we proceed as follows: 

For h^ < 0, let h^ = -h'^, or 


h = jh* = jv4b-a^ 


and 


r 2 * 


9 = ^ 


Expressing f and g in polar form. 


f = ; 9 = |9|e^® 


where 



(26) 


(27) 


(28) 


T! • O 


(29) 



( 30 ) 


= (x^-2x/|gr sin | + /rgT)(x^+ 2x/IgT sin + /\^) 
= (x^- 2xq^ sin + q^)(x^+ 2xq^ sin |- + q^j 


Forming a partial -fraction expansion, namely, 

x 2 C^x C 2 X 

~3 ~ ~9 A 7 *9 ft ? 

X +ax +b X -2xqsin^+q x +2xqsin|-+q 

and solving for and C 2 yields 


(31) 


-C 


^ 1 
^ 4q sin |- 


(32) 


Thus, 


^^+ax^+b 4q sin ) jx‘ 


1 ^x^-2xqsin |-+ q^ 


dx 




J A 7 

X +2xq sin 2 + q 


dx 


< 0 


(33) 


The two integrals in (33) can be evaluated using [5, p. 68, Eq. (2.175-1)], 
namely. 



xdx 


+ gx + Y 


^£n(x^ + 6x + Y) -^tan’^ ’ 


A > 0 (34) 
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where 

2 

A ^ 4y - 8^ “ 4q^ -^±2q sin 
* 4q^ cos^ 7 ^ ® 


Finally* combining (33) and (34) gives the desired result 



h^ < 0 (36) 


or, since from (29) 


sin |- = cos ; cos = sin |- ; sine = sina (37) 

then 



h^ < 0 (38) 


Applying (25) and (38) to the specific evaluations of a and b 
in (19) and f, g, h, and q of (20), 
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where 


and 



( 40 ) 


The corresponding definite integrals with limits of zero and Infinity are 
identical and are simply given by 



16 
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Applying (45) to the case of Interest gives 



where, as before, and R“ are defined In (40). Finally, Integrating 
(46) between 0 and [note the upper limit cannot be Infinity since 
the first term of (46) would become unbounded] and Introducing the 
notation 





■r^ 
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Noting that 



(49) 


Jx +ax +b J X +ax + 


The second integral in (50) can be obtained from a generalization of 
(25) and (38). In particular, multiplying (25) by A and using [5, p. 67, 
Eq. (2.161-1)], we get 


Ax^t B_ ^ 1 r(A9 i B l t -1 /JLU tan'V— (]; < 0 (61) 

x^ax^b ,'g Ui /r W 

— ■ 


Similarl>, following the procedure which produced (38), we get 
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Letting A « a and B 
the desired result. 

Since, to evaluate o^' 


b and substituting (51) or (52) Into (50) gives 
21 


1-way’ 


we need the definite Integral 


version of (5) with limits 0 and u>|^, letting a and b be as defined 
In (19), we get, after some simplification 



n 

tan-' /- 



V 

/F Wl 


R > 2 (53a) 


and 



R < 2 (53b) 


Finally combining (50) and (52) with (16) and recognizing the further 
simplifications 


and 


RR~^ - 1 



1 - RR 



(54) 


gives the desired result 
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(56) 


When the loop has c = 0.707 damping (I.e., r = 2 or R = r-2 = 0), then 
(23), (42), (48) and (55) simplify immensely. In particular, 

2 . •'3 

I ~ 2 

1 -way 4 u3^ 


2 ^ ^2 
^ 1-way 
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The ^irst three relations in (57) agree with similar results* 
obtained in [2]. 


4.0 CALCULATION OF MEAN-SQUARE PHASE NOISE FOR 
TWO-WAY COHERENT COMMUNICATION LINKS 

Here we consider the evaluation of mean-square phase noise com- 
ponents whose forms resemble the first two terras of (9). Subsection 4.1 
of this section treats the evaluation of the first of these two compon- 
ents (i.e., the in-band phase noise in the vehicle system) while subsec- 
tion 4.2 considers the second term corresponding to the out-of-band 
vehicle system phase noise. 

For simplicity of notation, we shall use the integers 1 and 2 
to replace, respectively, the subscripts V and R previously used to denote 
vehicle and reference systems. Furth‘;rmore, it is noted that the vehicle 
and reference system second-order carrier-tracking loops have in-band and 
out-of-band transfer functions of the form given by (10) and (11), with 
arbitrary damping actors r^ r 2 and radian natural frequencies w^ 2 * 

4.1 In-Band Mean Square Phase Noise 
Consider the evaluation of 


o 


2 

I 


2-way 





n - h2(jw)i^ 


d(u 


(58) 


2 2 

where |H^(j(j)| and |l-H 2 (ju)| are both given by (12) with the appro- 
priate subscript appendages and S^(u) is defined as before in (14). Sub- 
stituting (14) combined with (1) into (58) yields 


0 


2 

I 


2 -way 


4 


3 



k=0 


2 -way 


(59) 


where 


* 

Eq. (11) of [2] (which includes many typographical errors) 
essentially expresses the sum of the first three results of (57). 
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with 


23 


A « (b2-b^)^ + (a2-«]) («2 " «i M 


Then using (17) and (45), we obtain the desired result, which is 



i = 1.2 (65) 


The parameters f . , g. and h. are related to a. and b. as in (18). 

Ill 1/^1 

Applying (64) to the evaluation of 1 2 -way’ ^i»*st relate 
a., b. , C, and D to the loop damping and natural frequency parameters, 
namely, 

a, = ' 1.2 

'■1 “ ' “n/ ' ‘’1 


V 



24 


Finally, then, ccwiblnlng (61-66) and substituting Into (60)* with k « 3. 
we obtain, after much simplification. 



where we have 1ntroduv.ed the parameter y to characterize the ratio of 
loop natural frequencies. I.e.. 



with weighting coefficients 





★ 

An Infinite upper limit is used here. 


and 
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(70) 


b. (k«2) 

2 

To evaluate 012 I 2 ^ ^ again begin b> considering the 

generic form of the Integrand and expanding It Into Its partial fractions, 
I.e. , 


Cx^ + Dx^ 

(x^ + a^x^ + b^) (x^ + agX^ + b2) 


A^x^ + B^ A2X^ + B2 

~A 2 * ~4 2 

X +a^x +b^ X +a2X +bg 


(71) 


where now 


A 


1 


B 


1 



D|b2-b]| + C|a2bi -aib2| 

_ 

b^|o(a2-ai) - c(b2-bi)] 

B 2 = 


(72) 


and A Is still as in (63). Using (51) and (52) to Integrate the two terms 
on the right side of (71), we immediately get the desired result, namely. 



( 73 ) 
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When (73) Is evaluated between zero and Infinity and the appropriate 
substitutions are made for f^» h^» q^» and In terms of r^ and 
o>p{» i * 1»2» we get the much simplified result 



■'V* \ V / 

1 

all R,, Rg (74) 

Finally, substituting (72) into (74) and relating C and D to the appropri 

O 

ate parameters in the integrand of o \2 |2-way* ®*9** obtain, 

after much simplification. 



. "2 

r 1 A + ^ A 1 


2-way 

/R^+2 /Rg+2 



(75) 


where 


A 


12 ■ 



tR,) • y^(3*rJ Rg ^ (/*2tR,) R,} 

^0 

I (l - (R^ ^ 2 - Y^) - [l ^ y^(r^ 2 )] Rg ^ y^(y^ ^ ^ 1 


and Aq Is still given by (70). 
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c. (k» 1 ) 

For k » 1, we must consider the Integral 



T 


Cx® + Dx^ 




T 


T 


X +a^x +b,J(x +a^x +b, 


dx 


whose Integrand has the same partial-fraction expansion as in (61), 
namely. 


Cx® + Dx^ 


A,x^ + B,x A,x^ + B,x 
' ’ + (77) 


|x^ + a^x^ + b^^^x^ + a2X^ + b2 j x^ + a^x^ + b^ x^ + agX^ + b 


except that now A^, B^, A 2 , and B 2 are given by 


" '^2 




0(^2 -t>i) + C(a2b^ -a^b2) 

A 

b, [ 0(62 - a,) - C(b 2 - b,)] 


(78) 


with A as in (63). Since the right-hand sides of (64) and (65) are now 
also equal to the above integral, then using ( 66 ), evaluation of (60) with 
an infinite upper limit and k = 1 yields, upon simplification. 



2 -way 



+ 


^21 ^^^ 1 ^ ” ^31 ^(^ 2 ^ 


(79) 



Y'‘[l-y"t(R,t2)(R2Y^-R,)‘ 

^0 


where 
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d. {k=0) 

For k = 0, we must consider the integral 


h 


Cx® + Dx^ 


x^ + a^x^ + (x^ + + b2^ 


dx 


whose integrand has tne same partial -fraction expansion as in (71), namely 


Cx® + Dx^ 


(x^ + a^x^ + b^ ^ (x^ + a2X^ + h2^ 


A, x^ + B, A, x^ + B, 

~T 

X +a^x +b^ X +agX +b2 


where, now. 


A^ = C - A 2 = 


- c[b,(b2- b,) +a,(a2bj -a^bjjJ tD^a^b, -ajbj^ 




b,|c(a2b, -a,b2)*D(b2-b,)J 


A 


(82) 
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Since the right-hand side of (74) also applies to the above integral, 
then, using (66), evaluation of (60) with an infinite upper limit and 
k = 0 yields, upon simplification. 



Analogous to what was true for the one-way case, the two-way 
mean-square phase noise expressions for (67), (75), (79) and (83) simplify 
immensely when i-e., both loops have 0.707 damping. Thus, 

letting = R 2 = 0 in (67), (75), (79), and (83) yields 



A simple check on the result in (83) can be had by realizing 
that this calculation was performed in [3, Chap. 3]. In particular. 
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a 


2 


10 


2-wa^y 




( 86 ) 


where W|^^; 1 » 1,2 denotes the two-sided loop bandwidth and K|^ 
Is defined in [3, Eq. (3-18)] with 

« **n (■‘2 ^ W|.i(3*H2) 


(•^1 

(87) 


Since loop bandwidth and natural frequency are related by 


Wli = 


“n1 0 ^ , *^n1 ^ ^0 

2/2 + R7 


Zi/rT 


; i = 1,2 


( 88 ) 


then using (68), we can rewrite (87) as 




Y 


/rJT? 

J «,*2 


Finally, equating (86) and (83) gives 




+ 2 . J R^+z) 



^20 


3 + R, 


(89) 


(90) 


Substituting (84) into the right-hand side of (90) and executing some 
routine algebra results in [3, Eq. (3-18)], thus proving the result 
in (90). 
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Out»Of»Band Mean-Square Phase Noise 

In this section, we consider the evaluation of 


2-way 


= |1 - St(w) 


£ \ I 2 . . 


where once again |1 -H.Caw)! ; i = U2 are both given by (12), with the 
appropriate subscript appendages and S^(w) is defined by (14). Substi- 
tuting (14) combined with (1) in (91) yields 


2 ^ 

“ S ®0k‘‘l 

2-way k*0 1 2-way 


where 


2-wiy 


/ "t ii-H2(wr 


, \ p ^ 


The integrand required to evaluate 0Q3 j 2 -way 
partial -fraction expansion as ( 61 ), namely, 


+ a^x^ + b J + agX^ + 


A,x^ + B,x A«x^ + B«x 
-T — 

X +a^x +b^ X +a2X 


with now 


A , -A = V’ ' . 

n-t -no . « 


B, = ^ 
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Thus, letting D « 0 In (64), we can apply the right-hand side of this 
result to (94) which, when combined with (95) and (66), gives an equation 
Identical to (67), namely. 



Once again for k = 2, we can make use of the partial-fraction 
expansion of (71), namely. 
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Thus, (73) and (74) apply directly to this case which, when combined with 
(99) and (66), give an equation Identical to (75), namely. 





( 100 ) 


where now 

*12 = ^ 
'^2 ■ 


( 101 ) 


c. (k=l) 

For k = 1, the upper limit of (93) cannot be allowed to 
extend to infinity since the integral would become unbounded. Never- 
theless, the partial -fraction expansion of (61) still applies, namely. 


^x^ + a^ x^ + b^ j^x^ + agX^ + bg) 


•j 3 

A^x + B^x A2X + B2X 

~d ? ~4 2 

X +a^x +b^ X +a2X +b2 


( 102 ) 


now with A^ and B^ given by (99). Thus, again letting D = 0 in the right 
right-hand side of (64) and evaluating the result between the limits of 
0 and (»^ gives 
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where f(x; h.^) is defined in (65). Finally, substituting (99) together 
with (66) into (103) and simplifying enables evaluation of (93) for k * 1, 
namely. 



and 
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A 


4 


'll 


|(l-yVR,(r%-Rj 


4r 


^1 






-§■ (l -/) + r^|V- 


31 


/ 4 -R^ 


tan 




. tan-V-^\ 


f2(*^»r) ="< 


12 


■ 1 / 2 r^\ 

\/^) 

J aiiWSY J R_L 

\2r^ + R - /r^-4 / \R - 




/K 


; R > 2 (106) 


Also note that 


lira fp(R,r) = f(R) 

r-H» ^ 


(107) 


where f(R) Is defined in (69). 
d. (k=0) 

Here again, the upper limit of (93) cannot be allowed to 
extend to infinity. In addition, before applying the partial -fraction 
procedure to the integrand required in (93) for k = 0, we must first 
reduce this integrand to one whose numerator is a lower order polynomial 
than its denominator, i.e.. 


X +a^x +b 


1 ) * ® 2 ^^ * 4 ) 


'x^ + Dx^ + Ex^ + F 


+ a < +b 


,)( 


X +a 2 ^ + 62 ^ 


1 


( 108 ) 
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where 

C » d] + 82 

0 ■ + bg + 8^82 

E « 8^b2 + 82b^ 

F = b^b2 (109) 

Applying the p8rtial -fraction expansion technique to the second term in 
(108) yields 



Evaluating (109) using (66) and substituting into (111) results, upon 
simplification, in 



012) 
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Using (51) and (52) on (110) nr^th limits of zero and we 
arrive, after much simplification, to tl .Inal result, namely. 



and R” are given by (40). Also note that 

lim g(R,a,6,r) = j IStsd . 
r-x» ^ /TTr 


all R 


015) 




5.0 
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CALCULATION OF MEAN-SQUARE PHASE NOISE FOR TWO-WAY LINKS 
WITH COHERENT TRANSMIT/RECEIVE OSCILLATOR SOURCES 


Until now, we have implicitly assumed that the receive and 
transmit phase noise sources, and which emanate from the refer- 
ence system TCXO (see Figure 2) can be considered as statistically inde- 
pendent random processes insofar as their effect on the ccmiputation of 
mean-square pi.ase noise in this same system [see (9)}. Since, in real- 
ity, and are derived fnrni the same oscillator source, the above 
assumption is justified by considering tte extreme situation wherein the 
round-trip delay of to the vehicle system and back is long relative 
to the correlation time of the process itself. As such, a sample of 
which leaves the TCXO arrives back at the reference system essentially 
uncorrelated with the corresponding sample of 

The other extreme is to assume the case where the round-trip 


delay time to the vehicle system and back is short relative to the cor- 
relation time of the phase noise process. Here, the simplest model is 


to assume ^TX * 


whereupon (5) is modified to 


♦e(s) = Hy(s) - *TX^^^ 

+ - H^(s)][i - Hr(s)] *^(s) 

- [l - H^(s)][^^(s) + (116) 


and, correspondingly, (9) becanes 


S- U) <1» 

“<X> V 

.» r 1 


(117) 
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Thus, we must now consider evaluation of a new type of integral whose 
generic form [see the first term in [117)] is 



1 



|6H^(jo,)-l|^ ll-HgMI^ S-(u,) du 


(118) 


where, as before, we use the subscripts "1" and “2“ to replace “V" and 
“R", respectively, and, for further notational simplicity, we drop the 
“f“ subscript on the frequency multiplication ratio (turnaround tran- 
sponder ratio), G^. 

From the closed-loop transfer function defined in (10), we 
iimediately find that 


GH^(s) - 



* “nl 

* “nl Al * ^ “nl^ 


(119) 


Thus, from (11) and (119), 


I(s)I(-s) = (gh^(s)-i)(i -H2(s))(gh^(-s)-i)(i -H2(-s)) 

( 1 - G )2 ^2 ^ 




s^-o) .^(r.- 2 ) + « A 

' - '"ni / 


ni ' 1 


8 




+ (1-G) 




( 120 ) 


Finally, letting s = ju 
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( 121 ) 


When (121) is substituted in (118) together with the definition 
of 57 ( 0 )) given in (14), then, as before, we can express the two-way mean- 

y 

square phase noise as a sum of four tenns, i.e.. 




2 -way 


( 122 ) 


where 



(123) 


2 

f'lote that the first term of (123) is identical to ogj^ 2 -way 

(93). A-:.!:.') troiii the (1-G)^ multiplication factor, the second term of 

2 

(123) is sinilar to oji^ 2-way difference being that 
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2 

the coefficient of ui in the numerator of the Integrand Is r^-2/(l-G) 
rather than r^. Because of this similarity, evaluation of the second 
term of (123) should yield results identical to (67), (75), (79), and 
(83), with the only change being the definitions of the coefficients 
In these expressions. In particular, one can show, after much tedious 
algebra, that: 


(k=3) 



(125) 
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mi 

A * ^ 


[Y%-R,)j 


ni 

A * 


^0 

|P*^1 . 2G 

♦ v^("2-y%) - - 



^21 

/{ 

A , = -V 

-(l-T*) 

^0 


''31 

(k=0) 

A = \ 



^0 

026) 

r)Y(y%-R,)[(R,Y^)(l^l!,) -l]j 

no 

1 - v'')(i 


^0 

•^l ^’^ 2 ) ” 

1 

''20 

F 

®GI 

inally, then, 

2 2 
' 2-way ■ > 

"o 

+ 0-6)^ a,,. 
2-way ' " 

027) 

2 

9 

‘ 2 -way 


k=0.1,2,3 (128) 


where the prime on the 'T' subscript denotes the fact that (67), (75), 
(79), and (83) are evaluated using (124), (125), (126), and (127), 
respectively. 
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6.0 SYSTEM MODEL WITH FREQUENCY SYNTHESIZER 

In this section, we generalize the simple two-way coherent 
comnunl cation system discussed In Section 2.0 and Illustrated In Fig- 
ure 1 to Include the other Internal frequency synthesizer oscillator 
sources and other TCXO paths. In particular, we now consider the two-way 
system Illustrated In Figure 4 whose equivalent linear baseband model Is 
given In Figure 5. As such. Figure 3 also applies as the block diagram 
representation of Figure 4 and, thus, the reference system phase error. 

Is still given by (2). However, and the equivalent phase 
noises at the receive and transmit outputs of the frequency syntheslzei , 
must now be expressed In terms of the Individual oscillator component 
phase noise sources within the synthesizer Itself. To do this, we must 
postulate a specific synthesizer structure. 

Figure 6 illustrates the frequency synthesizer used in the 
Payload Interrogator (PI) of the Shuttle orbiter, and Figure 7 is its 
equivalent linear phase noise model. Using routine control loop analy- 
sis methods, it is straightforward to show that 


and 


*^TX ^ ^ ^ 


where 


*TCX0^^^ 

+ [i - *1^^^ 


= 2 


® '*■ DE ^RX^^^ ^TX^^^ " ^ E ^TX^^^ ‘ ^RX^^|]|'^TCX0^^^ 

r '•'1^^^ * ®v ’*'VCX0^^^ 

[l - 


(130) 


(s) 

H,,(s) N 




Htx(s) 


(s) 

4 N 


+ ^F 

^ ai r 1 


^ n" ''tx^^^ 


(131) 




Figure 4. A Simple Block Diagram of a Two-Way Coherent Communication System Employing a Frequency Synthesizer 


4 ^ 





Vehicle System 


Reference System 


Figure 5. An Equivalent Linear Baseband Model of the System Including Oscillator Phase Noise Sources 
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►Rx Output 



Command 

Subcarrier 
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Figure 6. PI Frequency Synthesizer 










Figure 7. Phase Noise Model for PI Frequency Synthesizer 



Also, 
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6 for L-band 
10 for S-band 


(132) 


Note from (129) and (130) that the TCXu phase noise <l>-p^xo trans- 

mit VCO phase noise 4*i each contribute to both the transmit and <^ece1ve 
synthesizer outputs. Furthermore, the first terms of the TCXO contribu- 
tions, namely, 16 4 *jqx 0* direct feedthrough paths of the 

TCXO previously accounted for In Section 2.0. 

Substituting (129) and (130) Into (2) and making the notation 
simplifications and asijmptlons leading to (5) gives the desired result 
for the reference system (PI) phase error in terms of all the Individual 
phase noise sources, namely. 


4>Js) = 


Sf Hy(s) jj - H|^(s)^ 2 ^ 

- 8 ^ Htj((s) [l - 

+ H^(s) [l - H,^(s)][l - H^j((s)] 4 . 2 ( 5 ) 

+ Gf [i - H^(s)][l - H^(s)] 4 iv(s) 

[1 ' Hj^(s)J [^6 + 2(^+4) Hj^j^(s)J 4 'tcxO^^^ 

\] ■ ■ ^RX^^j] 

- - Hp(s)J j^f^(s) + 4.^(s)J 


( 133 ) 
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The frequency division constants E» F» N and 0 of Figure 6 
depend on the operating mode of the vehicle transponder (lUS) and the 
channel number, CN, assigned to the particular frequency selected In the 
frequency band corresponding to this mode. In particular, the following 
table applies: 


Frequency 

Division 

Multiple 

STDN 

DSN 

SGLS 

E 

240 

240 

256 

F 

222 

221 

205 

N 

CN + 519 

CN - 431 

40CN - 35460 

0 

1280 

432 

1280 

CN 

1-808 

850 - 882 

900-919 

<5f 

240/222 

240/221 

256/205 


From the above, we observe that, for all three modes, 

Gf (I) = 1 (134) 


Also, from (133), the worst-case (largest) phase error would occur for 
the maximum value of N which corresponds to the highest channel number. 
Thus, 



max 


/1327\ _ 
\1280/ “ 

1.0367 

STDN 

(m\ = 

1 432/ 

1.044 

DSN 

/1300\ _ 
\1280/ 

1.0156 

SGLS 


(135) 


Before proceeding with the evaluation of mean-square phase 
noise corresponding to >l>g(s) of (133), we point out that we have inten- 
tionally written the pair of contributions from 'f'|(s) as 

separate terms in this equation. The reason for this is to identify the 
fact that in each case one of the contributions (the first of the pair) 
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originates at the transmit output of the frequency synthesizer, travels 
through the vehicle, and returns to the reference system phase error out- 
put, while the second contribution of each pair enters the reference sys- 
tem directly frcyn the receive output of the synthesizer without ever 
passing thi'^ugh the vehicle system. As was done In Section 5.0, we shall 
assume In what follows that the round-trip to the vehicle system and back 
Is long relative to the correlation time of either phase noise process 
Itself; thus, the two phase noise components In each pair are essentially 
uncorrelated. 

We now turn to deriving expressions for the necessary 3-way and 
4-way mean-square phase noise components associated with the phase error 
process In (133). Furthermore, to simplify matters. In this section, we 
shall Ignore altogether the terms corresponding to the TCXO phase noise 
since, typically, these contributions are much less dominant than those 
arising from the other phase noise sources. Thus, analogous to (9), we 
get 



• CD 



7.0 
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CALCULATION OF THREE-WAY AND FOUR-WAY 
MEAN-SQUARE PHASE NOISE COMPONENTS 

Here we consider the evaluation of mean-square phase noise 
components whose forms res«nb1e the first three terms of (136). Sub- 
section 7.1 treats the evaluation of the third term, while subsections 
7.2 and 7.3 consider the second and first terms, respectively. 

As was done in Section 4.0. we shall, for simplicity of nota- 
tion, use the integers 1 and 2 to replace, respectively, the subscripts 
V and R previously used to denote vehicle and reference systems. Fur- 
thenw)re, we shall use the subscripts 3 and 4 to replace, respectively, 
the subscripts TX and RX previously used to denote transmit and receive 
synthesizer output terminals. Also, as before, all of the second-order 
carrier-tracking loops (vehicle system, reference system, transmit syn- 
thesizer, receive synthesizer) have in-band and out-of-band transfer 
functions of the form given by (10) and (11), with arbitrary damping 
factors r. ; i=l ,2,3,4 and radian natural frequencies i*l ,2,3,4. 

7 . 1 Out-of-Band* Three-Way Mean-Square Phase Noise 

Consider the evaluation of 


3-way 




= 27 1 iH^(jo.)r 11 -H2(j«.))r ii-H3(ju)r 


S^(w) 


du> (137) 


2 2 

where !H^(ju>)| and ll-H^(jw)! ; 1=2,3 are given by (12), with the appro- 
priate subscript appendanges and S^(u) is defined as before in (14). 
Substituting (14) combined with (1) into (137) yields 


2i A 

3-way 


3 

V 


‘^Qk^ 


3 -way 


where 


* 


The term "out-of-band" 
third tracking loop. 


( 138 ) 


is used here with respect to the 
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(139) 


a. (k=3) 

To evaluate °q 3^|2 wa y* clear from (60) that we iraist first 

consider the evaluation of an Integral of the form 



Cx^+Dx® 


x^ + a.x^ 



dx 


where 


a. 

1 



^<0. 


m 


R. 

1 




; 1 = 1,2,3 


C = 



a,„,^(R,+2) ; 0 = (140) 


The procedure employed is, as before, to expand the integrand into its 
partial fractions, namely. 


_ 7 , r, 5 
Cx + Dx 


3 

n 

i=l 


TT + b 




^ A ^ + B ^ 

y 13 *^13 

^4 2 

i=l X + a.x + b^ 


(141) 


Then, integrating the right-hand side of (141) using (17) and (45) gives 
the desired result, namely. 
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Cx^ + Ox® 


4 2 

X +a^x‘ 


♦ 



dx * 



042) 


O \ 

where f^x; h^'^^ is given by (65). Before evaluating (142) between the 
limits of 0 and infinity, we digress to discuss the determination of the 
six unknown coefficients 8 ^ 3 ^ i*1.2,3. , 


the relation 


Multiplying both sides of (141) by TT (x^ + a^x^ + bi) gives 

i=l ’ 


5 (a^^x- 


M3' 




i=l 


) TT (x^ 

' £=1 

m 


+ a^x^ + bjj) 


= Cx^ + Dx^ 


(143) 


Equating coefficients of like powers of x on both sides of (143) results 
in a set of six linear equations whose solution gives the desired unknown 
coefficients. 

Rather than obtain specific ^losed-form expressions for these 
coefficients aj was Jone previously in the one-way and two-way results, 
it is more convenient here to describe the solution in matrix notation. 
Specifically 3 t denote the unknown coefficient vector 




'13 

^3 

^23 

^23 

S 3 

S 3 


(144) 


and Fj the forcing function vector 
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(145) 


then the solution to (143) as described above can be written as 


[NlVj = Fj 


(146) 


or, equivalently. 


*3 ' t"] f3 


where, makiij use of (140), the 6x6 matrix [M] is characterized by 


[M] 


1 2 • • • • 1(1-^ g 

^21 ^22 • • • • ^26 


"*61 '^’62 * • • * *"66 


(148) 


with 
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» 1 ; i«l,3,5 



iDj.* “ j~2»3...j6 (149) 

Using the first two relations of (149) in (147), i.e., the ele- 
ments of the first row of [M], we immediately see that 

3 

A-3 = 0 (150) 

i=l 

Furthermore, from the second row of [M], we have that 

3 r n ^ 

I ®i3-(“n(S)^3 = 1 “ ‘'5'> 
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or, equivalently. 



(152) 


Solving for from (150) and (152) and 

substituting into (142) gives the alternate form 



It is now convenient to evaluate (153) between the limits 0 and infinity 
which, upon algebraic simplification, results in 



where f(R) is defined in (69) and the parameter characterizes the 
ratio of loop natural frequencies relative to that of the third loop, 
i .e. , 
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4 “ni 




'n3 


i'1.2 


055) 


Finally, then. 



b. (k=2) 


To evaluate a 


Q2 


3-way 


we again begin by considering the generic 


form of the integrand and expanding it into its partial fractions, i.e.. 


Cx® r Ox® 


TT(x 

i=1 


^ + a^x^ + b 


,) 


y ^2^ »Bj2 
i=l + + 


(157) 


where a^, b^. ; i=l,2,3 and C, D are as previously given in (140). From 
(71) and (73), the integral of (157) is easily seen to be 



Cx® Dx 
4 

x + a^x 


6 


? 



1 


1=1 


f (x;A^2»5i2) 


dx 


(158) 
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where f(x; A^,B^ )is given in (73). When (158) Is evaluated between 0 
and infinity and the appropriate substitutions are made for f^, h^, 
q^, and a. in terms of r^ and i*l,2,3, we get the much simplified 
result 



Cx® + Dx® 

TT (x^ + a.x^ + b.) 
i=l ^ ^ 


dx 



(159) 


Finally, then, 



(160) 


Evaluation of the coefficients ^^2’ '*®1»2,3 proceeds as before 
by multiplying both sides of (157) by 



+ a^x 


+ b 




and eauating ccefi'icients of tike powers of x. The resulting set of 
linear equations has the solution 

V 2 = [2 (161) 


where now 
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and [M] is still described by (148) and (149). 

In the interest of brevity, we shall merely present the results for 
the cases k*l and k=0. With only minor changes, the procedure to obtain 
these results follows that just described for k=3 and k=2, respectively. 
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where 




(165) 


with 




1 

c 

1 


»„,^(R,*2) 

®1k 


4 

“nl 

*2k 

“2k 

• F = F • F = 

’ -1 L2 * -0 

0 

0 

*3k 


0 

_ ®3k _ 


0 


(166) 


7.2 I n-Band Three-Way Mean-Square Phase Noise 

In this section we consider the evaluation of 


3-way 


^ |l-H2(ju.)!^ lH3(ju.)l^ S^(.^) du. 


(167) 


,2 .. 


where once again [1 -H^(j(.i)| ; i=l,2 are both given by (12), with the 
appropt i«te subscript appendages and S-(a)) is defined by (14). Substi- 
tuting (14) combined with (1) in (167) yields 


3 -way 


3 

V 

z- 

k=0 


V 2 
1 ‘'Ik 


(168) 


3 -way 


where 
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C = “n/ '"i "n3^ '■3 ' 

“ ' »nl^'-I*“n3^''3 ' * “n3^(“3 * 

E -mV <’'“> 


Since the partial fraction approach taken in subsection 7.1 
also applies to the above integral, we can essentially write down the 
desired results by inspection. Thus, rather than go through the develop- 
ment for each value of k, as was done previously, we shall merely present 
a summary of the results in compact notation. Thus, 
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For k«1 .3 



where i=l»2 is defined in (155) and the unknown coefficent vector 
defined in (166) has a solution of the form in (165) with, however, 
forcing vectors given by 



(173) 


The parameters C, D and E in (173) are defined in (170), thus completing 
the solution. 
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7.3 Four-Way Mean-Square Phase Noise 

In this final subsection, we consider the evaluation of the 
four-way mean-square phase noise component necessary to compute the first 
tenn of (136), namely. 


4 -way 




lH^(joi)|^ |1-H2(jw)|^ |H3(jo.)|2 Il-H^(jo))i^ S-(u.) do. 


k=0 


4 -way 


(174) 


where 


. ..2 


X.2 


4-way * 

D 


II 

U T |l-H 2 (j(o)r |H 3 (ja.)r |l-H 4 (ja.)l^ da, 




21 2 ^ 2 \ 21 2 ^^ 2 \ 
‘"nl ("’nl ^’'l^ )‘"n3 (“n3 ^*^3“ ) 


2 . 2\ 8-k 

0, 


4 . 2T 2 


■= — do. 


(175) 


TT w + w M .^[r. - 2) + 0, 

[jni 1 ^ m I 


Comparing (169) with (175), one wou^d intuitively expect results 
which are similar in form to those of (171) and (172). Indeed, this is 
the case. In fact, without belaboring the details, one obtains 

For k=l ,3 



(176) 
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For k»0.2 


4-way 



W!V'"e the are now defined relative to i.e.. 


A ni 




n4 


1=1 .2.3 


(177) 


(178) 


The unknown coefficient vector which now has eight components, I.e., 


(179) 


still has a solution of the form in (165), where the forcing vectors are 
the same as those in (173) with the addition of two zeros to the bottom 
of each and the 8x8 matrix [M] is characterized by 


h 


'Ik 

hk 

'2k 

*2k 

'sk 

^3k 

'4k 

^4k 


[M] 


m^l m ^2 ' • • "’]8 

ru^i ^^^22 • * • ^28 


'"81 '"82 '"88 


(180) 
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with 


■"li 

"’ll 


m^i 


"*31 


"«4i 




5i 


"•ei 


1 ; 1»1 ,3.5.7 

0 ; i=2,4.6.8 
4 

y w .^Ri 

nj J 

j=l 

j^i 

4 4 


I .S “nj Vnk\ * 2 “„j’ 

^ j I, Ji 


j.k^i 

m 

4 4 


j=i 


J= 
j.k}«i 
J7k 


y y t. + 1 T /r- 

^ A. nj j nk .'J, no o 

i=i k=i j-j 


4 4 


4 4 4 


° 1 S S “nj\k^ ♦ II S S “njS“nk^“ 


■ <•> s 

k n£ 


>i=l, 3,5,7 


j=l K-| 
j.kj^i 
jj'k 

4 4 4 


j=] k=l £=1 
j, k, £^i 

OJ‘kj‘£ 


III “„j^v„k\t' 


j=l k=1 £=1 
j,k,£^i 
j^k^£ 


m. 


7i 


m,. 

ol 


"•ji 


= TT 
j=i 
J7i 

= 0 


<ij 


nj 


m 


i=2,4,6,8; j=2,3,...8 


(181) 
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APPENDIX VII 


MEAN-SQUARE ERROR APPROXIMATION TO PHASE NOISE SPECTRA 


Here we wish to derive a formula for approximating a given 

-ic 

phase noise power spectral density S(f) by a power series in f , keN, 
in the least mean-square error sense (LMSE). 

S(f) is either directly available from measurements or is given 
in a piecewise linear form. The approximating function S(f)is assumed 
to be 


S(f) = h(j + h^f^ + h2f + h3f"^ (1) 

where h.; i = 0,1 ,2,3 are unknown coefficients to be determined. 

* ^ 

The specific form of S(f) in (1) is dictated by physical con- 
siderations. The f^ term corresponds to white phase noise, f~^ corre- 

_2 

sponds to white flicker noise, f corresponds to vdiite frequency noise, 

_3 

and f serves to model flicker frequency noise. The physical signifi- 
cance of (1) then is that it represents the spectral superposition of 
the n., rementioned random components, assumed independent, in the fre- 
quency and of interest. 

The spectra to be approximated usually exhibit a flat portion 
for high frequencies. We can then assume that this value, which is 
directly measurable, determines h^. Henceforth, we will deal with the 
evaluation of hj^, k > 1. 

The frequency band of interest will be denoted by (a, 6). We 
shall also assume, as is coiranon practice, that g » a. Indeed, custom- 
arily, a is on the order of tens of hertz while g is on the order of 
mega' i tz. 

Since hg is known, let us define 

S'(f) = S(f) - hg (2) 


The function to be minimized is 



2 



Since E is a quadratic function of h^; n = 1,2,3, it follows that 

lim E = +«> 

n=l,2,3 


Therefore, the extranum of E, to be determined from the condition 


9E 


= 0 . 


n = 1,2,3 


» 


clearly corresponds to a minimum. 

i!}' fferentiating (3), we get 


or 


Cl 



0 




df ; n=l,2,3 (4) 


0 ,.,- ;■ 1 1 , • ‘hat, for n > 1 , 



LI _L 

n-1 ,n-1 

I? P 



3 


and inserting the assumption « a"\ we eventually derive from (4) 
the following system of equations: 


h, + -iir ho + h, = W, 


a "1 " ^ "2 " ^ "3 


^5 '>1 * ^ "2 ^ ^ *-5 = «3 


where we have defined 


6 c. 


W ^ i n=U2.3 

" \ f" 


(5) 


If we further define 


h’. 


. h 
4 JL 
n n 


W. 


n=l,2,3 


4 n-1 ,, 

- a W„ 


( 6 ) 


the system I is written in matrix notation as 
1 1/2 


1/2 1/3 

1/3 1/4 


1 



=r 

1 


1 

1 


h*2 

II 

W*2 

i 

ro 

! 


_"‘3_ 


II 


The solution of II is of the form 


h* 


ji 

A 


(7) 
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where 


A = 


end 




W ' U ' U ' 

•1^ viQ 

240 ■ W*7T 

(8a) 

-wr 4 W 2 ' W 3 - 

“Sir ^ 45 • 

(8b) 

72 ■ 12 12 

(8c) 


From (6), (7) and (8), it eventually follows that 


a|g W^' - 36 W^' + 30 W 3 * 

12 

- 3 W^‘ + 15 W^' - 1 

30 

■ - 6 W 2 ' + 6 W 3 ' 


(9) 


where, from (5) and ( 6 ), 


W ' = 




df; n=l,2,3 


( 10 ) 


Once the W ' is evaluated as in (10), then the h follows direct »y from 
n n 

(9). However, except for the unlikely case that S(f) is given in a 
closed form, evaluation of (10) is not apparent. To circumvent this, we 
propose the following: assume that S(f) is given by a dB-amplitude/log- 

frequency plot, as measured directly in the laboratory. We approximate 
the plot by a series of N piecewise linear segments, each valid in the 
region f f 5 ^i+i’ i = U---tN, where F-j = a and = g. Let us also 
denote by the value of the spectrum at F^ , i.e., S- = S(F^). Obviously, 
the quality of approximation to S(f) increases as N increases. The 
piecewise linear approximation in the log/log scale implies that S(f) is 
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assumed to have the analytic representation 


/f \ 

s(f) ' Sj(f) ' 


for 


^i 1 ^ 1 


In (11), the exponent -P^ represents the slope of the line, viz. 


p slope of line in dB/decade 
“ 10 


Direct substitution of (11) into (10) gives 



( 11 ) 


( 12 ) 


(13a) 


(13b) 


and 
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Note In the above that. In general, the P^'s need not be Integers. 


Comment 

The above theory solely derives the LMSE estimate S(f) of S(f). 
However, no care has been taken to account for the fact that S(f) should 
be nonnegative as it represents a spectral density. Indeed, on some 
occasions, the results provide a negative S(f) for some frequency regions. 
In these cases, the coefficients of the LMSE estimate can serve as good 
starting values which, after some heuristic small perturbation, produce 

A A 

a good new estimate S*(f). This new S*(f) should also conform with the 
aforementioned physical constraint of positiveness. The above comments 
were used in Example 2 following. 

Applications 

We now apply the previous results to the following two 

examples. 

Example 1 

Consider the TCXO phase noise spectral density of Figure 1. 

The coefficients of S(f) were calculated to be 



= 0.31622776 ^ 

10’^^ 


= 0.39822958 x 

O 

1 



-7 

h- 

= -0.14756626 x 

10 ' 

^7 

= 0.31740942 x 

O 

1 
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Table 1 following compares certain values of 
corresponding ones of S^f). 


(f) with the 


Table 1 


Frequency 

f(Hz) 

10 Hz 

100 Hz 

1 kHz 

10 kHz 

100 kHz 

1 MHz 

(db) 

-75 

-105 

-125 

-135 

-145 

-155 

m 

-75 

-104.7 

-123.8 

-134 

-144 

-155 


Example 2 

Next we consider the VCXO phase noise spectral density, shown 
in Figure 2. It was approximated by three piecewise linear segments. 
Table 2 summarizes the characteristic values considered for the imple- 
mentation of the previous algorithms. 


Table 2 


n 

1 

2 

3 

4 

F„ (Hz) 

1 Hz 

1 kHz 

8.1 kHz 

1 MHz 

s„ (db) 

-14 

-104 

-115 

- 


-3 

-1.4 

0 

- 
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The coefficients were found to be 

hjj * 0.31622777 x 10*^^ 

h^ = -0.3000000 X 10'^ 

hg = 0.10874600 x 10*^ 

h3 = 0.39684600 x lo'^ 

Table 3 provides a comparison between S(f) and S(f) at distinct 
points for this case. 


Table 3 


Frequency 

S(H2) 

1 Hz 

10 Hz 

100 Hz 

1 kHz 

10 kHz 

100 kHz 

1 

1 

-14 

-44 

-74 

-104 

-115 

-115 

-115 

s(f) 

(db) 

-14 

-43.9 

-73 

- 99.2 

-118 

-115.4 

-115 


paECECireG .^AGE f^OT FIuMCO 
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CALCUUTIONS AND FINAL RESULTS ON THE PHASE NOISE ERROR VARIANCE 
FOR THE S-BAND PI AND TRANSPONDER FOR SOLS AND STDN 


The noncoherent result is 


*g(s) * 128 (1-Hj^5)(l-Hpj) + ,^^(l-H,)(l-Hpj).H2.Hjy5 

- ♦i'(l-H,)(l-Hpj) + (l-H^) Hjyj (1-Hpj) 

- 10 4»pj (1-Hpj) - 0^^ n-Hpj) (1) 


where ij*^' and 4>i are two Independent processes with the same statistics. 
Furthermore, the spectra o^ ^jus* ^PI* ®^2 (VCXO spectrum) 
as given in Figure 1. Since the spectrum in Figure 1, Sj (f), is mea- 
sured with respect to a 2211 MHz carrier, the following is true: 


(a) 

(b) 


(c) 


n28ij» 


I US 


(f) 


Sj(f) 


i|»pj is measured at 1.84 MHz, hence 

^I0ii»pj^^^ " (Mtt) ^I^^^ 

9 is measured at 31 MHz, hence, 
•^2 

c. 


The common spectrum Sjj(f) (see Figure 2) of the processes ' , i >2 

is measured at 2211 MHz. The method of extracting the shape of Sjj(f) 
from the TCXO output spectrum measurement is described elsewhere (see 
Addendum). Incorporating the above conments in (1) yields 

S^(f) = Sj(f) |2.658xl0''^ il pi(f)l^ + n-Hpi(f)l^| 

+ S„(f)||l-H^(f)|^ |H2(f)|^ il‘Hpj(f)|^+ |l-H^(f)l^ 

X |l-Hpj(f)|^ + 6^ ll-H2(f)|^ M-Hpi(f)l^| 

( 2 ) 



dBc/Hz 
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He now make the following simplification. For the SGLS mode* the 
bandwidth of Hjy^(f) is 3774 Hz, while that of H 2 (f) Is 560 Hz, or about 
1/7 of the latter. Hence, we can safely ignore Hjy^(f) In the fourth* 
order products of (2) and use the approximation 

The above substitution simplifies matters significantly because we now 
have to calculate two three-product Integrals of the form 


and 


f 

f 


r''|H2(f) (1-H,(f))(l-Hpj(f))|^ df 


(3a) 


f"‘'|Hiys(n (1-H2(f))(l-Hpj(f))i^ df (k=0,l,2.3) (3b) 


instead of a three-product and i o*iir-product Integral. 
The evaluation of the variance 


a 


2 


f 


Sjf) df 




requires approximating the spectra Sj(f) and Sjj(f) by a series of the form 


I 


i=0 




-1 


The least square approximation method is described in Appendix VII of 
this report and gives the following results 


Sj(f) « 0.158x10 


■11 0.15x10 


-7 


0.544x10 


0.1984x10 


-1 


(4a) 
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and 


Sjj(f) 


0.5 


(4b) 


Note: The spectra show) In Figures 1 and 2 are one*s1ded, while Sj(f) 
and Sjj(f) are two-sided expressions. Therefore, 3 dB should be sub- 
tracted first frcxn the figure plots before approximation Is attempted. 
The evaluations of the one-product Integral 



- Hpj(f)|2 df 


and the two two-product Integrals 





|1 - Hpj(f)|2 df 


and 



n - 


M - Hp,(f)|^ df; 


k=0.1,2.3 


are performed by the program PHNOISE. Next we describe the evaluation of 
and I 2 . 


1.0 EVALUATION OF and I 2 

According to the theory developed, the evaluation of and I 2 

requires the prior computation of the matrix M, as given by (148) and 
(149) of Appendix VI of this report. This was done by the program M, and 
the results are as follows: 
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I 


M * 


and 



M = 


1 I 0 I 1 

0.4939x10^ I 1 10.327x10® 

0.250x10^^ 10.4939x10^ |0.625xl0^^ 
0. 3923x10^® |0.250xl0^^ |0.3092xl0^® 
0.1814x10^® 10.3923x10^® jo. 983x10^® 

0 10.1814x10^®! 0 


1 I 0 I 1 I 
0 . 3367 x 10 ® 1 1 jo. 307 x 10 ® j 

0 . 1086 x 10 ^® 1 0 . 3367 x 10 ® | 0 . 203 xl 0 ^®| 

0.558x10^^ jo. lOSCxIO^® 10.623x10^^1 

0 . 178 x 10 ^® 10 . 538 x 10 ^^ i 0 . 199 xl 0 ^®| 
0 | 0 . 178 xl 0 ^® I 0 I 


1 0 

1 1 

1 0 

1 1 

jo. 356x10® 

1 1 

10.327x10® 

10.180x10^^ 

jo. 356x10® 

|0.625xl0^^ 

|0.628xl0^^ 

jo. 180 x 10 ®^ 

10.3092x10^® 

'jo. 178x10®^ 

j0.628xl0®^ 

|0.983xl0^® 

1 0 

jo. 178x10®^ 


0 

1 1 

1 0 

1 

jo. 336 x 10 ® 

1 1 

0 . 307 x 10 ® 

jo. 22 x 10 ^® 

10 . 336 x 10 ® 

0 . 203 x 10 ^® 

jo. 739 x 10 ^^ 

jo. 22 x 10 ^® 

0 . 623 x 10 ^^ 

jo. 362 x 10 ®^ 

10 . 739 x 10 ®^ 

0 . 199 x 10 ^® 

1 0 

10 . 362 x 10 ®^ 


Since we are interested in the k=2 term only (see ( 4b ) in connection with 
( 2 )} , we have to evaluate the following column vectors: 



and 


0 0 . 313 x 10 ® I 0 . 983 x 10 


11 


I 0 I 0 I o]^ 


‘2- 


0 I 0 . 28486 x 10 ® I 0 . 202865 x 10 


.15 


I 0 I 0 I oj^ 


The solution of the system [M].[V 2 ] = [F 2 ] is accomplished 
through the subroutine SLEQ for the solution of a system of linear equa- 
tions. Despite the fact that the matrix M is ill-conditioned (i.e., 
the eigenvalues differ by many orders of magnitude), the program produces 
results of satisfactory accuracy. They are 
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Applying (160) of Appendix VI for the proper transfer functions and k=2, 
we find that = 0.11504x10 ^ and I 2 - 0.3034x10”^. These numbers rep- 
resent the integrals and I 2 as defined by (3), not the parameter 
o2 2~\nay Appendix VI. To evaluate the latter, one should 

multiply with h 2 = 0.5 as given by (4b). 

The transfer function parameters (natural frequency f^, damping fac- 
tor 5 ) which were used in the computation of the integrals are given in 
the following table: 



C 


— = 0.707 


0.82 


1.0 


5.0 
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2.0 SUMMARY OF RESULTS AND COMMENTS 

With respect to the calculation of the one-product and two-product 

Integrals appearing In (2 ), we should mention that some of them, notably, 

2 2 

the parameters oq and for both of them, are evaluated with the upper 
limit of Integration being some finite f^ Instead of -hb. This parameter 
was chosen to be f^^ - 6 MHz (I.e., 12 MHz two-sided bandwidth) since that 
was found as the lowest cutoff frequency of the IF filters. Also, the 
gain Is of the value 256/204 « 1.2488. 

The final deviation found Is a. = 8.93® RMS (S6LS mode). In 
the STDN mode, the bandwidth of H^y^ Is 755 Hz, I.e., of the same order 
as the H 2 bandwidth. In this case, the approximation of the four- 
product integral by a three-product one Is not valid anymore, and precise 
results can be obtained only by using the appropriate formulas. Never- 
theless, It Is straightforward to see that the STDN variance is going to 
be slightly less than the SGLS variance found because I^jy^l^l; we vir- 
tually substituted H^y^ with Its upper bound in the previous calculations. 

We can therefore claim that the 8.93® found represent an upper bound for 
both SGLS and STDN modes. 

A final comment pertains to the value of a^. By straightforward 

9 

Integration of the phase noise spectrum of the TCXO output (see appropriate 
memo), one gets a deviation of approximately 3.5®. Adding the same amount 
at the PI receiver demodulator plus additional contributions from the PI 
and lUS VCXO's (of considerably lower value, however, see Figurel), one 
concludes that the value found indeed lies within the expected range. 
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ADDENDUM 


Let us denote by Q(f) the phase noise spectrum measured 
at the output of the PI TCX0» as plotted in Figure 3. From the PI dia- 
gram* It follows that 


hcxo.ot'* ■ lH2(f)l^ 

♦ S (f) |1 Gf + 256 Sg^(f) (6) 

We now entail the reasonable assumption that Sej(f) = 0, i.e., ivS con- 
tribution is negligible. This assumption is verified if we plot the 
product 256 Sej(f) against Q{f) and realize that it is several 
decibels below the latter. Furthermore, we assume that the noise pro- 
cesses ip^(t) and possess the same spectral characteristics. 

Denoting this common spectral density by Sjj{f), we get, from 

(5), that 


c /rx - ^TCXO.O^^^ 

.. ,2 i2 . M ...2.2 


n-H^r IH2T + h-H^r v 


( 6 ) 


Here, 


|H; 


f^ + R. f ^ + f ^ 
2 


n-H,r 




; i=l,2 


where equals 


2 


2 and = 4^^ 


2 



100 Hr 


1 kHz 10 kHz 100 kHz 1 MHz 10 MHz 

Frequency 

Figure 3. Measured Transmitter Single-Sideband Phase Noise Spectral Density 


A I 
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Substituting the previous expressions in the denominator D 
of (6), we get 



Through (6), (4) and the Q(f) of Figure 3, the spectrum Sjj(f) of 

Figure 2 was derived. 



